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READER'S  REFERENCE 


This  report  is  Volume  III  of  IDA  Study  S-341  "Determination  of 
Winds  and  Other  Atmospheric  Parameters  by  Satellite  Techniques." 

The  volumes  describing  the  study  are  the  following: ' 

Volume  I  Summary 

Volume  II  Potential  Needs  for  Determination  of  Atmospheric 

Parameters 

Volume  III  Techniques  for  Determining  Winds  and  Other 
Atmospheric  Parameters 
Volume  IV  Physics  of  the  Atmosphere 
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ABSTRACT 


Techniques  which  have  potential  use  for  explicit  determination 
•of  atmospheric  parameters  by  observations  from  satellites  are  eval¬ 
uated  for  inherent  accuracy  of  measurement.  The  discussion  describes 
the  instrumental  techniques  for  measuring  densities  of  air  and  its 
constituents  by  atmospheric  emission,  atmospheric  absorption  of  solar 
radiation,  deceleration  of  low  altitude  satellites,  occultation  satel¬ 
lite  to  satellite  microwave  transmissions,  backscatter  by  dust  of 
pulsed  laser  signals  and  wind  by  the  motion  of  clouds  determined  by 
time  sequenced  pictures  from  satellites  and  the  displacement  of  bal¬ 
loons  tracked  in  position  by  satellites.  The  principle  of  operation, 
uncertainty  of  measurement,  coverage  in  volume,  and  status  of  develop¬ 
ment  of  each  of  techniques  considered  are  detailed. 

In  addition  to  the  methods  for  explicit  measurement  of  atmospheric 
parameters,  an  implicit  method  of  determining  winds  at  all  altitudes 
by  machine  computation  of  the  equations  of  continuity  of  momentum, 
mass  and  energy,  using  satellite  derived  inputs  of  densities,  tempera¬ 
ture  and  heating,  is  suggested.  Experience  in  application  of  the 
method  for  purposes  of  numerical  prediction  of  weather  for  more  than 
a  day  ahead  is  described.  Special  problems  in  extrapolating  the  ex¬ 
perience  in  tropospheric  simulation  for  the  determination  of  winds 
of  the  upper  atmosphere,  for  which  input  data  are  potentially  avail¬ 
able  by  measurements  from  satellites,  are  named.  Estimates  are  given 
indicating  that  the  upper  atmosphere  simulation  is  within  the  expected 
capacity  of  next-generation  computers  such  as  ILLIAC  IV. 
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INTRODUCTION 


The  primary  objective  of  the  Institute  for  Defense  Analyses  Study 
S-341,  "Determination  of  Winds  and  Other  Atmospheric  Parameters  from 
Satellites,"  is  the  review  of  various  instrumentation  techniques  oper¬ 
ated  from  satellites  and  the  evaluation  of  their  potential.  This  re¬ 
port  undertakes  to  evaluate  the  most  promising  of  these  techniques  for 
the  determination  of  winds  and  other  atmospheric  parameters. 

The  nature  of  the  atmosphere,  which  is  the  medium  to  be  examined 
by  the  measurement  techniques,  is  complice  .  i.  Its  complexities  and 
variabilities  have  been  described  in  "Physics  of  the  Atmosphere," 

Volume  IV  of  this  study.  Requirements  for  knowledge  of  atmospheric 
parameters  to  which  a  number  of  Department  of  Defense  potential  mis¬ 
sions  are  sensitive  are  described  in  "Potential  Needs  for  the  Deter¬ 
mination  of  Atmospheric  Parameters,"  Volume  II  of  the  study.  These 
requirements  of  the  potential  missions  are  the  standards  for  per¬ 
formance  of  the  instrumentation  schemes  described  herein. 

Techniques  which  have  potential  use  for  explicit  determination  of 
atmospheric  parameters  from,  satellites  are  evaluated  for  inherent  ac¬ 
curacy  of  measurement  in  this  volume.  In  particular,  the  discussion 
describes  the  instrumentation  techniques  and  the  characteristic  un¬ 
certainties  of  measurement  for  the  determination  of  atmospheric  param¬ 
eters  by  inversion  of  atmospheric  emission,  by  the  atmospheric  absorb- 
tion  of  solar  radiation,  by  the  deceleration  of  low-altitude  satellites, 
by  the  occultation  of  satellite-to-satellite  transmissions,  by  the 
measurement  of  backscatter  by  dust  of  pulsed  laser  signals,  by  the 
motion  of  clouds  determined  from  TV-type  (APT)  pictures  from  satellites, 
and  by  the  motion  of  balloons  tracked  in  pcs? tier  by  satellites. 
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In  summary,  the  ways  to  determine  wind  and  density  from  satel¬ 
lites  which  have  been  considered  are  outlined  in  Table  1. 

Of  the  several  techniques,  two  ways  of  determining  density,  by 
measurement  of  atmospheric  emissions  and  of  the  atmospheric  absorption 
of  solar  radiation,  and  one  way  of  determining  wind,  by  the  computa¬ 
tional  simulation  of  the  atmosphere  with  density,  temperature  and 
heating  inputs  derived  by  satellite  measurement,  are  of  major  impor¬ 
tance  to  the  success  of  military  objectives  and  are  unlikely  to  be 
accomplished  by  civil  agencies.  Two  additional  ways  of  determining 
density,  microwave  satellite- to-satellite  occuitation  and  satellite 
drag,  are  important  to  meet  the  needs  of  military  objectives,  but 
not  as  crucially.  They  represent  significant  backup  potential  for 
military  needs.  The  first  may  not  and  the  second  probably  will  not 
be  developed  by  civil  agencies.  One  other  way  of  determining  density, 
laser  backscatter,  is  still  in  an  immature  research  phase  which  merits 
encouragement.  Two  other  ways  of  determining  wind  include  the  deter¬ 
mination  from,  pictures  of  clouds,  characterized  by  limited  altitude 
coverage  and  marginal  accuracy,  and  the  determination  from  the  motion 
of  satellite  cracked  balloons,  which  have  particular  military  vulner¬ 
ability.  Both  of  the  latter  ways  may  perhaps  be  adequately  developed 
by  civil  agencies. 

Also  shown  in  Table  1  are  the  ranges  of  altitude  and  the  uncer¬ 
tainty  of  determination  for  which  the  way  may  be  effective,  and  the 
military  mission  for  which  knowledge  cf  the  atmospheric  parameter  is 
necessary.  For  each  such  military  mission,  the  tolerable  fractional 
uncertainty  of  the  atmospheric  parameter  is  also  shown  to  be  a  number 
equal  or  greater  than  the  uncertainty  of  determination. 
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TABLE  1.  WAYS  OF  DETERMINING  V7INDS  AND  DENSITY  FROM  SATELLITES 


DETERMINATION  OF  DENSITIES 

RANGE  OF 

DETERMINATION 

NEED 

TECHNIQUE 

PARAMETER 

ALTITUDES, 

UNCERTAINTY, 

( APPLICATION/TOLERABLE 

km 

% 

UNCERTAINTY) 

ATMOSPHERIC 

DENSITY 

0-60 

) 

WEATHER  PREDICTION 

10% 

EMISSION 

BALL.  MISS.  TARGET 

10% 

SOLAR  ABSORPTION 

DENSITY 

90-300 

10 

SAT,  ORBIT  PRED. 

10% 

SATELLITE  DRAG  BY 

DENSITY 

10 

SAT.  ORBIT  PRED. 

10% 

ACCELEROMETER 

MICROWAVE 

DENSITY, 

4-20 

0.1* 

WEATHER  PREDICTION 

10% 

OCCULT  ATI  ON 

AIR 

DENSITY, 

WATER 

1 

LASER  BACK5CATTER 

DENSITY,  O 

0-30 

15-100 

WEATHER  PREDICTION 

10% 

DENSITY,  1 

15-100 

DUST 

DENSITY, 

WATER 

5 

DETERMINATION  OF  WINDS 

COMPUTATIONAL 

WIND 

ALL 

25-100 

.  .  » f HER  PREDICTION 

10% 

SIMULATION,  WITH 

ALTITUDES 

BALL.  MISS.  TARGET 

30% 

DENSITY,  TEMPERA¬ 
TURE,  HEATING, 
INPUTS 

BM.D  DISCRIMINATION 

40% 

PICTURES  OF  CLOUDS 

WIND 

1-18 

25 

WCATHER  PREDICTION 

10% 

SATELLITE-TRACKED 

WIND 

1  AND 

3-25 

WEATHER  PREDICTION 

10% 

BALLOON 

10-15 

•BASED  SOLELY  ON  THEORETICAL  HYPOTHESES 


Table  2  lists  the  characteristics  of  satellite  vehicles  and  the 
orbits  required  for  the  application  of  the  several  measurement  tech¬ 
niques.  Comparisons,  rough  and  tentative,  with  existing  satellite 
systems  are  made  in  Table  2.  Confirmation  of  the  comparisons,  and 
further  elaboration  of  system  cetails,  together  with  cost  effective¬ 
ness  comparisons  with  competitive  systems,  are  desirable  objectives 
of  a  further  study. 


Table  3  summarizes  the  results  of  this  review  cf  measurement 
capability,  expressed  as  a  practical  minimum  fractional  uncertainty 
of  the  capability  for  determining  atmospheric  parameters.  Also  shown 
tor  comparison  is  the  tolerable  maximum  in  the  same  altitude  region 
required  for  the  missions  discussed  in  Volume  II . 


i 

i 
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TABLE  2.  TENTATIVE  SATELLITE  SYSTEMS  TO  DETERMINE  ATMOSPHERIC  PARAMETERS 


TABLE  3.  FRACTIONAL  UNCERTAINTIES  OF  MEASUREMENT  TECHNIQUES,  COMPARED  WITH 
NEEDS  FOR  APPLICATIONS  AND  VARIABILITY  OF  CLIMATIC  AVERAGES 


note  FIGURES  WITH  RESPECT  TO  APPLICATIONS  ARE  THE  TOLERARLE  FRACTIONAL  UNCERTAINTY S  EXPRESSED  IN  THE  NUMERATORS 

UNCERTAINTIES  Of  THE  CLIMATIC  AVERAGES  IN  THE  DENOMINATORS.  FIGURES  WITH  RESPECT  TO  MEASUREMENT  TECHNIQUES 
ARE  THf  INSTRUMENTAL  UNCERTAINTIES  IN  THE  NUMERATORS.  UNCE RTAINTlES  Of  THE  CLIMATIC  AVERAGES  IN  THE 
DENOMINATORS.  CHECKS  INDiCAIC  THAT  THE  PARAMETER  IS  SIGNIFICANT  OR  mEASURAUE  RUT  THE  RELEVANT  UNCERTAINTIES 
HAVE  NOT  KEN  EVAIUATEO. 


The  most  significant  conclusion  drawn  from  the  review  of  instru¬ 
mentation  techniques  is  that  satellite-borne  instrumentation  may  re¬ 
motely  measure  all  atmospheric  parameters  considered,  with  accuracy 
and  coverage  adequate  to  meet  the  needs  of  the  important  defense  mis¬ 
sions  considered,  'with  the  exception  of  two  parameters:  winds,  be¬ 
cause  of  inadequacy  in  coverage  and  in  accuracy;  and  electron  density, 
because  of  inadequacy  in  coverage. 

In  addition  to  the  explicit  measurement  of  atmospheric  winds, 
this  report  suggests  that  implicit  determination  of  winds  may  be  made 
by  computational  simulation  for  which  the  input  data  are  the  explicit 
measurements  by  satellite  of  atmospheric  parameters  other  than  winds. 
Justification  for  the  suggestion  is  not,  however,  quantitatively 
shown . 


INVERSION  0?  ATMOSPHERIC 


By  inversion  of  satellite  radiometer  measurements  of  infrared 
emission,  in  the  wing  of  a  collision  broadened  spectral  line,  effec¬ 
tive  height  profiles  of  temperature,  density  and  pressure  may  be  de¬ 
termined  . 


A  characterization  cf  the  technique  is  given  in  Fig.  1,  sum¬ 
marizing  discussions  given  in  the  following  descriptions ,  in  turn,  of 
the  casic  theory,  temperature  and  height  uncertainties,  determination 
of  ozone  density  by  ultraviolet  emission,  determination  of  HgO  and  02 
density  by  microwave  emission,  and  systems  considerations. 

PRINCIPLE 


The  principle  of  the  technique  was  initially  suggested  by  King 
( 1958) ,  elaborated  by  Kaplan  (1959),  and  further  described  by  Wark 
and  Fleming  (1966). 

At  a  single  frequency  v,  the  radiative  transfer  equation  is 


d  I(  v,  CO 


(!•!) 


where  I(v,S)  is  the  radiance  at  v  in  the  direction  9  from  the  local 
vertical,  k(v)  is  the  mass  absorption  coefficient  of  the  absorbing 
gas,  p(z)  is  the  density  of  the  gas,  B(v,T)  is  the  Planck  radiance 
at  v  and  the  temperature  T,  and  z  is  the  height  ir.  the  local  vertical. 


Assuming  that  the  instrument  observes  ir.  a  narrow  cor.e  in  the 
local  vertical  so  that  everywhere  in  the  cone  sec;  -  1,  and  trans¬ 
forming  the  variable  z  to  pressure  p  by  means  of  the  hydrostatic 
equation 


pdz  = 


r\  ra 


(1-2) 


PRINCIPLE  OF  OPERATION 


•  EMISSION  =  f(T,  p,  K,  X)  OF  EMITTING  GAS 

•  MASS  ABSORPTION  COEFFICIENT  (K.)  KNOWN  IN  WINGS  OF  LINE  (  I.E.,  15,1,  COj  ) 

•  PATH  LENGTH  (X)  KNOWN  BY  GEOMETRY 

•  TWO  LINES  OF  C02  AND  02  YIELD  DENSITY  ( p  )  AND  TEMPERATURE  (  T )  AS  FUNCTION  OF  ALTITUDE 

•  ADDITIONAL  LINE  OF  H20  YIELDS  WATER  VAPOR  DENSITY 

UNCERTAINTY  OF  MEASUREMENT 

•  DENSITY:  1% 

•  TEMPERATURE:  1% 

COVERAGE  IN  VOLUME 

•  ALTITUDE:  RANGE  OF  0-30  KM  BY  NADIR  OBSERVATIONS;  30-60  KM  BY  HORIZONTAL  OBSERVATIONS 

•  HORIZONTAL  GRID:  3-200  KM 

STATUS  OF  DEVELOPMENT 

•  I5p-LINE  SENSORS  OF  0-30  KM  VERTICAL  RANGE;  3-200  KM  HORIZONTAL  RESOLUTION  BY  NIMBUS 

•  DATA  PROCESSING:  NIMBUS 

•  DESIGN  NOT  YET  OPTIMIZED 

RECOMMENDATION 

•  INVESTIGATE  FEASIBILITY  OF  HORIZONTAL  VIEW  EMISSION  SATELLITE  FOR  T,  p,  AT  30-60  KM  FOR 

NEEDS  OF  BALLISTIC  WEAPON  TARGETING. 

•  INVESTIGATE  FEASIBILITY  OF  VERTICAL  VIEW  EMISSION  SATELLITE  FOR  MEASURING  3  HR  AND  300  KM 

AVERAGE  VERTICAL  EMISSION  TO  DETERMINE  TOTAL  DENSITY,  TEMPERATURE  AND  WATER  VAPOR 
DENSITY  BY  GEOSYNCHRONOUS  SATELLITE  FOR  NEEDS  OF  NUMERICAL  WEATHER  ANALYSIS  AND 
BALLISTIC  WEAPON  TARGFTING  IN  ALTITUDE  RANGE  0-30  KM.  THIS  MAY  BE  DEVELOPED  ADEQUATELY 
BY  CIVIL  AGENCIES. 


FIGURE  f.  Measurement  of  Atmospheric  Emission  for  Density  at  0-60  km 
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where  q  is  the  mass  nixing  ratio  of  the  absorbing  gas  and  g  is  the 
gravitational  acceleration,  Eq.  1.1  may  be  written  as 


I(v  , o )  =  €  (v  )  B  V  ,T(f„) 


,t(fJ 

'w' 


P)dp 


r' 

*1 


B (v  ,T(p) )  e 


-f/kCv.I 


(1.3) 


)dp 


k(v  ,p)dp 


where  €(v)  is  the  emissivity  of  the  lower  bound  (i.e.,  at  p  )  and  the 
subscript  denotes  the  lower  boundary.  The  first  term  is  the  component 
of  radiance  arising  from  the  surface  ,  and  the  second  term  is  that 
component  arising  from  the  atmosphere  itself. 

It  is  convenient  to  simplify  by  introducing  the  fractional 
transmittance 


t(v  ,p)  =  e 


-n 


k(v,p)dp 


(1.4) 


of  the  beam  between  the  level  p  and  the  effective  top  of  the  atmosphere, 
so  that  Eq.  1.3  becomes 


Kv,0)  =  €(v)  B(v  >T(po))t(v  ,Pq) 


-/ 


t(v,Pq) 


(1.5) 


b(v  ,?(p) )  dT (v ,p) 


Equation  1.5  expresses  radiance  for  a  single  wavelength.  An  instrument, 
however,  distinguishes  only  finite  bandwidths.  The  measured  instrument 
radiance  in  normalized  rorm  for  a  spectrometer  is  t.ne  convolution 
product 


j*  i 


L 


where  w(v)  is  the  spectrometer  slit  function  and  v^,  v2  are  its 
limiting  frequencies. 

A  filter  radiometer,  or.  the  other  hand,  involves  an  inner  product 
rather  than  a  convolution  product. 

Equation  1.6  is  complicated.  However,  if  the  spectral  interval 
v,  to  v2  is  small,  then  B(v,T)  varies  little,  is  linear  in  the  interval, 
and  may  be  replaced  by  its  value  B(v,  T)  at  a  properly  defined  mean 
frequency  v  ,  so  that  it  may  be  factored  out  of  the  integral  -with 
respect  to  v. 

One  may  further  define  the  mean  effective  transmittance  and  its 
partial  derivative  as 
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BP 


v)  dv 

oP 


(1.7b) 


v)dv 


In  the  spectrometers  discussed  the  spectral  intervals  are  narrow 
(5  cm-1  to  20  cm-1)  and  the  slit  functions  w(y)  are  symmetrical  and 
triangular,  so  that  one  may  set  v*  equal  to  the  mean  frequency  v, 
rewriting  Eq.  1-6  as 


I(v,o)  =  B(v,T(pq))  t(v,p0) 


) 

5(v,T(p)) 


dr(v,p) 


(1.3) 


The  purpose  of  the  inversion  is  to  solve  for  the  function  B(v,T(p)), 
knowing  the  radiance  I(v,o)  and  the  transmittance  given  by  Eq.  1*7. 
The  Planck  radiance  B(v,T(p))  has  a  different  form  at  different 
frequencies.  For  example,  ir.  the  microwave  region  the  Rayleigh  Jeans 
approximation  holds 


B(v,T) 


2\/ckT 


Rayleigh  Jeans,  for  microwave  (1.9) 


and  in  the  region  of  the  4.3  micron  carbon  dioxide  band  Wien's  law 
holds  : 

hey 

x  2  ~  krr 

B  (v  ,T )  =  2hv"'c  e  '  *;Wein's  Law,  for  4.3|j,  (1.10) 


In  the  vicinity  of  the  15  micron  carbon  dioxide  band,  Wark  suggests 
the  approximation 

B  (v  , T(p) )  =  a  (v)  B  (vr,T(p) )  +  3  (v  )  for  15u.  line  (1.11) 
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where  subscript  r  indicates  a  fixed  reference  frequency  which  is 
chosen  near  the  middle  of  the  range  of  frequencies.  This  approxima¬ 
tion  holds  only  over  a  limited  part  of  the  spectrum. 

With  the  several  limitations  stated  above,  the  Eq.  1.8  for  the  15 
micron  band  may  be  written  as 


I(v,o)  - 


-  B(ur,T(t))  dt 

J  r\ 


(1.12) 


where  v  is  written  as  v  for  convenience,  and  t  is  a  general  independent 
variable  to  which  pressure  may  be  transformed,  usually  t  =  log  p. 

"  „ion  1.12  is  a  Fredholm  equation  of  the  first  kind 


g(v)  -p 


k(v,t)  f(t)  dt 


wherein 


k(v,t)  =§1 


f(t)  =  B(vr,T(t))  (1.13) 

g(v )  =fl(v,o)  -  b(v)1/q.(v) 


In  the  atmosphere,  the  transmittance  of  Eqs.  1.4  and  1.7  is  a  deter¬ 
mining  factor  in  the  radiance  to  be  observed  from  a  satelli  wc  ci ~c  any 
frequency  and  depends  upon  the  spectral  interval  observed.  Through 
most  of  the  atmosphere  up  to  about  50  km  the  broadening  of  spectral 
lines  is  dominated  by  collisions;  above  that  level  Doppler  broadening 
of  the  lines  dominates.  One  finds  in  the  meteorologically  interesting, 
collision-broadened  atmosphere  that  lines  have  Lorentz  shapes  for 
which 


(1.14) 
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PRESSURE  (mb) 


necessary  that  the  weighting  functions  overlap  a  little,  but  not  too 
much,  in  order  to  obtain  enough  data  to  make  possible  adequate  defini¬ 
tion  of  the  temperature  profile.  With  no  overlap,  the  solution  is 
meaningless  except  that  it  represents  a  mean  weighted  value  of  tempera 
ture  over  the  whole  atmosphere.  With  too  much  overlap,  the  errors  of 
measurement  may  exceed  the  ability  of  the  weighting  functions  to 
discriminate  between  the  levels  of  the  atmosphere  from  which  the  radi¬ 
ation  is  arising. 

Once  the  form  of  the  indicial  function  f(t)  has  been  determined 
(from  the  alternative  Zqs.  1.9,  10  or  11),  one  may  reduce  the  integral 
Eq.  1.13  to  a  system  of  linear  algebraic  equations  by  applying  an 
appropriate  numerical  quadrature  formula.  If,  for  example,  the 
indicial  function  is 

n 

f(t)  e.  cpj(t)  (1.16) 

j=l 

where  cpj(t)  are  a  set  of  orthogonal  functions,  then  with  Eq.  1.16, 

Eq.  1.12  becomes 


n 


3<V  =  £ 


5=1 


,t)  cpj  (t)  dt 


t  =  1, 


n  (1.17) 


or  in  matrix  notation 


Inverting  yields 

c  =  A"1  g  (1.18) 

where  A  is  the  inverse  of  the  matrix  A  whose  elements  are  given  by 
the  integrals  of  1.17,  g  is  the  vector  exp 


r ess  ion  of  the  observed 


quantities  g(  ,  and  c  is  the  vector  of  undetermined  coefficients, 
which  then  evaluated  are  used  to  determine  B(T)  =  f(t),  and  from 
thence  temperature  T(t). 

The  solution  of  Eq .  1.18  is  unstable  if  more  than  three  equations 
are  considered,  due  to  the  following  sources  of  error:  the  approxi¬ 
mation  of  the  Planck  function,  the  numerical  quadrature  process  and 
roundoff. 

The  computational  instability  has  been  made  tractable  by  use 
of  the  solution 

c  =  (AT  A  +  y  H)'1  (AT  g  +  y  h)  (1.19) 

<•**/  's/ 
where  A  is  the  transpose  of  A,  y  is  the  smoothing  parameter,  H  is 

the  smoothing  matrix  and  h  is  a  bias  vector  which  is  a  gross  estimate 

of  the  solution  vector  c.  Fleming  and  Wark  (1966)  find  that  y  =  1°  5 

is  pragmatically  effective. 

TEMPERATURE  UNCERTAINTY 

Wark  and  Fleming  (1966)  have  made  estimates  of  error  of  temper¬ 
ature  determinations  of  a  six  channel  (of  5  cm  1  bandwidth)  spectrometer 

-1  -2  -1  -1 

measuring  radiances  of  about  50  ergs  sec  cm  ster  cm  .  For 

errors  of  1,  2, and  4  percent  in  the  measure  of  radiance,  corre- 

-1  -2  -1  -1 

sponding  to  a  =  0.25,  0.5,  1.0  and  2.0  ergs  sec  cm  ster  cm  , 
of  a  set  of  110  random  values ,  the  errors  of  the  computed  Planck 
function  B,  and  of  the  computed  temperature  T  associated  with  B  by 
the  relation 

(1.20) 

are  given  in  Table  4.  In  the  table,  the  columns  marked  "maximum"  are 
for  the  fifth  largest  value,  the  units  for  a,  B  ,  &B  and  Max  AB  are 
ergs  sec  *  cm  z  ster  cm  1  and  the  units  for  T,  AT,  and  Max  AT  are 


TABLE  4.  TEMPERATURE  ERROR  RESULTING  FROM  RANDOM  ERROR  OF  OBSERVATION 

(AFTER  WARK  AND  FLEMING,  1966) 


c 

V 

5  ±  3 

T  ±  Li 

Max . 

£B 

Max . 
L_jt 

KEY  WEST 

0 

io"4 

47.59  +  0.96 

224 .65  +  1.01 

0.96 

1.01 

0.25 

0.025 

47.59  ±  1,80 

224.65  ±  1.90 

2.50 

2.63 

0.5 

0.050 

47.60  ±  2.01 

224.66  ±2.12 

2.85 

3.00 

1.0 

0.150 

47.61  ±  3.04 

224.67  ±3.20 

4.63 

4. 88 

2.0 

0.400 

47.63  ±  4.42 

224.69  ±4.66 

6.36 

6.70 

LITTLE 

ROCK 

0 

10  "4 

47.39  ±  0.93 

224.43  +  0.98 

0.93 

0.98 

0.25 

0.025 

47.39  ±  1.21 

224.43  ±  1.28 

1.69 

1.79 

0.5 

0.050 

47.40  ±  1.54 

224.44  ±  1.63 

2/55 

2.69 

1.0 

0.150 

47.41  +  2.13 

224.46  ±  2.25 

3.62 

3.82 

2.0 

0.400 

47.43  i  2.81 

224.48  x  2.97 

4.38 

4.63 

SEATTLE 

0 

10 

43.84  t.  0. 69 

220.60  ±  0.76 

0.69 

0.76 

0.25 

0.025 

43.84  ±  1.60 

220.60  ±1.77 

2.47 

2.73 

0.5 

0.050 

43.85  i  1.89 

220.61  ±  2. OS 

3.27 

3.61 

1.0 

0.150 

43.85  +  2.45 

220.62  ±  2.70 

4.47 

4.93 

2.0 

0.400 

43.88  ±  3.05 

220.64  ±3.37 

5.25 

5.79 

CHURCHILL 

0 

10 

35.13  ±  0.86 

210.36  ±1.08 

0.86 

1.08 

0.25 

0.025 

35.13  ±  1.71 

210.36  ±2.14 

2.24 

2.81 

0.5 

0  .050 

35.14  ±1.96 

210.38  ±  2.46 

2.89 

3.63 

1.0 

0.150 

35.15  ±  2.79 

210.39  ±  3.50 

4.65 

5,84 

2.0 

0.400 

35.17  ±  4.14 

210.41  ±5.20 

6.25 

7.84 

From  the  table  it  may  be  seen  that  errors  or  radiance  equal  to 

-1  -2  -1  -1 

or  greater  than  about  1%  or  0.5  erg  sec  cm  ster  cm  lead  to 
results  with  temperature  error  greater  than  l'T  or  about  2  °K  . 
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W.  L.  Smith  (196.8)  has  made  an  estimate  of  the  errors  of  the 
atmosphere's  temperature  versus  pressure  distribution  as  determined 
from  satellite  radiation  measurements,  using  as  an  instrument  the 
Satellite  Infrared  Radiation  Spectrometer  (SIRS)  which  has  flown  in 
the  NIMBUS  III  meteorological  satellite. 

SIRS  detects  emissions  in  5  cm-'1'  bands  at  the  five  CCL  channels 
near  15  microns  (669,  677.5,  692,  699,  and  706  cm  )  as  well  as  two 
channels  in  the  far  wing  (at  714  and  748  cm- ^ )  and  an  atmospheric 
window  at  899  cm  \  The  uncertainty  of  temperature  determination 

by  SIRS  is  shown  in  Fig.  4  for  the  case  of  both  clear  and  cloudy 

-1  -2  -1 

conditions,  and  for  errors  of  0.01  and  0.25  erg  sec  cm  ster  cm 

-1  -2 

in  the  determination  of  radiance  values  from  50  to  70  erg  sec  cm 
ster-^"  cm-\  The  temperature  uncertainty  is  seen  to  be  less  than 
3°K. 


FIGURE  4.  Uncertainty  of  Temperature  Determination  by  SIRS 
(after  W.L.  Smith,  1968) 
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HEIGHT  UNCERTAINTY 

W.  L.  Smith  (1968)  also  describes  a  method  of  estimating  the 
pressure  versus  height  relation  from  the  radiometric  measurements, 
employing  empirical  orthogonal  functions  developed  by  Holmstrom  (1963). 
By  the  Holmstrom  technique,  the  geopotential  height  is  represented 
by  the  rapidly  converging  series 


Z(p) 


•fc 


cr  *R(p) 


(1-21) 


Holmstrom' s  studies  of  eight  days  data  indicate  that  for  meteorological 
studies  in  the  pressure  range  100  to  1000  mb,  the  empirical  ortho¬ 
gonal  functions,  cpR(p),  are  statistically  stable  in  being  nearly 
invariant  with  respect  to  geographical  position  and  time,  over  the 
space  and  time  intervals  usually  considered  by  numerical  weather 
prediction.  Prom  the  hydrostatic  Eq.  1.2,  and  the  gas  lav?,  used  with 
Eq.  1.21,  it  follows 


T(p) 


az(p)  2  V  c  v^R(r) 
htnv  '  RL  R  S-tnp 
k=l 


k 

=  U  Cr  Xr(?) 

In  matrix  notation,  Eqs.  1*21  and  1.22  may  be  written 
Z  -  I  c 


wherein 


*  -  V5 

*  =  (xjk) 


(1.22) 
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Combined,  these  yield 


— ^  T*  “  I  '■'w  T1  — *  '■>«*  — ♦ 

^  =  X(S  $)  S*Z  =  DZ  (1.23) 

which  states  the  geopotential  height  profile  is  uniquely  related  to 
the  temperature  profile.  Then,  to  the  extent  that  Eq.  22  is  valid, 
the  height  distribution  can  be  specified  solely  from  the  temperature 
profile.  W.  L.  Smith  (1968)  estimate  of  uncertainty  of  height 
determination  by  SIInS,  for  which  the  temperature  uncertainty  is  given 
in  Fig.  4,  is  shown  for  similar  case  in  Fig.  5. 

Assuming  that  the  empirical  orthogonal  functions  evaluated 
statistically  at  the  time  are  invariant,  the  uncertainty  of  height 
due  to  the  inversion  computation  alone  is  seen  to  be  less  than  60 
meters,  or  expressed  fractionally  less  than  0.002  to  0.01. 


GEOPOTENTIAl  HEIGHT  DEVIATION  (METERS) 

FIGURE  5.  Uncertainty  of  Height  Determination  by  SIRS 
(After  W.L.  Smith,  1968) 


By  the  equation  relating  one  uncertainties  of  density  anc  height 


c 


(1.24) 


wherein  H  is  the  scale  height,  about  7  km,  the  fractional  uncertainty 
of  density  is  about  1)4. 

Since  the  linear  spatial  resolution  on  the  ground  of  the  SIRS  is 
about  225  km,  the  interval  between  clouds  cannot  be  resolved  by  a 
single  observation.  However,  according  to  Suomi  (1969),  the  NASA  ATS 
satellites  have  demonstrated  that  time  and  space  can  be  traded;  that 
is ,  a  90%  cloud  cover  could  be  effectively  sounded  to  the  ground  in 
one  hour  of  observation.  By  means  of  such  a  tradeoff,  a  1C  to  15  km 
resolution  radiometer  makes  possible  altitude  coverage  all  the  way  to 
the  ground.  A.  W.  McCulloch  and  W.  L.  Smith  (1552)  have  proposed  a 
high  resolution  temperature  sounder,  with  7.5  km  ground  resolution, 
by  which  the  intervals  between  clouds  may  be  resolved,  thereby  per¬ 
mitting  resolution  c?  smaller  temper. mure-  and  height  uncertainties 
described  as  "clear,  cy  =  0.2?”  in  Figs.  4  and  5,  namely  less  than 
2  °K  end  -C  meters  ever  most  of  the  altitude  range,  0  to  16  km.. 


As  nay  be  judged  by  inspection  of  Figs.  3,  4,  and  5,  the  deter-' 
minaticn  of  temperature  and  pressure  or  height  by  the  inversion  of 
nadir  directed  radiometer  measurements  of  the  15  micron  CC^  line  is 
effective  only  between  the  earth's  surface,  or  the  top  of  the  cloud 
layer  if  uninterrupted,  anc  the  height  of  the  10  mb  surface,  approxi¬ 
mately  30  km..  J.  C.  Gille  (195S)  has  suggested  the  inversion  cf 


horizon -directed  radiometer  measurements  of  the  1 
determine  the  temperature,  ar.c  the  pressure  cr  he 
pause,  from.  30  to  60  km.  Since  the  absorption  is 
density-path  length  integral,  the  horizon -directs 
sensitive  to  the  atmosphere  a 
directed  measurements,  anc  the  weighting 
shape ,  but  much  narrower,  than  those  of 
of  uncertainty  is  base;  on  the  icilov.dr.c 


igher  altitudes  t 
ric?  : ur.cticri 


5  micron  00^  line  to 
ight,  in  the  strato- 
c report ional  to  the 
d  measurements  are 
han  are  the  nacir- 
~  are  similar  in 
J .  C.  Gille' a  estima 
ations .  Assumina  a 


satellite  at  ■-0C  V.ir.  altitude .  the  vertical  -dimension  for  a  1  milli- 
ra.liar  field  cf  viov:  is  3j prcxir.iarely  d.--  kit.  Over  this  height  the 
curvature  o:  tat  radiance  profiles  is  snail,  so  that  changes  nay  be 
easily  resolved  it.  the  field  cf  vie.:,  tille  considers  a  cooled  or¬ 
biting  optical  system  and  detector  capable  c-f  measuring  a  full  seal 

“2—1  -2  -1 

signal  of  7  w  nr  ster  x  to  an  uncertainty  of  0.01  w  m  ster  ,  or 

C.14T.  Point inc  errors  nav  be  the  liir.it inc:  factor  in  determining 


resolution  because  of  the  large  radiance  profile  slope,  about  0.2 

-2-1-1 

v:  m  ster  J'  km  .  If  the  pointing  error  is  one  mi  Hi  radian,  the  re¬ 
sulting  uncertainty  of  tangent  height  v.iil  approximate  0,4  km  yield! 

_  2  _  i 

error  of  about  C.08  w  m  ster  ,  or  about  15.  The  two  error  source 
(in  radiance  measurement  and  pointing)  may,  therefore,  contribute  to 
a  total  measurement  error  of  about  lb.  The  contribution  of  error  du 
to  the  computational  inversion  may  be  expected  to  be  the  same  as  in 
the  case  of  the  nadir-directed  radiorr.etry . 

Mcr-.ee  et  al.  (xSo9)  have  analyzeo  the-  errors  or  temperature  de¬ 
termination  from  horizon  radiance  profiles  of  the  16  micron  CO,  line 

•  c 

for  the  altitude  range  23.3  to  55.5  km.  McKee  et  -1.  find  that  the 
error  sources  shown  in  Table  5  are  realizable,  and  contribute  errors 
in  temperature  determination  as  indicated. 


TABLE  5.  EFivTRS  CF  TEMPERATURE  FROM  HORIZON  RADIAIOE  PROFILES 

(McKee  et  al.,  1SS8; 


ERROR 

AT  ALTITUD 

20  KM 

AT  ALTITUDE  50  KM 

FOR  SIGNAL 

FOR  SIGNAL 

TYPE 

MAGNITUDE 

w  m"2  ster"’ 

AT 

°k 

w  m'2  ster-1 

mm 

SCALE  ERROR 

(0.01  SIGNAL) 

5.92 

0.5 

1.10 

1.16 

BIAS  ERROR 

0.01 

0.03 

0.12 

W  m'2  ster"1 

NOISE  ERROR 

0.01 

2.30 

0.93 

w  m'2  ster"1 

TANG.  HEIGHT 
ERROR 

0.5 

1.35 

3.80 

wm"2  ster"’ 

PRESSURE 

ERROR 

(0. 1  PRESSURE  ) 

1.52 

5.79 

DENSITY  OF  OZONE  BY  ULTRAVIOLET  EMISSION 

Other  radiation  bands  rr.ay  be  used  similarly.  At  the  low  alti¬ 
tudes,  be low  30  km,  two  interesting  constituents  of  the  atmosphere 
which  are  highly  variable  in  mixing  ratio  are  water  vapor  and  ozone . 
As  is  shown  in  Fig.  6,  the  absorption  coefficient  of  water  vapor  in 
the  infrared  is  about  two  orders  of  magnitude  greater  than  that  of 
carbon  dioxide,  which  is  nearly  constant  in  mixing  ratio.  Radiometer 
me  a  s  ur  emeries  in  either  of  the  two  bands  shown  in  Fig.  6  may 
yield  data  on  the  mixing  ratio-height  distribution  of  water  vapor. 
Hov.’ever,  since  the  ozone  band  absorption  coefficient  is  of  the  same 
approximate  size  as  the  CO^  band,  and  the  mixing  ratio  is  far  smaller 
in  the  altitude  range  of  interest,  still  another  band  of  radiation 
should  be  considered  --  that  of  the  Hartley  band  in  the  ultraviolet  -- 
for  the  determination  of  the  ozone  mixing  ratio-height  distribution. 


wavelength  X(  microns  ) 

100  SO  30  20  16  12  10  8  6  5  4. 


0  200  600  1000  1400  1800  2200 

WAVE  NUMBER  V  (CM'1  ) 

FIGURE  6.  Generalized  Absorption  Coefficients  for  the  Major  Absorbing  Bands 


HO,  CO„  and  0_  (After  Elsusser  and  Culbertson,  1960) 

A  ^  O 


S.  Twomey  (1961)  has  described  the  problem  of  manipulating 
ultraviolet  spectral  measurements  in  the  -wavelength  region  2300  to 
3000  A  for  the  determination  of  ozone  mixing  ratio-height  distribution 
from  the  solar  radiation  backscattered  from  above  the  median  ozone 
level  (abouc  25  km).  Twomey  anticipated  that  the  error  of  ozone 
determination  can  be  made  less  than  1%.  Herman  and  Yarger  (1969) 
have  elaborated  the  concept  to  indicate  the  method  should  be  effective 
over  the  range  2  to  40  km. 

DENSITY  OF  H20  and  02  BY  MICROWAVE  EMISSION 

A  third  region  for  application  of  inversion  of  radiometer 

measurements  if  the  atmospheric  microwave  spectrum,  notably  the  weak 

g 

resonance  of  water  vapor  near  22  GHz  (1  GHz  =  10  cps)  or  1.35  cm 
wavelength  and  the  strong  complex  of  oxygen  resonances  near  60  GHz 
(0.5  cm).  Atmospheric  emission  at  these  frequencies  can  be  received 
with  microwave  radiometers ,  and  the  received  radiance  expressed  in 
terms  of  the  brightness  temperature  ?B  may  be  measured.  The  accuracy 
with  which  T^  may  be  measured  is  0.2  to  1.0  °K.  Since  the  Planck 
function  B  is  essentially  linear  in  this  spectral  region,  as  indicated 
in  Eq.  1.9,  the  inversion  of  Eq.  1.8  is  straightforward. 

In  the  nearly  opaque  regions  of  the  spectrum  the  temperature 
profile  may  be  obtained  in  a  manner  exactly  analogous  to  that  in  the 
infrared  spectral  region.  In  the  nearly  transparent  regions  water 
vapor  and  perhaps  ozone  may  be  inferred  in  some  circumstances , 
especially  oversea  areas.  Because  a  microwave  radiometer  looking 
down  may  see  water  vapor  or  ozone  in  absorption  or  emission, 
dependent  on  the  brightness  temperature  of  the  surface,  the  highly 
emissive  land  areas  are  characterized  by  a  brightness  temperature 
nearly  the  same  as  the  atmospheric  temperature.  Over  water  areas, 
the  surface  emissivity  (6(v))  of  Eq.  1.3  and  the  surface  brightness 
temperature  (effect  of  the  first  term  of  Eq.  1.3)  are  low,  and  the 
atmospheric  emissions  (the  second  term  of  Eq.  1.3)  are  unambiguous. 
Wtters  of  MIT,  cited  by  C-ille  (1968),  has  quoted  the  inversion 


uncertainties  associated  with  zero  and  1%  errors  of  microwave  radio¬ 
meter  measurements  to  be  2  to  3.7  °K  (about  1%)  for  temperatures 

_2 

at  altitudes  below  15  km,  and  to  be  0.09  to  0.11  gm  cm  (about  3 
to  30%)  for  integrated  water  vapor,  and  that  the  height  resolution 
is  comparable  to  that  of  the  15  micron  infrared  channels;  that  is 
the  emission  by  oxygen  in.  a  spectral  channel  of  0.2  GHz  width  comes 
mainly  from  an  atmospheric  layer  about  10  km  thick. 

SYSTEMS  CONSIDERATIONS 

Although  systems  studies,  involving  cost-effectiveness  compari¬ 
sons  with  physically  feasible  alternatives,  are  beyond  the  scope  of 
the  study  effort  reported  in  this  volume,  nevertheless  some  analogies 
to  proven  systems  can  be  made.  The  techniques  for  measurement  of 
atmospheric  emission,  as  probes  of  the  altitude  region  up  to  30  km, 
have  been  largely  proven  by  NIMBUS  and  TIROS  satellites,  the  history 
and  characteristics  of  which  are  described  in  some  detail  in  Appendix 
A  of  this  volume.  The  700  lb  NIMBUS  and  300  lb  TIROS  satellites  orbit 
sun- synchronously  at  altitudes  of  4CC  to  600  nmi ,  making  measurements 
continuously  along  a  swathe  about  50C  nmi  wide  along  the  orbital 
projection  on  the  earth's  surface.  Eight  such  satellites,  with 
ascending  nodes  at  45°  intervals,  are  required  to  obtain  measure¬ 
ments  over  all  the  earth  at  three  hour  intervals. 

By  means  of  a  system  cf  similar  sensor  complexity  with  the 
addition  of  a  directionally  stabilized  pointing  control,  in  similar 
sun- synchronous  orbits,  the  atmosphere  can  be  probed  for  density  and 
temperature  at  all  altitudes  up  to  about  50  km.. 

SUMMARY 

In  summary,  inversion  of  radiometric  measurements  from  an  earth 
orbiting  satellite  yields  the  distribution  with  respect  to  altitude, 
of  temperature,  total  particle  density,  pressure,  and  density  of 
water  vapor,  carbon  dioxide,  molecular  oxygen  and  ozone,  with  frac¬ 
tional  uncertainty  between  one  and  two  percent,  as  indicated  in  more 
detail  in  Table  6. 


TABLE  6.  INVERSION  OF  ATMOSPHERIC  EMISSION 
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One  technique  by  which  the  vertical  distribution  of  the  temper¬ 
ature  aiid  density  of  the  atmosphere  may  he  determined  remotely  by 
satellite-borne  instrumentation  is  the  observation  of  the  variable 
absorption  of  solar  radiation  as  the  satellite  crosses  the  termi¬ 
nator,  that  is,  as  the  line  of  sight  from  the  satellite  to  the  sun 
is  occulted  by  the  solid  earth  and  its  atmosphere.  By  suitable  choice 
of  lines  or  wave  bands  of  the  solar  spectrum  which  are  selectively 
absorbed  by  the  earth  atmosphere,  the  density  and  temperature  of 
neutral  components  of  the  atmosphere  may  be  determined  in  the  upper 
atmosphere  (e.g.  ,  from,  about  30  to  300  km.  altitude)  where  ether 
techniques  arc  often  inapplicable. 

A  characterization  of  the  technique  is  given  in  Fig.  ?,  which 
summarizes  considerations  reviewed  in  the  following  paragraphs,  deal¬ 
ing  :n  turn  with  basic  theory,  characteristics  cf  the  solar  radiation 
wave  bands,  detector  characteristics  and  vehicle  and  system  consider¬ 
ations  . 

RADIATIVE  TRANSFER 

The  principles  of  radiative  transfer  which  govern  the  technique  of 
determining  density  by  observing  the  absorption  of  solar  radiation  at 
various  depths  in  the  atmosphere  differ  from  those  described  in  the  dis¬ 
cussion  of  the  inversion  of  atmospheric  emission  in  that  the  emission  of 
the  atmosphere  is  taken  to  be  negligible  by  comparison  with  that  of  the 
solar  source.  Ir.  consequence,  the  spectral  radiance  from  the  sun 
detected  at  the  satellite  may  be  described  as 

Cs(X)  =  Eo(X)  expj ~J  cn(z)F^-p,  njcz 


(z.i) 


PRINCIPLE  OF  OPERATION 


•  ABSORPTION  =  f(p,  K,  X)  OF  ABSORBING  GAS 

•  MASS  ABSORPTION  COEFFICIENT  ( K.)  KNOWN  IN  X-RAY  AND  EUV 

•  PATH  LENGTH  (X)  KNOWN  BY  GEOMETRY 

•  DETERMINES  DENSITY  (p) 

UNCERTAINTY  OF  MEASUREMENT 

•  DENSITY:  10% 

•  TEMPERATURE:  10% 

COVERAGE  IN  VOLUME 

•  DENSITY  AVERAGED  OVER  VOLUME  7  KM  x  300  KM  x  1 .6  KM 

•  TWO  VERTICAL  PROFILES  FROM  90  TO  300  KM  ALTITUDE,  PER  ORBIT 

•  NINE  SATELLITES  FOR  DAILY  SAMPLE  IN  1800  KM  GRID  SPACING 

STATUS  OF  DEVELOPMENT 

•  PRINCIPLE  DEMONSTRATED  BY  ROCKET  EXPERIMENTS 

•  DATA  ANALYSIS  BY  COMPUTER  DEMONSTRATED 

RECOMMENDATION 

•  INVESTIGATE  FEASIBILITY  OF  SATELLITE  SYSTEM  FOR  NEEDS  OF  LOW  ALTITUDE  SATELLITE  ORBIT 

PREDICTION,  USEFUL  FOR  RADIO  COMMUNICATION. 


ATMOSPHERIC  ABSORPTION  OF  SOLAR  RADIATION 

FIGURE  7.  Measurement  of  Absorption  of  Solar  Radiation  for  Density  at  90-300  km 


'wherein  E  is  the  spectral  radiance  at  a  given  wavelength  X,  at 
altitude  h 

Eq  is  the  spectral  radiance  of  the  sun  outside  the  absorbing 
atmosphere 

o 

■z  is  the  absorption  cross  sectror.  (cm  ) 
n  is  the  number  density  of  the  absorbing  gas  (cm"-*),  at 
altitude  z 

X  is  the  wavelength  cf  monochromatic  radiation 
F  is  the  Chapman  function  (of  a,  z,  H  and  3) 
a  is  the  radius  cf  the  earth 
3  is  the  solar  zenith  angle 
H  is  the  atmospheric  scale  height 
z  is  the  altitude  in  the  atmosphere 

The  Chapman  function  F^a  2 ,  3  ^  takes  into  account  that  the  incidence 
of  radiation  is  not  normal.  In  rhe  relation  of  Ec.  25,  if  the  latter 
were  written  for  a  plane  stratified  atmosphere,  F ,  3^  becomes 
sec  3.  It  has  values,  in  the  case  of  the  spherically  stra'tified 
atmosphere  of  the  earth,  that  are  tabulated  in  Swidsr  (1954)  (e.g., 

F  =-  39  ar.d  sec  3  =  «  at  3  =77).  The  power  detected  by  a  radiometer 
is  represented  by  the  integral  over  the  inner  product  of  the  incident 
spectral  radiance  and  the  spectral  efficiency  of  the  detector,  over 
a  band  wavelengths  to  which  the  radiometer  is  sensitive. 

If  the  radiant  flux  is  characteristically  that  of  a  black  body, 
then  the  incident  radiance  over  a  wave  band  lira  ted  by  ,\9  detected 
by  an  ior.  chamber  sufficiently  excited  to  pass  current  I  within  limits 


zv 


wherein 


T  =  ion  chamber  current 

i  =  photon  detection  efficiency,  within  limits 
Xmin’  'V.5x 

E(\15\2)  =  energy  flux  in  the  band  of  wavelengths  X^>  \2 

x  =  ion  pairs  formed  in  the  ion  chamber  per  erg  of 
radiant  energy 
X  =  wavelength 

a  =  aperture  area  of  detector 
q  =  electron  charge 
c  -  velocity  of  light 
h  =  Planck’s  constant 
C,  =  1.439  cm  °K 

T  =  temperature  of  black  body  source 


For  the  case  when  in  Eq.  2.2,  the  bands  of  energy  flux  to  be  determined 
and  those  of  the  detector  sensitivity  are  the  same,  then  the  flux 
is  proportional  tc  the  current  in  the  detector 

(2.3) 

yqyjai  J 

wherein  c  is  an  av* -raged  value  of  photon  detection  efficiency. 

Considering  now  the  flux  at  an  altitude  h,  over  a  wave  band 
(limited  by  \1 ,  \9),  as  the  integral. of  Fq .  2.1 


din  E(h) 


1  d  E(h) 
E  cz 


( X  )  ex  pf"'|  <j  (x  >z)n(z 


E0(x)  exp (-[ <.(X,z)n(z)F(a_p,s)dzja>. 

'•'h  (2.5) 


(2.5) 


If  ever  the  narrow  band  cf  sensitivity  of  the  radiometer  Eq(\ )  is 
taken  as  that  of  a  black  body,  then  Sq.  2.5  may  be  written,  witr.  z 

in  lieu  cf  k,  as 


x2  exo  -/  a (X ,z)n(z)az 

f  .  .  i 

J  t(X,z)  5  /  c2/Vf  T 


din  E(z) 


=  n(z)  F(z) 

az  \ 


1  \ 


-  (X  ,z)r.(z)az^c 


which  the  density  may  be  evaluated  by 


.  .  din  E(z)/dz 

n(z)  =  ~r - 

c(z)  F(z) 


(2.5) 


(2-7) 


* git.  o  ( z')  i s  3.  mean  v  si.u 
brackets  ir.  Eq.  1.6. 


lue  evaluated  as  the  quantity  in  square 


?.lso  the  density  over  a  snort 
evaluated  in  terms  of  scale  height 


’nort  interval  of  altitude  may  be 


_  k  T(z) 


as 


n  -  n  e 
o 


-z/H 


(2.S) 


so  -chao 


H(z)  =  - 


/dtn  n(z)\ 

\  dz  ! 


(2.10) 


where in  T  is  temperature 

z  is  altitude 
k  is  Boltzmann's  constant 
g  is  the  acceleration  due  to  gravity 
rr.  is  the  mean  molecular  we ight  of  the  atmosphere. 

With  the  description  of  density  given  as  in  Eq.  2.2,  the  temperature 
may  be  determined  by  the  derivative  relation 


XU)  .  - 1  (^IJ1  (su,)-1 


(2.11) 


Ir.  short,  from  a  height  profile  of  radiance,  determined  from 
radiometer  measurements  by  Eq.  2.2,  the  height  profile  cf  density 
may  be  determined  by  Eq.  2.7,  the  height  profile  of  scale  height  by 
Eq.  2.10,  and  if  the  mean  molecular  weight  m(c)  be  known,  the  temp- 
ature  profile  may  be  determined  by  means  of  Eq.  2.11. 

The  certainty  cf  an  individual  data  point  associated  with 
determination  of  each  of  the  parameters  for  which  profiles  in 
height  are  made  may  be  described  by  the  following  relations. 

Assuming  the  measurements  of  detector  current  I  and  detector 
efficiency  are  independent ,  the  variance  in  determination  of  flux 
is  described  bv 


« ■  W  •  (f 


(2.12) 


Similarly  assuming  independence  of  the  determinations  of 
radiance,  absorption  cross  section  and  Chapman  constants,  and  that 


32 


the  independence  of  radiance  and  radiance  gradient  is  described  by 
correlation  pF ,  the  variance  of  the  determination  of  density  is 
giver,  by 


I .  J.  Ui  L  -k. » Iki  j  1 1 ;  !l ;  1. 1 .  J  Ui  UL_iM ,  lu  U  ki  jJlijUilll  1 1 


ALTITUDE,  IN  KM 


As  shown  by  S 'wider  (19&4),  the  Chapman  function  from  a  satellite 
at  altitude  200  km,  observing  the  solar  radiation  passing  through  a 
tangent  height  of  60  to  200  km,  takes  values  from  about  30  to  about 
90. 

The  absorption  cross  sections,  for  the  principal  constituents 
of  the  atmosphere  in  the  altitude  region  60  to  200  km  (e.g.,  0,  0  , 
and  N, )  are  summarized  for  wavelengths  1  to  2000  A  in  Fig.  11  from 
ok  ha  published  by  Hinteregger,  et  al.  (1964  ). 

The  following  are  examples  of  diagnostic  wavebands,  for  which 
computations  of  effective  optical  depth  are  given  in  Table  7.  Shown 
in  Fig.  12  is  the  absorption  of  solar  flux  in  the  bands  10  -  20  A, 
near  the  Hydrogen  Lyman  Alpha  (1215.7  A)  line,  and  in  the  Schumann- 
Runge  continuum  (1375  -  1500  A).  The  first  of  these  absorbs  each  of 
the  three  principal  constituents  0,  0^,  and  N?.  Under  the  assumption 
that  near  the  turbopause  diffusion  has  not  appreciably  affected  the 
proportions  of  oxygen  and  nitrogen,  the  one  waveband  alone  is  diagnos¬ 
tic  oi  the  neutral  particle  density  in  the  region  105-135  km.  The 
latter  two  are  diagnostic  only  of  the  density  of  molecular  oxygen  (O2), 
to  cover  the  region  00  -  115  km  nd  125  -  145  km,  respectively.  The 
data  of  Fig.  12  were  determined  at  Eglin,  Florida  (Lat.  30  A°N  at 
07(jy,  r'9  CST  on  21  February  1966.  Zenith  angle  was  83°,  and  Chapman 
function  F  V. 6.  Measured  from  a  satellite,  the  altitudes  for  which 
the  three  examples  are  diagnostic  may  be  perhaps  two  scale  heights 
higher,  or  for  total  particles  (N2,  O2  and  O),  the  10  -  20  A  band  is 
diagnostic  in  the  altitude  range  130  -  160  km,  anc  for  0 2  the  ultra¬ 
violet  bands  arc  diagnostic  in  the  altitude  range  from  about  98  km 
to  the  turbopause  r-.-ar  112  km  and  from  about  152  to  198  km,  respec¬ 
tively. 

A  'waveband  that  is  diagnostic  of  0„  alone  at  altitude  'arhaps 
two  sea]*,  heights  above  that  of  the  Schumann-Rvnge  continuu  . j  that 
encompassing  the  line  oi  Hydrogen  Lyman  Beta  (1025.7  A).  Absorption 
oi  II  Ly  p,  in  vertical  incidence  has  been  described  by  L.  A.  Hall 
wt  a]  ,  ,  Vj'.'U,  as  iii  i  ig .  1  h 


3C 


TOUl  C*OSS  SECTION  (CM1) 


I LOMETER5 


The  absorption  cross  section  by  CL  for  H  Ly  g  ,  3^  =-  1.6(10) 

2  ^  °2 
cm  ;  by  0  and  by  N2 ,  the  absorption  cross  section  is  negligibly 

small.  From  a  satellite,  the  H  Ly  p  is  diagnostic  of  the  density 

of  in  a  region  about  three  scale  heights  above  that  shown  in 

Fig.  13,  or  from  about  130  to  160  km.  Density  solutions  obtained 

simultaneously  from  measurements  of  the  absorption  of  10  -  20  A  and 

H  Ly  Beta  yield,  if  diffusive  separation  is  assumed  negligible,  the 

hoicrht  distribution  of  H,. ,  0,  and  0  in  the 
-  / 


altitude  range  130  to 


PHOTON  FLUXES,  <MX  ,t  X,h)  IN  PH  CM'2$EC~' 


I _ I _ I _ I _ I  I  l  '  i  1 _ I _ | _ l_L 

100  500  1000  5000 
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FIGURE  13.  Flux,  F;.tinction  Coefficient,  and  Absorption  of  H  Ly  B  1025.7  A 
(After  L.  A.  Hall  et  al.,  1962) 
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it  the  same  altitude  range  for  'which  the  Schumann -Runge  continuum 

is  niacr.ostic  for  C-,  the  abscrution  of  23  -  33  A  radiation  is  ab- 

/ 

sorbed  by  both  and  C„.  Here  aaair. ,  if  it  can  be  assumed  that  cif - 
fusive  separation  is  negligible,  density  solutions  obtained  simul¬ 
taneously  from  measurements  of  the  two  tar.cs  yield  the  height  distri¬ 
bution  of  hg!  G0  and  0.  The  assumption  of  negligible  diffusive 
separation  is  a  weak  one,  as  may  be  seen  by  inspection  of  Fig.  8. 

In  the  altitude  range  21?  to  248  km,  simultaneous  density  solu¬ 
tions  from  measurements  of  the  absorption  of  the  lines  at  787.7A(0  IV), 
790A(0  IV)  and  236.3  A  (He  II)  may  be  used  to  determine  unambiguously 
the  height  distributions  of  i^,  Co  and  0.  Also  in  the  altitude  range 
237  to  263  km,  simultaneous  density  solutions  from  measurements  of 
the  absorption  of  the  lines  at  256.3  A  (He  II),  368.1  A  (Mg  IX),  and 
584.3  A  (He  I)  may  be  used  tc  determine  the  height  distributions  of 
N25  O2  ana  0.  Probably  representative  of  such  measurements  are  the 
experimental  data  reported  by  L.A.  Hall  et  al.  (1962)  as  described 
in  Fig.  14  for  vertically  incident  790  A  radiation  and  in  Fig,  15 
for  vertically  incident  368.1  A  radiation. 

Summarized  in  Fig.  16  are  the  altitude  regions  for  which  the 
density-height  distribution  of  principal  atmospheric  constituents 
may  be  determined.  In  these  same  regions,  by  use  of  the  temperature- 
density  relation  of  Eq.  2.11,  the  temperature-height  distribution  may 
also  be  determined. 

DETECTOR  CHARACTERISTICS 

The  detectors  described  by  L.A.  Hall  et  al.  (1962)  as  employed 
by  H.  E.  Hinteregger  and  his  group  at  AFCRL  for  measurements  in  the 
region  of  wavelength  250  to  1300  A  are  monochroma ters  equipped  with 
a  2  meter  concave  grating  of  7500  lines  per  inch  and  an  entrance 
slit  of  width  50  microns.  The  angle  of  incidence  of  radiation  on 
the  grating  is  86  degrees.  The  narrowest  slit  employed  v.'as  0.75  mm 
wide,  with  foreshortened  width  of  53  microns  at  the  short  wavelength 
end  of  the  scan,  and  about  100  microns  at  the  long  wavelength  end. 

The  wavelength  increments  covered  by  the  slit  were  1.5  A  and  3  A  at 


41 


10 


50 


100 


PHOTON  ABSORPTION  RATE  DENSITIES,  o(X,AX,h)  IN  PH  CM*3SEC 


FIGURE  14,  Flux,  Extinction  Coefficients,  ond  Absorption  at  790  A 
(After  L.  A.  Hall  et  al . ,  1962) 
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PHOTON  ABSORPTION  RATE  DENSITIES,  o  (X  ,AX  ,h)IN  PH  CM'3SEC'' 

FIGURE  15 .  Flux,  Extinction  Coefficients,  and  Absorption  at  368.1  A 
(After  l.  A .  Hal  I  et  a  I . ,  1 962) 
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FIGURE  16.  Density  by  Atmospheric  Absorption  of  Solor  Radiation 

che  short  and  long  wavelength  limits,  respectively.  Dispersed  radia¬ 
tion  passing  t. '.rough  the  exit  sl-;t  in  any  position  was  intercepted 
by  a  photomultiplier  with  a  3.5  inch  long  sens i conductive  cathode. 

The  photomultiplier  output  was  input  to  a  pulse  amplifier  to  counting 
circuits,  with  overall  gain  so  adjusted  that  a  variation  in  gain  over 
a  factor  two  would  cause  a  change  in  counting  rate  of  only  5  percent. 
The  count  accumulation  was  read  out  continuously  ir  32  steps  in  the 
0  to  5  volt  d.c.  range  of  the  telemetry.  Readability  of  the  flight 
record  had  a  high  speed  upper  limit  of  two  megacycles.  However, 
since  the  chosen  entrance  slit  height  of  8  mm  limited  the  expected 
counting  rate  for  the  iie  II  resonance  line  to  about  5  (10)4  sec  , 
dead-time  corrections  were  not  required  for  most  cr  the  observed 
counting  rates,  with  the  exception  of  the  Hydrogen  Lyman  alpha  line 
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with  counting  race  5 .7  (10) J  sec  .  Hall  ct  al.  (1952)  estimate  that 
the  "estimated  error"  c;  the  absolute  photoelectric  yield  of  the 
tungsten  reference  electrode  is  ±10%  in  the  wavelength  range  536 
to  1216  A,  and  ±15%  tc  ±2551  in  the  range  250  to  536  A.  The  "estimated 
error"  in  quantitative  comparison  of  the  response  of  different 
rocket  instruments  is  much  smaller  than  the  errors  listed  above. 

An  alternative  to  the  grating  monochromator  and  photocell 
combination  described  above  is  the  use  of  narrow-band  photon- 
counters  or  ionization  chambers.  The  technique  of  preparing  such 
detectors  has  been  described  by  Chubb  and  Friedman  (1955)  and  by 
Grobecker  (1967). 

The  described  ion  chambers,  typically,  have  an  outside  diameter 
of  1  inch  and  are  1%  inches  long.  Each  chamber  has  an  end  window 
acting  to  filter  out  short  wavelength  radiation,  and  each  has  a 
gas  filling.  The  ionizing  potential  of  the  gas  filling  sets  the 
upper  limit  of  the  detector’s  'wavelength  response.  Windows  for 
narrov;-band  ion  chambers  in  the  soft  x-ray  region  are  ma.de  by 
coating  Mylar  film  with  sufficient  thicknesses  of  metal  and  using 
the  absorption  characteristics  of  the  filling  gas  as  part  of  the 
window  function.  For  windows  of  ion-chambers  in  the  ultraviolet 
region,  at  wavelengths  greater  than  1200  A,  solid  crystals,  such 
as  LiF,  3a F2 ,  etc.  are  used. 

The  fractional  uncertainty  of  measurements  of  flux  using  ion- 
chambers  is  about  the  same  as  those  made  using  monochromator  and 
photocell  combinations ,  approximately  ten  to  twenty-five  percent, 
as  obtained  by  rocket  borne  measurements  described  by  Grobecker  (1969). 

VEHICLE  AND  SYSTEM  CONSIDERATIONS 

The  orbit  characteristics  of  the  satellite  vehicle  for  the 
instruments  used  to  measure  the  absorbed  solar  radiation  may  require 
a  detailed  study  for  optimization  of  coverage,  spacing  of  profiles 
and  the  like.  A  notable  distinction  of  the  solar  absorption  measure¬ 
ments  from  the  atmospheric  emission  measurements  is  that  the  former 
may  be  accomplished  only  over  vertical  profiles  at  each  of  two  discrete 


positions  along  \ he  path  of  the-  Satellite,  as  compared  to  a  continuous 
swathe  u.orr7  Un-  orbit  of  the-  satellite-  maV  ip-j  omission  measurements . 

Of  course,  both  types  of  measurement  may  conceivably  be  mace  from 
the  same  satellite  vehicle;  however,  the  data  output  of  solar  absorp¬ 
tion  profiles,  not  continuous  along  the  orbit  path,  is  much  smaller 
than  that  of  atmospheric  emission  profiles. 

In  example  of  the  coverage  obtainable  is  illustrated  in  tig.  17. 
One  set  of  five  satellites  having  orbits  (numbers  1  tc  5  cf  the 
figure)  chosen  to  intersect  the  terminator  'with  a  large  angle  between 
the  orbit  path  and  the  terminator,  at  latitudes  evenly  spaced  between 
the  equator  and  the  pole,  are  launched  with  a  common  ascending  node. 
h  second  set  cf  four  satellites  (numbers  6  to  9  inclusive  in  the 
figure),  are  launched  in  orbits  having  similar  inclinations,  but 
with  a  common  ascending  node  displaced  90°  from  that  of  the  first 
ret.  If  the  satellite  altitude  is  about  150  nmi  (278  km),  the  period 
is  about  lb  hours.  Since  the  satellite  transits  the  terminator  twice 
in  oath  orbit,  each  satellite  measures  absorption  profiles  at  52  loca¬ 
tions  per  day.  Such  a  .onrination  o:  nine  satellites  yields  density- 

height  profiles  cf  one  observation  per  day  at  a  grid  spacing  or  about 
” 


Simple  geometric  considerations  lead  to  the  following  character¬ 
istics  of  the  volume  of  atmosphere  probed  by  the  technique.  The 
angle  of  the  earth's  horizon  below  the  horizontal  plane  of  a  satellite 
at  altitude  cf  200  km  is  about  3  degrees.  The  path  length  from  the 
satellite  to  the  point  of  earth  tangency  of  the  line  of  sight  is 
about  345  km..  The  distance  subtended  by  the  half  degree  diameter 
of  civ.  sun  at  a  distance  of  i-'S  V.m  is  about  2.75  km,  which  is  the 
approximate-  value  of  the  vertical  ar.d  meridional-horizontal  dimensions 
of  the  atmospheric  volume  observed.  The  dimension,  in  the  earth- 
sun  line,  of  the  observed  atmospheric  volume,  or  the  path  length 
within  one  atmospheric  scale  height  of  the  nominal  tangent  altitude, 
is  about  410  km.  In  short,  the  observations  determine  the  average 
density  in  a  volume  of  the  atmosphere  which  is  about  345  km  in  the 
direction  from  che  satellite  to  the  sun,  and  of  dimensions  about 
410  km  along  the  sun  lines,  and  2.75  km  normal  to  the  sun  line. 


OCCULT ATICK  CF  SATELLITE -TC—  SATELLITE  TRANSMISSION 

k  technique  for  one  determination  of  density  in  the  lower  tropo¬ 
sphere  is  the  coherent  microwave  occultation  system  described  by 
B.  Lusignar.  et  ai. ,  1256  and  1269,  along  the  .  :,es  demonstrated  by 
NASA  Mariner  spacecraft  in  the  measurements  of  the  atmospheres  of 
Mars  and  Venus  (A.  Xliore  st  al . ,  196 S  and  196"'), 

A  characterization  of  the  techniques  is  given  in  Fig.  16,  which 
summarizes  considerations  reviewed  in  the  following  paragraphs  dealing 
in  turn  with  the  basic  theory,  the  estimates  of  error  ana  system  con¬ 
siderations  . 

The  occultation  effects  have  beer,  desc: 'bed  by  Lusignan  et  al., 
1569,  in  terms  of  a  basic  geometry  illustrated  in  Fig.  15.  Two  satel¬ 
lites  in  the  same  orbit  are  spaced  so  that  the  radio  waves  from,  one 

to  the  other  pass  through  the  atmosphere .  In  the  atmosphere  the  radio 
waves  propagate  with  a  velocity  given  by 

v  “  c/r.  (5.1) 

where  c  is  the  speed  of  light  in  vacuum,  and  n ,  the  refractive  index 
is  given  by 

n  =  1  +  10“'  N  (3.2) 

The  refractive  modulus  N  is  proportional  to  the  air  density  and  water 
vapor  density,  as  described  in  Volume  II,  and  given  by  the  relation 

•<  pw 

N  =  0.223  0  4.  0.3^4  pw  +  1. 75(10)-'  (3.3) 

where  Pa  is  the  density  of  dry  air  in  err.  cm" ^ ,  p...  is  the  density  of 

water  vapor  and  T  is  the  temperature  in  degrees  Kelvin. 

The  changes  in  refractive  index  have  two  affects  cn  the  radio 
waves:  <3n  apparent  increase  in  path  length  due  tc  retardation  cr 
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PRINCIPLE  Of  OPERATION 

•  BENDING  AND  PHASE  CHANGE  OF  MICROWAVE  IN  OCCULTING  ATMOSPHERE  =  f(p,  T,  X) 

•  PATH  LENGTH  (X)  DETERMINED  BY  GEOMETRY 

|  •  DETERMINES  DENSITY  (p)  AND  TEMPERATURE  (T) 

P 

UNCERTAINTY  OF  MEASUREMENT  (Theoretical,  with  no  experimental  authentication) 

I  •  DENSITY  OF  AIR:  0.1% 

•  DENSITY  OF  WATER  VAPOR:  1% 

•  TEMPERATURE:  10% 

COVERAGE  IN  VOLUME 

:  •  AVERAGE  OVER  VOLUME  500  KM  x  1.3  KM  x  150  M 

!  •  SEVEN  SATELLITES  IN  GROUP  DETERMINE  DENSITY  AT  6  ALTITUDES  FROM  2.5  TO  20  KM 

•  TWENTY-ONE  SATELLITES  MAKE  DAILY  SAMPLE  ON  1000  KM  GRID 

STATUS  OF  DEVELOPMENT 

•  SYSTEM  SPECIFICATIONS  STUDIED;  ERRORS  ESTIMATED 

•  EVALUATION  OF  COMPUTATIONAL  REQUIREMENTS  INCOMPLETE 

RECOMMENDATIONS 

•  INVESTIGATE  FEASIBILITY  OF  A  7-SATELLITE  DEVELOPMENTAL  EXPERIMENT  TO  MEET  NEEDS  OF 

NUMERICAL  WEATHER  PREDICTION,  RADIO  COMMUNICATION  AND  BALLISTIC  WEAPON 
TARGETING.  THIS  MAY  BE  DONE  BY  CIVIL  AGENCIES. 


TRANS 


FIGURE  18.  Measurement  of  Microwave  Satel li te-to-Satel I ife  Occupation  for 
Density  at  4-20  km 


FIGURE  19.  Basic  Geometry  of  Occupation  Measurement 
(After  B.  Lusignari  et  al.,  1969) 


The  mathematics  of  the  inversion  have  been  described  in  detail 
by  G.  Fjelbo  and  V.  R.  Eshleman  (1958)  and  R.  A.  F'hinr.ey  and  D.  L. 
Anderson  (1969). 

Lusignat  et  al. (1969)  describe  a  computer  ray  tracing  program 
based  on  Fermat's  principle  that  the  path  followed  by  an  electro¬ 
magnetic  ray  through  a  refractive  medium  is  such  that  the  integral 
Juds  is  a  minimum.  Referring  to  Fig.  19,  the  integral  equation  for 
radius  r,  the  distance  from  the  earth’s  center  to  a  point  on  the  ray 
path,  as  a  function  of  the  angle  9  ,  measured  from  the  center  of 

T 

the  ray  path,  is  given  by 


f  'r(s) 

r(6 )  n(r(c  )) 

r  r. 

'o 

—  *  1  _ 

where  r  is  the  altitude  of  the  ray  at  the  tangent  point  %  •=  0. 

t  ~ 

Equation  3,4  is  integrated  by  computer  to  determine  the  actual  path. 
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The  integral  J'nds  along  the  path  ds  =  (r  d6  +  dr)  gives  the  total 
phase  path  length  of  the  ray. 

Lusignan  et  al.(1969)  suggest  that  the  experimental  behavior 
may  be  crudely  perceived  from  an  analytical  solution  assuming  spherical 
symmetry,  small  bending,  and  constant  scale  height  in  the  atmosphere. 
The  assumption  of  scale  height  impxies  that  the  density  p  is  an 
exponential  function  of  height  h: 


Neglecting  water  vapor  terms  in  Eq.  3.3,  the  refractive  modulus 

N  is 


N  -  Hse  ^  (3.6) 

where  the  subscript  s  implies  a  value  at  the  earth’s  surface. 

The  following  approximate  equations  f ollow  from  the  assumptions : 


Total  bending  angle 

e(h)  =  (2  x  10"6)Nse"h/H(rtR/2H)?l  (3.7) 

Total  excess  path  length 

0(h)  =  AR(h)  +  AL(h)  (3.8) 

Excess  path  length  due  to  retardation 

&R(h)  -  (2  x  10‘6)Nse"h/H(TrRH/2)^  =  He(h)  (3.9) 

Excess  path  length  due  to  bending 

lb(h)  =  R  e2(h)/ 4  (3.10) 

Cl. 

Difference  in  altitudes  of  closest  approach  of 
bent  and  unbent  rays 

h  -  h  =  Re (h)/2  (3.11) 

o  a 


where 


■;r  =  surface  refractivity 
R  =  distance  from  the  earth’s  center  to  the  ray 
H  =  atmosphere’s  scale  height 
h_  =  closest  approach  of  an  unbent  ray 

U  J 

h  =  closest  approach  of  the  bent  ray 

The  radio  system  of  the  satellites  measures  the  path  length 
within  about  one  half  wavelength  (or  3  cm).  A  continuous  radio  wave, 
transmitted  from  the  master  satellite  ■  '  the  repeater  satellite,  is 
retransmitted  phase  coherently  back  to  c’ne  master  where  the  received 
and  originally  transmitted  frequencies  are  compared.  A  difference 
in  transmitted  and  received  frequencies  is  a  "Doppler  frequency," 
recorded  on  a  counter  in  the  master.  Since  at  a  frequency  of  5  GHz, 
each  wavelength  is  6  cm  long,  the  system  has  an  uncertainty  of  one 
way  path  length  of  3  cm.  The  total  effective  path  length  change  is 
indicated  as  a  function  of  altitude  in  Fig.  20,  from  which  it  may 
be  inferred  that  the  uncertainty  of  3  cm  represents  a  fractional 

_  C  _  X 

uncertainty  of  5(10)  ~  near  the  surface  to  about  3(10)  "  at  20  km. 

Measurement  of  a  3  MHz  signal  modulating  the  5  GHz,  with  an 
uncertainty  in  phase  of  4  degrees,  yields  a  separate  indication  of 
one-way  path  length. 

In  passage  through  the  atmosphere,  the  radio  waves  are  influ¬ 
enced  by  the  air  in  a  volume  indicated  approximately  in  Fig.  21. 

The  cross  section  normal  to  the  ray  path  is  determined  by  the  Fresnel 
zone  which  in  the  horizontal  direction  has  a  dimension  of  about  1.3  km. 
In  the  vertical  direction  the  gradient  in  the  atmosphere,  contributing 
to  the  change  in  path  length,  shrinks  the  effective  Fresnel  zone  to 
about  150  meters.  Along  the  path,  the  atmospheric  effects  decrease 
mainly  as  a  result  of  the  radio  waves  leaving  the  atmosphere,  with 
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horizontal  distance  of  ii GO  km. 

As  indicated  roughly  by  tiro.  5.7  to  5.1i  inclusive ,  as  the 
atmospheric  pressure  increases,  the  total  path  length,  the  bending 
angle  and  the  altitude  of  closest  approach  also  increase.  Per  a 
given  satellite  spacing,  Lusignan's  (1169)  ray  tracing  program  is 
used  in  a  recursive  search  procedure  to  find  the  ray  that  reaches 
the  satellites.  Once  this  ray  is  identified,  its  starting  height 
is  recorded  and  phase  defect  calculated.  Figure  22  shoos  the  total 
phase  defect  as  a  function  of  surface  atmospheric  pressure  for  several 
satellite  spacing  distances.  For  the  illustration  of  Fig.  22,  the 
pressure  decays  with  altitude  according  to  the  scale  height  variation 
of  a  standard  atmosphere .  Figure  25  snows  the  height  of  closest  ap¬ 
pro  an  or  t.cc  sure  rays,  satoulitc  soacinos  u i .  ;  svr'ooo  measures. 

Lu  signer.  et  ai .  (1969)  have  analyzed  the  specifications  of  the 
system,  to  1  ■:  sue:-.  that  ur. ta  iut  its  arc  loss  than  about  C.101  ir. 
a';. solute  value  of  versify.  -c  yet,  cnere  has  beer.  riC  experimental 
authenti  -aticr.  of  Lucijv.nu  tutor  rial  error  estimates. 

Such  f  razticral  ur.ue:  tsinties  translate  tc  about  C.E  n  ir.  hori- 
sor.f a2  spacing  and  about  24  m.  ir.  -i  vert  1"- a]  direction.  The-  contribu¬ 
tions  to  horizontal  unuortaintc  arc  estimate  :  or.  theoretical  grounds 
by  I/u  sign  an  ot  al.  ( 19  Cr)  to  be-  the  following :  measurement  of  phase 
defect,  0.05  m;  ionesphe  r  it  ref  r  act  ion  r.egl  igible  :  radio  oscillator 
instability  negligible;  air  drag,  O.OC-  m,  solar  pressure ,  0.2  ir. ; 
orbit  oblateness ,  C*.2  m;  weather  f  luctuatior.s  averaged  over  ore  or 
two  r.cuths ,  less  than  0.S  r.. 

The  contributions  to  vertical  uncertainty  are  estimated  to  be : 
orbit  period  negligible;  orbit  oblateness ,  G.C2  m;  eccentricity  of 
orbit,  less  than  C.5  m;  and  oblateness  of  the  earth; ’s  surface,  less 
than  54  meters. 

The  oblateness  of  the  earth’s  surface,  known  with  uncertainty 
less  than  24  meters,  is  such  that  mean  sea  level  varies  about 
from  the  spherical  from  the  poles  to  the  equator.  Ir.  spacing 


20  km 


several  satellites,  the  oblateness  must  be  considered  so  that  the 
measurements  may  be  made  at  a  useful  set  of  altitudes.  Shown  in 
Fig.  24  is  the  variation  in  minimum  altitude  measured  from  sea  level 
of  each  of  a  set  of  six  satellites  which  pass  from  pole  to  equator. 

Set  1,  spaced  to  keep  its  rays  above  the  atmosphere,  suffers  the  full 
20  km  variation.  For  rays  in  the  atmosphere,  the  decreasing  altitude 
of  the  wider  spacing  is  offset  by  increased  bending  and  rise  in  ray 
altitude,  so  that  variations  from  pole  to  equator  for  sets  2  to  6 
are  only  5  to  .10  km. 

The  precise  tracking  information  is  derived  by  measurements 
between  the  individual  satellites,  not  from  ground  tracking  systems 
which  have  additional  sources  of  error.  However,  determination  of 
the  absolute  ground  projected  position  of  each  pair  of  satellites 
depends  upon  ground  tracking,  with  uncertainty  of  several  kilometers 
acceptable  for  the  needs  of  numerical  weather  prediction. 

The  effect  of  water  vapor  in  the  lower  troposphere,  given  analyt¬ 
ically  by  Eq.  3.3,  is  shown  in  the  example  of  Fig.  25.  The  total 
refractive  modulus  of  dry  air  is  .  'wn  as  function  of  altitude,  as 
well  as  that  of  an  atmosphere  with  the  partial  pressure,  of  'water 
vapor  at  the  surface,  of  10.2  m  bars.  While  water  droplets  in  clouds 
considerably  attenuate  radio  waves,  they  add  less  than  1  K  unit  to 
the  refractive  modulus  at  altitudes  2  km  and  above  and  consequently 
do  not  seriously  affect  the  phase  path  length. 

The  iter  vapor  content  of  the  atmosphere  may  be  sensed  by  a 
second  occultation  frequency  at  an  exact  multiple  of  the  primary 
5  GHz  signal  near  the  water  vapor  absorption  near  20  GHz.  Whether 
the  20  GHz  signal  has  a  different  refractive  index  and  consequently 
a  different  phase  path,  than  the  5  GHz  signal  is  currently  being 
evaluated . 

water  vapor  may  affect  the  occultation  system  by  offering 
multiple  paths  for  rays  between  satellites ,  with  scintillations 
strong  enough  to  disrupt  the  phase  lock  cf  the  satellite  receivers. 


t  is  expected  that  scintillation  errors  or  tr.is  kmc  will  affect 
r.iy  rays  in  the  lowest  2  km  of  the  atmosphere . 
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FIGURE  25.  Total  Refractivity  as  a  Function  of  Altitude 
( After  6.  Lusignan  et  al 1969) 


Radio  systems,  penetrating  the  clouds,  have  potential  for  sup¬ 
plying  meteorological  data  not  available  from  visual  and  infrared  sys¬ 
tems.  In  Fig.  26  is  shown  the  average  occurrence  of  clouds  at  various 
altitudes  between  latitudes  41°  and  77°  over  Russia  during  the  years 


in  summary,  a  pair  of  satellites  in  the  same  earth  orbit,  but 

Q 

separated  about  60  in  phase  will  have  a  radio  path  between  them  that 
intersects  the  atmosphere  continuously.  With  a  coherent  transponder 
between  satellites  at  a  frequency  of  5  GHz ,  changes  of  2  err.  in  the 
phase  path  can  be  recorded,  giving  theoretically  a  sensitivity  to 
changes  in  atmospheric  density  of  better  than  0.001  in  fractional 
uncertainty  of  density.  With  a  set  of  sever,  polar  orbiting  satellite 
to  sense  density  at  five  levels  in  an  altitude  range  2.5  to  20  km,  in 
each  of  three  systems,  samples  or  density  may  be  taker,  over  a  hori¬ 
zontal  worldwide  grid  of  spacing  less  than  1CGQ  km,  twice  daily. 


SATELLITE  DECELERATION 


The  effect  of  aerodynamic  drag  or.  the  orbit  characteristics  of 
satellites  becomes  important  for  military  applications  that  require 
orbits  with  a  low  perigee  altitude  (say,  100  nmi ) .  The  aerodynamic 
drag  at  satellite  velocities  at  such  altitudes  has  a  significant  ef¬ 
fect  on  the  life  cf  the  satellite  in  orbit,  which  can  expressed  as 
a  function  of  the  ballistic  parameter  of  the  satellite  ^ w/CpA ) ,  the 
orbit  eccentricity,  and  the  perigee  altitude.  For  a  typical  orbit 
eccentricity  of  5x10" ^  and  a  given  ballistic  parameter,  the  life  cf 
the  satellite  in  orbit  decreases  by  nearly  three  orders  of  magnitude 
as  the  perigee  altitude  drops  from  400  nmi  (i.e.  ,  TIROS  satellite)  to 
100  nmi.  This  drastic  effect  of  aerodynamic  drag  cn  the  life  of  satel¬ 
lites  at  low  altitudes  has  also  important  implications  when  determining 
the  position  of  a  satellite  in  low  altitude  erbit. 


Measurement  of  the  deceleration  of  the  satellite  by  means  of 
accelerometers  within  the  satellite  itself  is  a  technique  for  in  situ 
determination  of  atmospheric  density  at  satellite  altitudes.  Charac¬ 
terization  of  the  technique  is  given  in  Fig.  28  which  summarizes  con¬ 
siderations  discussed  in  the  following  paragraphs. 

The  aerodynamic  drag  force  on  the  surface  of  the  satellite  (D)  is 
a  function  of  the  ambient  density  (p),  the  satellite  velocity  relative 
to  the  ambient  medium  (V  ),  the  drag  coefficient  of  the  satellite  (Cn) 
and  the  frontal  area  of  the  satellite  (A),  i.e., 


D  =  \  p  Vf  C-A 


(4.1) 


where  \  p  denotes  the  dynamic  pressure,  and  is  the  drag 
coefficient  for  free -molecule  flow.  Since  the  ambient  density  at 
altitudes  above  80  nmi  depends  also  on  the  variability  of  solar 
phenomena,  the  accuracy  in  the  determination  of  the  life  and  position 
of  a  satellite  in  low  altitude  orbit  depends  cn  the  accuracy  in  the 


PRINCIPLE  OF  OPERATION 


•  DECELERATION  OF  SATELLITE  IN  ORBIT  =  f  |(P)a)n  ond  (v-  cd-  A^vehicleJ 

•  "DRAG  FREE  SATELLITE",  WITH  ACCELERATION  CONTROLLED  PROPULSION,  FLIES  IN  ORBIT 

nIPERTURBED  BY  DRAG 

•  VELOCITY  (V),  DRAG  COEFFICIENT  (Crf)  AND  DRAG  AREA  (A)  OF  VEHICLE  ARE  KNOWN 

•  DETERMINES  DENSITY  (p) 

UNCERTAINTY  OF  MEASUREMENT 

•  DENSITY:  10°b 

•  RELATIVE  DETERMINATIONS  HAVE  ERRORS  LESS  THAN  3% 

COVERAGE  IN  VOLUME 

•  DENSITY  MEASURED  BY  SATELLITE  IN  SITU 

STATUS  OF  DEVELOPMENT 

•  ACCELEROMETER  DEMONSTRATED  IN  OV1-I5,  OV1 -16  AND  OTHER  SATELLITES 
.  "'DRAG  FREE  SATELLITE""  PROPULSION  SERVO  DEMONSTRATED  IN  LABORATORY 

RECOMMENDATION 

•  DESIGN  AND  TEST  SATELLITE  INSTRUMENTED  WITH  ACCELEROMETERS  FOR  NEEDS  OF 

LOW  ALTITUDE  SATELLITE  ORBIT  PREDICTION,  USEFUL  FOR  R/V  DISCRIMINATION 
AND  RADIO  COMMUNICATIONS 


FIGURE  28.  Measurement  of  Satellite  Deceleration  for  Density  at  100-300  km 
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knowledge  of  the  variability  of  the  ambient  density  at  a  civer. 
altitude  above  80  nmi.  The  variability  of  the  at.bient  density 
such  altitudes  can  be  determined  from  measurements  of  the  drag  force, 
trie  satellite  speed,  the  crag  coefficient  and  the  frontal  area  of 
the  satellite.  The  accuracy  in  the  determination  of  each  of  these 
parameters  is  described  below. 

The  aerodynamic  drag  force  c.n  the  satellite  can  be  measured 
instantaneously  by  using  a  lev.’  -g  accelerometer  capable  of  detecting 
accelerations  of  the  order  of  1C  ^  g's.  A  requirement  for  such  low 
decelerations  levels  bring  about  a  need  to  calibrate  che  low  -g 
accelerometer  in  orbit.  This  has  been  accomplished  by  using  an 
accelerometer  that  is  supported  by  electrostatic  forces.  The 
accelerometer  is  calibrated  by  mounting  it  on  a  stable  platform 
provided  by  a  nen-spinning  satellite  such  as  the  Agena  vehicle.  The 
calibration  of  the  accelerometer  requires  determination  of  the  scale 
factor  and  bias  of  the  instrument  as  a  function  of  time  during  the 
measurement  of  the  aerodynamic  drag.  This  calibration  is  accomplished 
by  mounting  the  accelerometer  on  a  rotating  table  within  the  satellite. 
Since  the  angular  speed  of  the  rotating  table  can  be  varied  from  zero 
to  about  one  revolution  per  minute,  it  is  possible  tc  obtain 
accelerometer  data  (during  every  2C48  seconds)  for  trie  following  four 
conditions:  (1)  rotation  at  a  fixed  angular  speed  (of  Q .  44-495  rpm), 

(2)  rotation  at  a  fixed  angular  speed  equal  to  twice  the  first  value, 

(3)  no  rotation  'while  the  sensitive  axis  of  the  accelerometer  points 
forward  (towards  the  nose  of  the  satellite),  and  (4)  no  rotation, 
sensitive  axis  aft.  The  scale  factor  is  determined  from  the  average 
accelerometer  output  of  the  two  rotating  modes,  while  the  bias  is  fixed 
from  the  accelerometer  output  of  the  two  non-rctating  modes  together 
with  the  scale  factor.  Results  described  in  Appendix  E  indicate  good 
stability  of  the  bias  and  scale  factors  in  the  calibration  of  the 
accelerometer  over  67  orbits.  The  sensitivity  of  the  low  -g  acceler¬ 
ometer  to  10  g’s  provide  measurements  of  the  drag  acceleration  (F/M) 
at  low  orbit  altitudes  to  within  lie  (T=c3ra,  lech). 

The  satellite  mass,  needed  for  accelerometer  measurements,  is 
known  to  fractions  of  a  percent. 


The  speed  of  the  satellite  relative  to  the  ambient  medium1, 
depends  on  knowledge  of  the  motion  of  both  the  satellite  and  the  air. 
The  principal  uncertainty  in  this  relative  motion  is  due  to  the 
presence  of  winds,  which  may  bo  of  the  order  of  about  300  ft/sec  at 
very  low  satellite  altitude  (say,  70  r.mi).  Hence,  the  uncertainty 
.  .  2 

in  V  is  approximately  1%  and  that  in  V_  about  2%  (DeBra,  1969). 

O  25 

The  frontal  area  of  a  non -spinning  satellite  may  be  measured 
with  uncertainty  less  than  1%. 

In  terms  of  the  above  listed  uncertainties,  the  relative  error 
in  CDp  may  be  written 


A (CpP ) 


(-4.2) 


and  is  less  than  three  percent. 


The  drag  coefficient  of  a  satellite  in  free-molecule  flow  depends 
on  the  accommodation  coefficient,  the  distribution  of  the  reflected 
gas  molecules  after  impact  and  reemission  from  the  outer-most  layer 
on  the  surface  of  the  satellite  (i.e.,  points,  chemical  deposits,  etc), 
and  the  body  shape.  The  accommodauion  coefficient  (a)  indicates  the 
degree  of  energy  transfer  between  the  impinging  gas  molecules  and 
such  outer  layer  of  the  satellite  surface,  i.e., 


a  = 


(4.3) 


where  denotes  the  average  kinetic  energy  of  the  incident  gas  mole¬ 
cules,  E  that  of  the  re-emitted  gas  molecules,  and  E  the  kinetic 
’  r  b 

energy  corresponding  to  the  temperature  of  the  outer  layer  of  the 
satellite  surface.  The  accommodation  coefficient  also  depends, 
therefore,  on  the  ratio  of  the  mass  of  the  incident  gas  atoms  to  the 
mass  of  the  atoms  in  the  outer  layer  of  the  satellite  surface.  The 


distribution  of  the  reflected  gas  molecules  is  characterised  by  its 
two  extreme  mechanisms,  i.e.  ,  specular  and  diffuse  reflections.  Specu 
lar  reflection  occurs  when  the  angle  of  reflection  of  a  particle  is 
equal  to  its  angle  of  incidence,  while  diffuse  re-emission  takes 
place  according  to  the  Knudseri  cosine  law.  Finally,  the  body  shape 
does  r.ot  introduce  significant  changes  in  the  drag  coefficient  for 
rree-molecule  flow  as  the  average  energy  of  the  reflected  particles 
E  is  controlled  by  the  temperature  of  the  outer  surface  layer,  i.e,, 
as  E^  approaches  or  Q.  approaches  1.0. 

Theoretical  and  experimental  results  for  satellite  applications 
indicate  the  following  trends  concerning  the  parameters  controlling 
the  level  of  the  drag  coefficient  fer  rree-molecule  flow  (Cock,  1965): 
(1)  the  accommodation  coefficient  at  the  surface  of  most  satellites 
probably  exceeds  0.8  at  low  satellite  altitudes  (say,  below  200  nmi ) 
at  all  times  of  the  day  and  for  all  levels  of  solar  activity;  (2)  the 
exact  distribution  of  the  reflected  molecules  is  unknown,  and  it  is 
probably  best  to  assume  diffuse  re-emission  until  further  experi¬ 
mental  data  are  available.  Fortunately,  for  accommodation  coefficients 
near  unity,  the  drag  coefficient  is  not  particularly  sensitive  to  the 
distribution  of  the  reflected  molecules  and  little  error  will  ensue 
from  assuming  an  incorrect  distribution;  (3)  the  effect  of  body  shape 
on  the  drag  coefficient  for  free-molecule  flow  is  also  negligibly 
small  as  the  accommodation  coefficient  approaches  unity;  (4)  the 
effect  of  altitude  or.  the  drag  coefficient  for  free-molecule  flow  is 
negligibly  small  ft-  satellites  in  low  altitude  orbits  (say,  less 
than  200  nmi);  and  (S')  the  value  of  the  drag  coefficient  for  free- 
molecule  flow  used  in  recent  years  (i.e.,  Cc  =  2.2)  represents  a 
rood  estimate  based  on  current  knowledge,  but  it  could  be  low  by  as 
much  as  10  to  15  percent. 

Experience  with  U.  S.  satellites  in  determining  atmospheric 

drag  by  accelerometer  measurements  is  described  in  Appendix  B.  The 
results  or  071-13  and  CV1-1&  satellite  experiments  have  been  describe 
by  Champion  (1969),  Champion,  Marcos  and  Hclsaac  (1969)  and  Champion. 
Marcos  and  Schweinfurth  (1969). 


The  foregoing  considerations  indicate  that  the  variability  of 
the  ambient  density  at  low  satellite-  altitudes  can  at  present  be 
measured  with  an  uncertainty  nearly  equal  to  that  cf  the  drag  coef¬ 
ficient,  i.e . ,  to  witnin  an  r.n.s.  error  of  1G  to  IS  percent  (DeBra , 
1969) . 

The  satellites  0V1-1S  and  0V1-1S  for  determination  of  density  by 
directly  measuring  crag  decelerations  have  been  launched  in  polar 
orbits  of  perigee  altitude  about  200  km.  From  the  initial  altitude, 
drag  decay  of  the  small  eccentricity  orbit  has  resulted  in  perigee 
altitude  reduction  to  less  than  ISO  km  in  30  to  60  days.  Satellites 
in  highly  eccentric  orbits  or  with  large  mass  to  drag  area  ratios  may 
nave  lifetimes  of  a  year  or  more.  A  large  number  of  such  satellites 
would  be  required  to  give  coverage  required  for  use  systematically  as 
regular  inputs  for  purposes  of  computational  simulation. 

A  more  complicated  variant  is  the  "drag-free  satellite,"  which 
has  a  self-contained  propulsion  unit  servo  controlled  by  the  satellite 
accelerometer  to  maintain  the  orbit  unperturbed  by  drag  decay.  Life¬ 
times  of  a  year  for  a  "crag-free  satellite”  in  a  circular  orbit  may  be 
practical.  If  operated  in  a  circular  polar  orbit,  12  "drag-free  satel 
lites”  would  provide  data  for  one  altitude  continuously  along  the 
orbital  paths  separated  by  a  maximum  distance  of  about  1600  km.  Even 
two  satellites  in  polar  orbits  of  large  eccentricity,  90°  separation 
of  ascending  nodes  and  90c  differences  of  perigee  latitudes  will  yield 
useful  data  (Champion,  1969). 
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DUST  BACKSCATTER  OF  LASER  TRANSMISSIONS 


INTRODUCTION 

Recent  developments  of  laser-  technology,  and  the  exploration  of 
the  potential  of  the  technology  to  develop  atmospheric  probes  of 
density,  temperature,  aerosol  content  and  winds,  have  made  exciting 
strides.  Although  the  techniques  are  not  as  yet  fully  developed, 
there  is  a  basis  for  believing  they  have  distinctive  potential  capa¬ 
bilities  which  deserve  further  development. 

The  special  characteristics  of  laser  atmospheric  probes,  now 
commonly  called  iidars,  are  the  demonstrated  ability  to  detect 
profiles  cf  aerosol  layers  in  t atmosphere  up  to  altitudes  of  about 
IOC  km,  the  small  divergence  of  she  laser  beam,  (i.e.,  one  milli- 
radian)  which  can  probe  an  area  of  transverse  diameter  of  ICO  m  at 
100  km  distance,  and  the  small  pulse  duration  of  about  30  nano¬ 
seconds,  corresponding  to  a  line-of -sight  resolution  of  about  10  m. 

A  characterization  of  the  technique  is  given  in  Fig.  29  which 
summarizes  considerations  reviewed  in  the  following  paragraphs  dealing 
in  turn  with  basic  theory,  non-coherent  detection,  coherent  detection, 
uncertainty  cf  measurements  and  system  considerations. 

THEORY 

Analysis  of  the  laser  signals  returned  from  the  atmosphere  is 
■based  on  the  theory  of  Rayleigh  scattering  of  light  by  molecules, 

Mie  scattering  by  aerosols  and  on  the  theory  of  Raman  scattering. 

In  theory  in  a  dust-free  atmosphere,  Rayleigh  scattering  leads  to 
information  describing  density  and  temperature  profiles.  Goody 
(1964)  has  summarized  the  theory  for  Rayleigh  scattering  in  the 
atmosphere . 


PRINCIPLE  OP  OPERATION 


•  LIDAR  (A  PULSED  LASER,  ANALOGOUS  TO  A  RADAR )  ILLUMINATES  AND  RECEIVES  BACKSCATTER  FROM 

AEROSOLS  AND  MOLECULES.  ANALYSIS  IS  BENEFITED  BY  NARROW  BEAM  (0.4  MR),  NARROW  PULSE 
WIDTH  (20  NANO  SEC)  AND  NARROW  BANDWIDTH  ( 20A ).  AEROSOL  AND  MOLECULAR  SCATTER 
MAY  BE  DISTINGUISHED  BY  RAMAN  EFFECT. 

UNCERTAINTY  OF  MEASUREMENT 

•  BACKSCATTER  BY  AIR  AND  AEROSOLS:  15%  AT  20  KM;  100%  AT  40  KM 

•  BACKSCATTER  BY  HjO:  5%  AT  4  KM 

COVERAGE  IN  VOLUME 

•  BEAM  OF  0.4  MILLIRAD  YIELDS  RESOLUTION  OF  120  M  AT  300  KM  RANGE 

•  RANGE  RESOLUTION  1.5  KM 

•  SCATTER  VOLUME  IS  120  M  x  120  M  x  1500  M 

STATUS  OF  DEVELOPMENT 

•  DEMONSTRATIONS  OF  GROUND  BASED  LIDARS  TO  100  KM  RANGE  WITH  5  JOULES  IN  PULSE  OF 

5  4  SEC  SHOWS  PARTICLE  SCATTER  SUPERPOSED  ON  MOLECULAR  SCATTER 

•  DEMONSTRATION  OF  WATER  VAPOR  BACKSCATTER  AT  GROUND  FROM  ALTITUDES  0-4  KM 

•  RAMAN  SCATTER  OF  LIDAR  HAS  NOT  BEEN  DEMONSTRATED 

RECOMMENDATION 

•  RESEARCH  TECHNIQUES  TO  DISTINGU'SH  AEROSOL  (MIE)  SCATTER  FROM  MOLECULE  (RAYLEIGH) 

SCATTER,  INCLUDING  RAMAN  SCATTER  BY  MOLECULES 

•  DEVELOP  HIGH  POWER  LASERS 


FIGURE  29.  Measurement  of  Laser  Backscatter  by  02,  H2C>,  Oust  for  Density  at  0-30  km 
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The  scattered  irradiance  I  at  distance  R  is  related  to  the 
irradiance  I  of  the  incident  beam  of  linearly  polarized  (in 
direction  $  )  light  scattered  in  direction  e  by  an  optically  small 
particle  is  given  by 


In  the  scattering  of  linearly  polarized  light  by  optically  small 
particles,  for  which  the  condition  2  n  a  m'  «  \  is  satisfied, 
where  a  is  the  particle  radius,  m'  is  the  refractive  index  of  the 
particle,  and  >  is  the  vacuum  wavelength  of  the  radiation,  the 
angular  scattering  cross  section  may  be  written 


aCM)  =  sir** 

n3 


(5.2) 


where  m  is  the  refractive  index  of  the  gas  rather  than  of  the 
individual  particles  and  n  is  the  number  nsity  of  particles .  The 
total  cross  section  c  is  found  by  integrating  over  all  angles. 

For  air  under  standard  conditions  and  light  of  wavelength 
6328  A,  c(8j4)  is  2.6(10)"'®'  cm5.  Also,  the  scattering  efficiency 


Qs  *  tT?  ~3<10>"la 

is  very  small  for  the  molecules  found  in  normal  air. 


(5.3) 


In  scattering  by  spheres  which  are  not  small  compared  with  \  , 

o 

the  scattering  is  strongly  forward  directed,  may  have  several  lobes, 
and  varies  less  strongly  with  wave  length  than  in  the  Rayleigh  case. 
The  theory  of  such  scattering  has  been  worked  out  by  Kie,  so  that 


f  (4  n  a) 


(5.4) 


For  a  water  droplet  of  radius  1  micron,  the  scattering  efficiency 
is  a  maximum,  so  that  for  wavelength  6326  A,  the  total  cross  section 


-is  1.3  (1G)~7  cm2.  The  aerosol  particles  of  principal  importance 
in  atmospheric  optics  are  those  with  radii  0 .1  <  a  <  1.0  microns. 

In  this  range,  according  to  Goody  (1964) 

a  =  IT  a2  (11.7(1G)4  a  -  0.67)  :  0.1  <  a  <  0.4  microns 

0  -■  n  a2  (  -  3.3(10)4  a  +  5.3):  C.4  <  a  <  1 .0  p  .crons  (5.5) 


A  third  form  of  scattering,  well  known  in  physics  but  not  yet 
exploited  in  atmospheric  instrumentation,  is  the  Raman  scattering 
by  molecular  vibration,  which  is  uniquely  characteristic  of  the 
scatterer  and  which  radiates  a  scattered  frequency  different  from 
the  incident  radiation.  Raman  scattering  represents  displacement 
from  a  critical  line  due  to  changes  in  polarizability  of  the  mole¬ 
cule  during  vibration-rotation,  or  rotation  alone  or  vibration 
alone.  For  example,  for  N2  ,  the  Raman  scatter  radiation  is  dis¬ 
placed  to  a  wavelength  longer  than  that  of  the  exciting  radiation. 

Tf  the  latter  is  8583  A,  the  Raman  vibration  wavelength  is  1.128 
microns.  Other  displacements  for  NT2  are  due  to  interaction  with  the 
rotation  or  combined  rotation-vibration  mode,  and  are  much  fainter  in 
amplitude  to  be  observed  in  the  far  infrared  (Herzberg,  1964). 

If  Rayleigh  scattering  only  occurs  in  the  atmosphere,  tempera¬ 
ture  may  be  determined  from  the  temperature  gradient,  using  the 
relation  between  scale  height  (the  distance  in  which  density  varies 
by  factor  e)  and  temperature. 

Noncoherent  Detection 

Using  Equations  5.1  and  5.2,  the  intensity  on  the  ground  of 
a  pulse,  scattered  from  a  height  range  dh  at  altitude  h,  containing 
Nq  photons  in  a  interval  t  seconds  is  given  by 


I 


9 


a  n  NqT2 
t  h2 


dh 


(5.6) 


where  T  is  the  fraction  of  light  transmitted  by  the  lower  atmosphere. 
Allowing  for  the  finite  duration  t  of  the  pulse,  and  integrating 
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ever  a  time  t 
diameter  D  is 


the  number  of  photons  incident  cn  a  mirror  of 


n  tc 


(5.7) 


If  the  quantum  efficiency  of  counting  of  the  detector  is  e,  the 
number  of  photons  actually  counted  is  given  by 


n 

P 


TT  Da  TC  N 


(5.8) 


The  relation  (8)  gives  the  count  of  photons  to  be  returned  by 
Rayleigh  scattering,  and  is  the  basic  equation  by  which  the  per¬ 
formance  of  a  lidar  equipment  may  be  evaluated  . 

Equation  5.8  gives  the  parameters  to  be  considered  in  writing 
specifications  for  a  lidar.  Those  which  may  be  varied  include 
wavelength  the  diameter  D  of  the  receiving  mirror,  the  inte¬ 
gration  time  t,  the  number  of  photons  transmitted  N  ,  and  the 
quantum  efficiency  e  of  the  photomultiplier  tube  used  as  detector. 
The  transmission  T  of  the  atmosphere  must  also  be  taken  into  account 
since  it  is  wavelength  dependent. 

Choice  of  wavelength  should  consider  the  inverse  proportion  to 
trie  fourth  power  of  the  wavelength  of  the  scattering  coefficient 
o  and  that  the  quantum  efficiency  o  of  the  cetectcr  is  greater  for 
shorter  wavelengths .  However  in  the  lower  atmosphere  the  trans¬ 
mission  T  decreases  as  the  'wave  length  gets  shorter.  The  optimum 
wavelength  seems  to  be  between.  0.4  and  C  .  5  microns.  With  recent 
experiments,  however,  the  choice  of  wavelength  has  been  dictated 
by  the  fact  that  the  ruby  laser,  operating  at  0,634  microns,  is  the 
most  powerful  and  efficient  of  the  high  power  pulse  lasers  working 
in  the  visible  spectrum.  The  performance  of  lasers  as  transmitters 
is  currently  developing  rapidly. 


The  limitation  on  the  size  of  the  receiving  mirror  used  for 
non-coherent  detection  under  standards  easier  than  those  of 
astronomy  seems  only  to  be  one  of  cost. 

The  number  of  photons  transmitted  is  directly  proportional  to 
energy  output.  Clearly  the  highest  available  energy  should  be  used. 
In  general,  the  basic  lasing  pulse  is  0.5  to  1.0  millisecond.  Such 
length  does  not  yield  range  resolution.  A  system  of  Q-switching, 
involving  a  rotating  prism,  causes  lasing  to  occur  in  giant  pulses 
of  about  20  nsec  duration.  By  such  a  means,  a  series  of  giant 
pulses  is  obtained  which  have  a  total  energy  of  about  5  joules  in 
an  overall  duration  of  about  5  microsec  at  a  repetition  frequency 
of  about  0.16  pulses  per  second 

The  integration  time  for  a  single  pulse  depends  upon  the  height 
resolution  required.  A  long  integration  times  yields  a  larger 
number  of  photons,  but  decreases  the  height  resolution.  A  common 
integration  time  is  1C  microsec,  corresponding  to  height  resolution 
of  1.5  km . 

The  quantum  efficiency  of  photomultiplier  tubes  depends  on  the 
type  of  cathode  and  the  wavelength.  The  S-20  response  normally  used 
does  not  allow  high  efficiency  at  0.6S4  microns.  Typically  the 
efficiency  is  3  to  6%. 

Noise  in  an  optical  sounder  can  arise  in  two  ways:  due  to  an 
external  source  such  as  sky  noise,  cr  due  to  an  internal  source  such 
as  receiver  noise.  Sky  noise  is  usually  wide  band,  originating  from 
scattered  sunlight  in  day  time,  and  from  starlight,  airglow,  moon¬ 
light  and  other  minor  sources  at  night,  and  may  be  reduced  by  use  of 
a  narrow  band  optical  filter  in  the  receiving  system.  The  bandwidth 
of  the  transmitted  signal  is  approximately  ICT0  micron.  Because 
the  insertion  loss  of  optical  interf erence  filters  increases  rapidly 
as  bandwidth  decreased,  it  is  practical  to  use  a  filter  not  narrower 
than  2(10)-3  microns  (20  A).  Sky  noise  is  further  reduced  by 
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The  beamwidtn  of  the 


restricting  the  beamwidtn  of  the  laser  system, 
laser  itself,  about  10  millirad,  is  further  reduced  by  use  of  a 
parabolic  transmitting  mirror.  A  20  cm  mirror  reduces  effective 
Ifeamwidth  to  C .4  millirad. 

Noise  from  the  photomultiplier  is  reduced  by  cooling  of  the 
photocathode,  by  use  of  a  shutter  to  protect  the  photomultiplier 
from  the  fluorescence  of  the  laser  after  the  pulse  and  from  the  over¬ 
excitation  of  the  cathode  by  intense  light  scattered  by  the  lower 
atmosphere . 

Temperature  regulation  of  the  laser  is  required,  to  avoid  wave¬ 
length  drift  as  the  laser  heats  up  from  the  energy  dissipated  within 
it.  In  general  temperature  must  be  stabilized  to  be  within  a  band 
of  about  1C°C . 

Examples  of  lidars  in  current  use  have  been  described  by 
Clerr.esha  et  al. ,  1967,  Colli s  one  Linda,  1961,  Sancford,  1967,  and 
others  . 

Coherent  Detection 

The  coherence  of  electromagnetic  radiation  of  a  laser  suggests 
immediately  its  exploitation  in  ways  analogous  to  radar  Doppler 
measurement.  J,  C.  Owens  (1969)  has  made  a  system  analysis  of 
optical  heterodyne  measurement  of  Doppler  shifts  as  a  method  for  the 
remote  determination  of  vector  wind  velocity,  by  bistatic  measure¬ 
ments.  Owens  (1969)  technique,  since  it  is  bistatic  and  requires 
that  the  unscattered  radiation  travel  an  optical  path  of  equal 
length  to  that  of  the  scattered  radiation,  is  clearly  not  useful 
for  monostatic  measurements  from  satellites.  However  his  compu¬ 
tation  of  the  maximum  range  attainable  in  wind  or  temperature 
measurement  using  aerosol  or  molecular  scattering  is  interesting. 
Owens  finds  that  the  signal- to- noise  power  ratio  of  scattering  from 
scatterers  of  number  density  n  is  given  by 
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n  X3  j  (5i§) 


(5.9) 


S 

N 


hvE 


wherein  S  is  power  of  received  signal,  N  is  power  of  noise,  e  is  the 
quantum  efficiency  of  the  detector,  hv  is  the  energy  of  a  single 
photon,  PQ  is  the  laser  transmitter  power,  B  is  the  receiver  system 
bandwidth . 


X  is  the  wavelength  of  radiation  and  d  is  the  diameter  of  the 
spherical  scattering  volume.  Setting  the  signal-to-noise  ratio  to 
unity,  and  using  the  approximation  for  the  diameter  of  the  scatter¬ 
ing  volume  to  be 

d  *  2  X  R  (5.10) 

D 


where  R  is  the  range  from  transmitter  to  scattering  volume  and  D  is 
diameter  of  the  receiving  aperture,  Owens  (1969)  finds  that  the 
effective  range  is  given  by 
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(5.11) 


Assuming  the  values  of  the  system  parameters  to  be  as  follows: 


Detector  Efficiency  s  = 

Transmitter  Power  PQ  = 

Wavelength  X  = 

pe c 6 1 v in g  Ap 0  nt  u ro  D  — 

Diameter 


0  .05 
50  mw 
5328  A 
15  cm 


Owens  finds  that  ranges  of  such  a  system  may  be  the  following:  for 


clear  air  conditions,  using  bandwidth  E  =  50  kHz,  the  range  is  55 
cm;  and  for  fog,  using  bandwidth  of  825  Hz,  the  range  is  175  m. 
Because  the  thermal  motion  of  the  molecules  is  much,  greater  in 
frequency  broadening  than  the  displacement  of  the  'wind  induced 
Doppler,  the  range  calculated  by  equation  5.11  is  1.8(10)““  cm, 
which  is  to  say  the  measurement  of  wind  by  Doppler  is  impossible. 


For  similar  reasons ,  Owens  finds  that  measurements  of  temperature 
based  on  Rayleigh  scattering  are  also  infeasible. 

OBSERVATIONS 

Coyer  and  V’atscn  (1966)  have  summarized  recent  results  of 
observations  of  the  upper  atmosphere  by  lidar  techniques.  Figure  30 
shows  Goyer  and  Wat son Ts  summary  of  the  ratio  of  the  total  back- 
scatter  cross  sections  (molecular  plus  aerosol)  to  molecular  back- 
scatter  cross  section  as  a  function  of  altitude,  as  measured  by  a 
number  of  scientists  including  Bain  and  Sandford  (I960),  Clemesha 
et  aa.  (IS06 ,  1967}-,  Fioccc  and  Cclumtc  £lS64b) ,  Ficccc  and  Grams 
(1964a,  1966),  Fiocco  (1965),  Kent  et  al.  (1967),  Kishikcri  et  al. 

(1965),  and  Grams  and  Fiocco  (1967).  All  have  detected  the  18  to  20 

km,  scattering  layer,  and  a  broad  secondary  scattering  layer  around 
24  km.  Three  studies  (Fiocco  et  al.,  1963,  Bain  and  Sahcfcrd,  1SG6 

and  McCormick  et  al..  1966)  have  reported  scattering  layers  as  high 

as  120  km,  which  others  have  challenged  or,  uhe  grounds  that  the  poor 
signal  to  noise  ratio  at  the  high  altitudes  could  not  show  definite 
evidence. 

Schotland  (1965)  has  mace  use  of  the  6942. 15a  water  vapor 
absorption  line  in  an  attempt  co  determine  the  water  vapor  absorption 
profile  ir,  the  lower  4  km  of  the  atmosphere.  Since  the  half-width 
of  the  line  is  only  0.1  A,  the  measurement  required  a  shift  ir.  the 
wavelength  of  the  laser  radiation  of  only  0.6  A,  easily  obtained  by 
the  temperature  control  of  the  laser  rod.  By  measuring  the  ratio 
R(z)  of  the  return  signal  intensity  in  the  absorption  line  to  that 
outside  that  line,  at  discrete  altitudes,  Schotland  obtained 
qualitative  water  vapor  profiles  shown  in  Fig.  31 . 

Uncertainty  of  Measurements 

Collis  arc  Ligda  (1966)  have  estimated  the  statistical 
variability  in  groups  of  ten  lidar  observations  between  20  and  40  km 
to  be  the  rollowing,  in  terms  of  standard  deviation  divided  by .the 
mean:  0.15  at  20  km,  0.25  at  25  km,  0.5  at  33  km,  0.7  at  35  km,  and 
1.0  at  40  km. 


FIGURE  30.  Ratio  of  the  Total  Backscatter  Cross  Section  (  Molecular  Plus  Aerosol ) 
to  Molecular  Backscatter  Cross  Section  as  a  Function  of  Altitude 
(After  Go/er  and  Watson,  1968) 


Schotland  (1965)  has  estimated  the  uncertainty  of  the  data  on 
water  vapor  in  the  atmosphere  up  to  4  km,  shown  in  Fig.  31  to  be 
between  3  and  5%. 

It  is  reasonable  to  expect  that  the  uncertainty  of  measurements 
will  depend  on  signal  to  noise  ratio,  which  will  be  smaller  at  long 
range  than  at  short  range. 

-  SYSTEM  CONSIDERATIONS 

The  high  power  laser  has  been  proven  as  a  tool  for  probing  the 
atmosphere  for  total  particle  number  density  with  fractional 
uncertainty  of  about  0.15  and  water  vapor  density.  It  appears  to 
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FIGURE  31 .  Water  Vapor  Sounding  Data  for  26  Marc 
(  R.  M.  Schotland  et  al  „ ,  New  York  Un 
{  From  Goyer  and  Watson,  1969) 


have  potential  as  a  probe  for  determination  of  density  of  other 
atmospheric  components  anJ  in  particular  aerosols.  Chief  limitation 
on  use  of  the  technique ,  in  satellites  other  than  the  incomplete 
state  of  development  of  components ,  appears  to  be  the  inefficiency 
(high  dissipation)  of  input  power. 
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MOTION  OF  CLOUDS  FROM  SATELLITE  PICTURES 

The  motion  of  clouds  is  to  some  extent  a  measure  of  the  wind. 
Such  motion  may  be  deduced  from  pictures  of  clouds ,  made  from  cameras 
in  satellites. 

A  characterization  of  this  technique  is  given  in  Fig.  32,  which 
summarizes  considerations  reviewed  in  the  following  paragraphs. 


PRINCIPLES  OF  OPERATION 

•  MOTION  OF  CLOUDS  IN  TIME  SEQUENCE  PICTURES  FROM  SATELLITES  DETERMINES  WIND 

UNCERTAINTY  OF  MEASUREMENT 

•  WIND  VELOCITY:  25% 

•  LIMITED  BY  IDENTIFICATION  OF  CLOUDS  WHICH  MOVE  WITH  AIR 

COVERAGE  IN  VOLUME 

•  AT  CLOUD  HEIGHT  (  1-18  KM)  ONLY 

STATUS  OF  DEVELOPMENT 

•  DEMONSTRATED  BY  TIROS,  NIMBUS,  ATS 

RECOMMENDATION 

•  SINCE  TECHNIQUE  HAS  LIMITED  UTILITY,  LEAVE  DEVELOPMENT  TO  OPERATORS  AND  OTHER  AGENCIES 


FIGURE  32.  Measurement  of  Cloud  Motion  From  Satellite  Pictures  for  Wind  at  1  -18  km 


A  series  of  successive  pictures  obtained  from  satellites  in 
geosynchronous  orbit  (i.e.,  in  the  equatorial  plane  at  about  19,300 
nmi  altitude)  can  show  in  a  global  scale  a  rather  continuous  notion 
of  given  cloud  elements  over  geographical  locations  under  sunlight 
conditions.  The  continuous  motion  of  given  cloud  elements  is 
detectable  under  sunlight  conditions  because  geosynchronous  satel¬ 
lites  remain  fixed  relative  to  the  earth.  The  motion  of  clouds 
then  becomes  of  interest  because  they  can  be  used  as  tracers  for  the 
determination  of  the  horizontal  wind  field  at  mainly  low  altitudes, 
since  preliminary  results  indicate  that  the  population  of  traceable 
clouds  has  a  distribution  of  about  60  percent  below  5000  ft  (850  mb), 

30  percent  between  300  ana  200  mb,  and  10  percent  in  other  altitudes 
(Fujita,  1969).  Proposals  for  the  early  GARP  program  would  utilize 
two  kinds  of  high  resolution  instruments  for  the  observation  of  the 
cloud  elements:  (1)  an  optical  scanner  similar  to  the  spin-scan 
cameras  used  in  the  NASA  ATS-1  and  ATS- 3  satellites  (Appendix  A), 
which  would  provide  daytime  cloud  tracking  with  a  horizontal  resolu¬ 
tion  of  2  to  3  km  with  an  accuracy  for  horizontal  displacement  of 
about  1  km;  (2)  a  scanning  infrared  radiometer  operating  in  the 
8-12  p.  window  region  to  measure  the  altitude  of  cloud  tops  with  an 
accuracy  of  ±0.5  km  (assuming  a  temperature  profile  and  a  uniform, 
opaque  cloud),  and  a  horizontal  resolution  of  8  to  12  km  with  an 
accuracy  for  horizontal  displacement  of  about  3  km.  This  IR  radio¬ 
meter  would  provide  pictures  with  a  spatial  resolution  range  (at 
nadir)  of  about  8  x  S  to  16  x  16  km  (depending  on  the  optics)  of 
the  disc  of  the  earth,  with  a  full  picture  available  every  30  minutes. 
This  s  .nsor  subsystem  would  yield  daytime  wind  data  from  cloud 
elements  that  are  larger  than  8  to  12  km,  but  which  have  a  detectable 
structure  of  smaller  scale.  Since  the  nighttime  wind  data  will  not 
be  as  good,  this  subsystem  is  limited  to  daytime  and  cloudy  conditions. 

Measurements  of  horizontal  wind  profiles  from  cloud  motion  have 
been  made  using  available  pictures  from  the  ATS-1  and  ATS- 3  geo¬ 
synchronous  satellites.  The  spin-scan  subsystem  on  these  satellites 
provides  pictures  with  a  resolution  of  about  4  km  (Appendix  A).  The 


"cloud  velocity"  is  obtained  by  the  measurement  of  a  cloud  displace¬ 
ment  as  derived  from  a  sequence  of  ATS  photos  (Johnson,  1967).  The 
suitable  cloud  elements  are  selected  by  avoiding  clouds  that  do  not 
move  with  the  wind,  i.e.,  (1)  orographic  clcucs,  (2)  eumolonimbus 
clouds  (3)  clouds  that  are  related  to  the  phase  velocity  of  hori¬ 
zontal-transverse  wave  motions,  (4)  clouds  with  a  rapidly  changing 
structure  at  their  boundaries,  and  (5)  clouds  in  the  area  of  jet 
streams  where  the  windshear  is  rather  large.  The  cloud  displacements 
have  been  measured  by  using  a  stereoscope  and  a  plastic  measuring 
device  called  a  "stereo  wind  graph."  Two  pictures  properly  positioned 
under  the  stereoscope  allows  one  to  measure  the  change  in  cloud 
position  accurately.  The  reduction  of  cloud  displacement  to  cloud 
velocity  requires  the  following  considerations:  (a)  the  scale  of 
the  photographs  as  obtained  from  the  relative  dimension  of  the 
earth’s  picture  image  diameter  to  the  earth’s  diameter;  (b)  recti¬ 
fication  of  the  geometry,  since  one  is  looking  at  the  projection 
of  an  oblate  spheroid  on  the  flat  plane  of  an  ATS  picture;  (c)  cor¬ 
rections  for  aspect  ratio,  which  is  the  ratio  of  oblateness  introduced 
into  the  printed  ATS  picture  from  the  deflection  of  the  electronic 
beam  in  the  display  print  system.  The  picture  elements  are  "squeezed" 
together  so  that  the  picture  is  printed  oblately;  and  (d)  print  paper 
distortion  introduced  by  paper  stretching  during  the  developing 
process.  The  resulting  cloud  velocity  values  can  then  be  compared 
with  radiosonde  wind  measurements  by  assuming  that  the  clouds  move 
with  the  wind  field.  Since  the  altitude  of  the  clouds  in  the  ATS 
pictures  is  unknown,  the  dew  points  determined  from  the  radiosonde 
soundings  were  used  to  fix  the  levels  of  the  clouds. 

A  comparison  of  the  direction  and  magnitude  of  the  velocity 
vectors  for  the  cloud  motion  derived  from  ATS  pictures  with  radio¬ 
sonde  wind  data  for  61  observations  has  been  made  by  M.  H.  Johnson 
(1968)  who  indicates  the  following  results:  (1)  the  standard  devia¬ 
tion  between  the  direction  of  the  cloud  motion  and  the  radiosonde 
wind  data  is  15  deg  for  wind  velocities  that  are  higher  than  about 
4  m/sec  (see  Fig.  33);  and  (2)  the  standard  deviation  between  the 


magnitude  of  the  cloud  velocity  and  the  radiosonde  wind  data  is  2.2 
m  sec  (see  Fig.  34).  The  fractional  uncertainty  of  wind  estimation 
varies  from  about  0.25  at  15  ir.  sec  ^  to  1.0  at  3  ir.  sec  ^ .  Comparable 
errors  have  been  reported  by  L.  F.  Hubert  and  L.  Whitney  (1969)  and 
Fujita  et  al .  (1968). 
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FIGURE  33-  Distribution  of  Directional  Differences  Between  Observed  and 
Calculated  Wind  Directions  (After  M.  H.  Johnson,  1967) 

The  use  of  computational  techniques  for  determination  of  divergence 
and  rotation  of  a  cloud  field  promises  an  improvement  over  the  tech¬ 
niques  based  on  comparison  of  the  cloud  positions  with  terrestrial 
references  . 

The  foregoing  results  from  ATS  picture  data  indicate  that  obser¬ 
vations  from  geosynchronous  satellites  of  the  motion  of  selected 
clouds  during  a  short  time  can  yield  the  horizontal  wind  velocity  at 
lew  cloud  levels  over  given  geographical  locations  under  sunlight 
conditions .  The  horizontal  wind  field  derived  from  Geosynchronous 
satellite  picture  data  by  itself  lacks  the  accuracy  required  for  the 


CALCULATED  WIND  SPEED 


long-range  forecasting  of  the  weather.  However  the  wind  data  derived 
from  cloud  pictures  appear  to  be  valuable  as  one  form  of  initial  con¬ 
dition  input  to  the  computational  simulation  of  'wind  fields. 


FIGURE  34.  Distribution  of  Speed  Differences  Between  Observed  and 
Calculated  Wind  Speed  (After  M.  H.  Johnson,  1967) 


TRACKING  OF  BALLOONS  FROM  SATELLITES 


The  measurement  of  the  displacement  of  constant-level  satellite 
tracked  balloons  is  a  technique  for  the  direct  measurement  of  wind, 
which  seems  to  nave  great  promise. 

A  characterization  of  the  technique  is  given  in  Fig.  35  which 
summarizes  specifications  and  performance  to  be  developed  by  the 
Global  Atmospheric  Research  Program,  an  international  collaborative 
effort  which  is  receiving  considerable  support  within  the  United 
States.  The  detailed  description  of  the  the  program  has  been  de¬ 
scribed  in  detail  by  Bolin,  1967. 

Components  of  the  balloon  system  have  been  partially  demon¬ 
strated.  The  first  satellite  system  tests  are  scheduled  for  NIMBUS-D 
in  1970-71.  However,  since  the  technique  appears  to  be  receiving 
great  attention  from  civil  agencies  and  the  system  of  balloons  also 
seems  to  have  great  military  vulnerability,  it  is  recommended  that 
the  major  tasks  of  the  development  be  left  to  the  civil  and  inter¬ 
national  agencies. 


PRINCIPLE  OF  OPERATION 

•  BALLOONS  AT  CONSTANT  DENSITY  LEVEL  ARE  TRACERS  FOR  WINDS,  WITH  DISPLACEMENT  MEASURED 


FROM  SATELLITE  BY  DIFFERENTIAL  DOPPLER  OR  EOLE  RANGE  RATE  MEASUREMENTS. 
UNCERTAINTY  OF  MEASUREMENT 


LOCATION  TECHNIQUE 

LOCATION 

ERROR 

WIND  ERROR 

FRACTIONAL 
WIND  ERROR 

RANGE/RANGE-RATE  SYSTEM 

1-2  KM 

0.15-0.3  MSEC'1 

1.5  -  3% 

DIFFERENTIAL  DOPPLER 

10  KM 

2.5  MSEC*’  IN  90  MIN 

25% 

1.3  MSEC'  IN  2.4  HR 

13% 

COVERAGE  IN  VOLUME 

•  10,000  BALLOONS  AT  200  MB  (10-15  KM  ALTITUDE)  WITH  LIFETIME  1  YEAR  YIELD  ONE  MEASUREMENT 

PER  DAY  ON  400  KM  GRID  IN  ZONE  60°N  -  60°S 

•  5,000  BALLOONS  AT  800  MB  (I  KM  ALTITUDE)  WITH  LIFETIME  20  DAYS  YIELD  ONE  MEASUREMENT  PER 

DAY  ON  430  KM  GRID  IN  ZONE  30°N  -  30°S 

STATUS  OF  DEVELOPMENT 

•  BALLOON  DEMONSTRATED  BY  U.S.  AND  FRANCE 

•  LOCATION  SYSTEM  EOLE  DEMONSTRATED  IN  AIRCRAFT;  DIFFERENTIAL  DOPPLER  NOT  DEMONSTRATED 

•  BALLOON  ELECTRONICS  AND  POWER  SUPPLY  DEMONSTRATED  BY  EOLE 

•  FIRST  SATELLITE  TESTS  NIMBUS  D  1970-71 

RECOMMENDATION 

•  SINCE  MILITARILY  VULNERABLE,  LEAVE  MAJOR  DEVELOPMENT  TO  CIVIL  AND  INTERNATIONAL 

ORGANIZATIONS 


SATELLITE-TRACKED  BALLOONS 


FIGURE  35.  Measurement  of  Displacement  of  Satellite  Tracked  Balloons 
for  Wind  at  1  km  and  10-15  km 


NUMERICAL  METHODS  AMD  COMPUTER  SIMULATION 


A  review  of  techniques  for  measurements  from  satellites  indicates 
that  wind  may  be  determined  explicitly  by  utilization  of  balloon¬ 
tracking  satellites  together  with  satellite  TV-like  observation  of 
cloud  notions.  However,  these  techniques  appear  to  be  marginal  as 
'udgea  both  on  theoretical  grounds  and  from  the  results  of  past  per¬ 
formance.  An  alternative  to  such  explicit  measurement  techniques  is 
the  implicit  determination  of  winds  by  computational  simulation  of 
atmospheric  notions,  and  a  discussion  of  such  an  approach  forms  the 
content  of  the  present  section. 

One  application  of  computers  in  a  simple  form  is  the  use  of 
temperature  data  as  determined  by  TIROS  satellites  in  the  mid-latitude 
field  for  the  estimation  of  velocities  utilizing  the  geostropic  wind 
equation  as  described  in  Volume  IV.  The  error  in  using  the  approxi¬ 
mate  geostropic  wind  equation  results  in  uncertaintites  of  the  order 
of  20%  in  the  wind,  a  figure  which  is  adequate  for  many  purposes  (for 
example  in  airline  flight  routing),  but  in  others,  such  a  numerical 
weather  prediction,  the  geostropic  approximation  is  clearly  quite 
inadequate . 

A  further  important  exploitation  of  computation  simulation  is 
the  prediction  of  cloudiness.  In  view  of  the  importance  of  albedo 
variation  in  we a the r  processes,  the  variability  of  cloud  distribution 
is  a  most  significant  contributor  to  the  variability  of  winds  in  the 
troposphere.  Experiments  in  the  prediction  of  cloudiness  have  been 
undertaken  by  Kasahara  and  his  associates  at  NCAR,  Boulder,  Colorado. 

The  global  simulation  of  the  atmospheric  circulation  in  the  upper 
atmosphere  for  the  implicit  determination  of  winds  using  a  full  set 
of  atmospheric  input  data  has  been  suggested  by  Muryatroyd  (1968) . 

The  success  of  such  simulations  for  the  lower  atmosphere  (the  tropo¬ 
sphere)  in  the  last  few  years  lends  credence  to  the  belie:  that  this  , 


Cl 


much  larger  tack  (which  must  include  the  tropospheric  simulation  as 
part  of  the  total)  can  indeed  be  accomplished.  To  indicate  the  nature 
of  the  problems  to  be  faced  in  such  a  project  the  remainder  of  this 
section  will  be  devoted  to  a  brief  summary  of  the  basis  of  the  numer¬ 
ical  programs  that  are  used  today  for  global  weather  and  climate 
simulations,  together  with  an  indication  of  the  types  of  difficulty 
that  have  beer,  overcome  in  their  implementation,  and  will  conclude 
with  an  estimate  of  the  scope  cr  the  corresponding  numerical  problem' 
for  simulation  of  the  complete  atmosphere  from  0  to  200  km. 

A  characterisation  of  the  technique  of  computational  simulation 
of  the  atmosphere  as  a  means  cf  determining  wind  at  all  altitudes  is 
given  in  Fig.  36,  -which  summarizes  considerations  discussed  in  the 
following  paragraphs  dealing  in  turn  with  differential  equations  for 
tropospheric  simulation,  finite  difference  forms  and  stability, 
nonlinear  effects  and  stability,  boundary  conditions,  characertistics 
of  current  tropospheric  numerical  simulations,  and  the  extension  from 
current  tropospheric  simulations  tc  simulations  of  the  global  upper 
atmospheric  circulation. 

Differential  Equations  for  Tropospheric  Circulation 

Until  the  advent  of  high  speed  computing  techniques,  attempts 
to  solve  the  mathematical  problems  of  atmospheric  motions  were  based 
on  gross  simplifications  cf  the  governing  equations,  several  of  'which 
have  been  discussed  in  "Hydrodynamic  Theory  of  Wave  Motion  in  the 
Troposphere"  of  Volume  TV.  Today,  with  the  aid  of  rapid  computation 
capability,  a  much  more  complete  and  satisfactory  set  of  fluid  dynamic 
equations  can  be  tackled.  The  commonly  accepted  formulation  that  is 
believed  to  be  quite  adequate  for  describing  tropospheric  conditions 
(i.e.,  weather)  consists  of  the  compressible  Navier-Stokes  equations 
with  the  following  assumptions: 

a.  The  planetary  atmospheric  layer  is  thin  compared 
with  the  earth’s  radius. 

b.  The  horizontal  component  of  the  Coriolis  force 
is  negligible. 

c.  The  hydrostatic  equation  holds. 


PRINCIPLE  OF  OPERATION 

•  CONTINUITY  EQUATIONS  (  MOMENTUM,  MASS,  ENERGY)  ARE  BASIS  FOR  WINDS  =  f(p,  T,  ENEkGY, 

POSITION)  IN  TROPOSPHERE,  AND  MORE  VARIABLES  IN  UPPER  ATMOSPHERE 

•  ALL  REQUIRED  INPUTS  ARE  DETERMINABLE  BY  SATELLITE  ON  300  KM,  3  HR  GRID 

UNCERTAINTY  OF  DETERMINATION 

•  WIND:  25%  (GEOSTROPHIC  APPROXIMATION  IN  MID  LATITUDE  TROPOSPHERE);  « 25%  (PRIMITIVE 

EQUATIONS)  FOR  TROPOSPHERE;  LARGER  FOR  UPPER  ATMOSPHERE 

COVERAGE  IN  VOLUME 

•  GLOBAL,  OVER  ALTITUDE  RANGE  0-300  KM  ON  300  KM  GRID  AND  IN  3  HOUR  TIME  STEPS 
STATUS  OF  DEVELOPMENT 

•  DEMONSTRATED  FOR  TROPOSPHERE  BUT  WITH  INADEQUATE  DATA  INPUT 

•  TROPICAL  TROPOSPHERE  ANALYSIS  LIMITED  BY  INADEQUATE  INPUT  DATA  ( ESPECIALLY  WATER  VAPOR) 

•  UPPER  ATMOSPHERE  ANALYSIS  LIMITED  BY  INADEQUATE  INPUT  DATA,  KNOWLEDGE  OF  PROCESSES 

AND  PRESENT-DAY  COMPUTER  CAPACITY 

RECOMMENDATION 

•  DEVELOP  SYSTEMS  TO  DETERMINE  ALL  MEASUREABLE  ATMOSPHERIC  PARAMETERS  FROM  SATELLITES 

«  RESEARCH  ATMOSPHERIC  PROCESSES,  ESPECIALLY  IN  UPPER  ATMOSPHERE,  FOR  WHICH  1973  GARP  MAY 
PROVIDE  LOWER  BOUNDARY  J 

•  DEVELOP  LARGE  CAPACITY  COMPUTERS,  I.E.,  ILLIAC  IV  f 

•  BEGIN  UPPER  ATMOSPHERE  COMPUTATIONAL  SIMULATION 


COMPUTER  SIMULATION 


FIGURE  36.  Computational  Simulation  of  Atmosphere  for  Wind  at  All  Altitudes 


In  the  meteorological  community  the  Wavier- Stokes  equations  with 
these  assumptions  are  commonly  called  the  "primitive  equations." 

Some  treatments  of  atmospheric  circulation  problems  do  omit,  some 
but  not  all,  of  the  horizontal  Coriolis  force  terms  (assumption  b) 
in  the  equation^,  but  the  logical  consistency  of  this  is  questionable 
(Eckart,  1960). 


Other  assumptions,  affecting  the  "forcing  terms"  in  the  equations, 
i.e.,  those  involving  body  forces  and  energy  transfers  arising  from 
viscous  and  external  sources,  are  also  made,  but  these  can  be  dis¬ 
cussed  independently  of  the  general  formulation  of  the  partial  dif¬ 
ferential  equations.  These  equations  are,  under  assumptions  8.1 
(compare  "Physics  of  the  Troposphere"  of  Volume  IV). 
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Hydrostatic: 


where  x,  y,  z  are  rectangular  coordinates  in  which  the  (x,  y) 
plane  is  locally  horizontal  to  the  earth’s  surface. 

v  =  (u,v,w)  is  the  fluid  velocity 

s  is  the  entropy 

is  the  density 

is  the  temperature 

is  the  pressure 


o 


P 

T 

P 


F  =  (Fx,Fy,Fz) 

q 

g 

R 

CP 

Y 


is  the  body  force 
is  the  rate  of  energy  input 
is  the  force  of  gravity 
is  the  gas  constant 

is  the  specific  heat  at  constant  pressure,  and 
is  the  specific  heat  ratio 


In  the  normal  procedure  for  a  finite  difference  numerical  solu¬ 
tion  of  these  equations  the  first  four  are  treated  as  "prognostic" 
equations,  i.e.,  the  time  derivatives  are  evaluated  in  terms  of  given 
(observed)  initial  values  at  time  t  of  all  the  dependent  variables, 
and  then  the  equations  are  integrated  over  a  short  time  step  At.  The 
result  are  estimates  of  p,  u,  v,  and  s  at  time  (t  +  At)  in  each  of 
the  horizontal  planes  selected  for  the  finite  difference  vertical 
grid.  Involved  in  this  process  is  knowledge  of  the  vertical  velocity, 
w,  not  normally  obtainable  from  meteorological  observations;  in  fact, 
w  can  be  derived  from  the  remaining  two  equations,  8.0  and  8.7,  com¬ 
bined  with  the  continuity  equation.  The  result  is  the  Richardson 
equation  that  couples  the  vertical  to  the  horizontal  motions: 


where  is  the  "top"  of  the  atmosphere,  the  upper  altitude  limit  of 
the  computational  grid  system. 

The  details  of  how  the  computation  proceeds  between  the  horizontal 
prognostic  estimations  (Eqs.  8.2-6.S)  and  the  "diagnostic"  equations 
(S.6-3.8)  vary  with  different  numerical  schemes,  and  need  net  be 
entered  into  bore,  since  no  fundamental  questions  involving  computa¬ 
tional  characteristics  are  involved.  However,  it  may  be  m.cr.tioned 
that  a  common  variation  is  to  use  the  pressure  p  instead  of  z  as 


independent  vertical  coordinate,  with  a  subsequent  simplification  in 
the  diagnostic  Eq.  8.8.  Another  standard  modification  of  the  equations 
as  presented  here  is  the  use  of  longitude,  X,  and  latitude,  9,  in  place 
cf  x  and  y  as  independent  variables.  This  does  not  affect  the  general 
discussion  in  any  way. 

The  procedure  for  solution  described  above  requires  specification 
of  the  stress  terms  Fx,  Fy,  and  the  energy  source  q.  In  Volume  IV, 
"Hydrodynamic  Theory  of  Wave  Motion  in  the  Atmosphere"  discusses  the 
external  sources  of  energy  in  some  detail.  Indeed  their  accurate 
evaluation  is  a  very  important  and  difficult  part  of  the  total  fore¬ 
casting  problem;  we  cite  in  particular  the  questions  cf  radiative 
transfer  and  latent  heat  in  phase  changes  of  water  as  being  particu¬ 
larly  significant.  However,  the  mathematical  aspects  which  we  dis¬ 
cuss  here  are  not  sensitive  to  these  questions,  rather  they  are  in¬ 
fluenced  by  the  viscous  and  thermal  conductivity  contributions  to  Fx, 
Fy,  and  q,  since  these  are  functions  of  the  dynamic  state  of  the 
atmosphere.  In  the  normal  formulation  of  the  Navier-Stokes  equations, 
these  factors  are  given  as 

(Fx,Fy,Fz)  =  ~  ||i  v:>  v  +  3  (V  •  v)}  (S.9) 

q  =  -1  |k  V  '(VT)  +  (6.10) 

where  u  and  k  are  the  (molecular)  viscoscity  and  thermal  conductivity, 
respectively,  assumed  constant;  the  first  term  in  8.10  is  the  thermal 
diffusion  contribution  to  q,  and  $  the  viscous  dissipation  function, 
is  given 
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The  role  of  the  quantities  8.9  and  8.10  in  atmospheric  circu¬ 
lation  is  not  clear.  They  are  certainly  important  in  the  Ekman  layer, 
but  in  the  finite  difference  computational  schemes  the  latter  is  a 
small  fraction  of  the  vertical  grid  size  and  the  molecular  diffusion 
and  dissipation  properties  are  incorporated  via  lower  altitude  boundary 
conditions,  which  have  not  been  discussed  as  yet.  For  similar  reasons 
(i.e. ,  the  large  dimensions  of  the  finite  difference  mesh,  both 
vertically  and  horizontally),  it  has  been  argued  that  the  forms  8.9 
and  8.10  are  not  appropriate  for  the  finite  difference  solution,  and 
that  the  dissipation  and  diffusion  effects  should  rather  be  expressed 
in  terms  of  an  assumed  random  turbulent  flow  eddy  motion  occurring  at 
a  sub-grid  scale.  3y  analogy  with  turbulent  boundary  flow  theory, 
it  is  assumed  that  this  turbulence  can  be  represented  in  terms  of  a 
Reynolds  stress  and  an  effective  eddy  viscosity  coefficient,  usually 
with  much  simplified  expressions  for  F  and  q  as  functions  of  the 
velocity  or  the  velocity  derivatives. 

Whether  such  a  representation  is  valid  physically  is  not  known, 
however,  from  a  computational  point  of  view  the  introduction  of 
dissipation  can  be  vital  in  order  to  achieve  computational  stability 
as  will  be  discussed  below.  Apart  from,  this  important  contribution, 
the  effects  of  the  diffusive  and  dissipative  aspects  of  F  and  q  are 
probably  not  important  for  short  period  forecasting.  For  longer  time 
periods  (currently  estimated  as  three  days  or  more),  it  is  not  yet 
known  whether  a  different  representation  is  necessary,  or  even  whether 
a  stochastic  formulation  for  the  sub-grid  scale  effects  is  possible 
at  all.  This  question  will  be  entered  into  in  more  detail  later  in 
this  section  when  nonlinear  effects  are  considered. 

Finite  Difference  Formulations  and  Their  Stability 

The  prognostic  system  of  Eqs.  8.2  -  8.5  appears  at  casual  in¬ 
spection  to  be  straightforward,  if  lengthy,  to  solve  by  finite  dif¬ 
ferences.  Computation  times,  however,  are  considerable,  as  can  be 
seen  from  the  characteristics  of  prcsent-cay  numerical  programs  listed 
in  Table  S,  and  result  in  a  limitation  cf  their  utility.  The  difficulty 


arises  both  from  the  large  number  of  spatial  mesh  points  required  by 
a  global  model,  and  from  the  restriction  in  the  finite  difference 
representation  to  a  maximum  time  grid  interval,  as  required  for  com¬ 
putational  stability.  In  fact,  the  question  of  computational  stability 
is  a  principal  issue  in  the  effective  employment  of  numerical  simula¬ 
tion  of  the  atmosphere  and  a  brief  description  of  the  nature  of  the 
problem  follows . 

The  fluid  dynamic  equations  for  the  horizontal  motion  8.2  -  8.5 
describe  a  nonlinear  system  in  which  the  nonlinearities  are  introduced 
via  the  advective  terms,  ~-v,  and  the  dissipation  function  §.  Since 
there  are  no  known  methods  of  establishing  the  stability  or  convergence 
properties  of  finite  difference  representations  of  nonlinear  equations 
of  this  degree  of  complexity,  the  main  approach  for  treating  them  has 
had  to  be  heuristic.  Thus,  linearized  approximations  are  studied,  and 
time  and  distance  mesh  scales  are  selected  that  ensure  computational 
stability  for  those  linearizations.  The  resulting  numerical  programs 
are  then  checked  as  to  their  nonlinear  stability  performance  by  trial 
computations .  If  as  is  common,  the  trial  shows  an  instability,  damping 
is  empirically  introduced  into  the  system  of  equations,  either  with 
artificial  viscosity  terms  or  by  an  increase  in  magnitude  of  the  vos- 
cosity  that  may  already  be  present.  As  discussed  earlier,  it  is  ex¬ 
pected  that  the  frictional  effects  thus  introduced  should  have  little 
effect  on  scales  larger  than  the  grid  size,  so  that  the  approach  can 
be  regarded  as  a  device  to  alter  the  mathematical  properties  of  the 
finite  difference  scheme  without  altering  their  physically  important 
properties.  There  is  clearly  the  danger  in  that  ad  hoc  approach  that 
in  fact  the  mathematical  and  physical  characteristics  are  more  closely 
linked  than  assumed,  and  this  possibility  will  be  discussed  in  greater 
detail  later. 

The  rationale  for  assuming  that  the  properties  of  the  linearized 
equations  car.  be  used  to  approximate  the  behavior  of  the  nonlinear 
equations  is  based  on  the  consideration  of  local  behavior.  Thus,  if 
the  correct  solution  is  continuous,  and  we  assume  that  the  difference 
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TABLE  8.  GLOBAL  TROPOSPHERIC  CIRCULATION  SIMULATIONS 


equations  are  solved  by  a  progressive  marching  method,  as  outlined 
earlier  in  this  section,  we  might  expect  a  local  linearization  to  have 
global  relevance  when  the  finite  difference  mesh  is  sufficiently  fine. 
Continuity  of  the  true  solution  is  in  fact  ensured  by  inclusion  of  the 
diffusive  effects  of  viscosity  and  thermal  conductivity. 

A  discussion  of  linear  stability  can  be  based  on  the  wave  nature 
of  the  linearized  equations  for  the  horizontal  motion  as  described  in 
"Hydrodynamic  Theory  of  Wave  Motions  in  the  Troposphere"  of  Volume  IV. 

Atmospheric  wave  motions  are  classifiable  into  these  categories; 

1.  Acoustic  waves 

2.  Gravity  waves  (8.12) 

3.  Inertial  waves 

Actually  acoustic  waves  have  been  eliminated  in  the  formulation  8.2  - 
8.7  of  the  equations,  by  virtue  of  the  hydrostatic  assumption,  and 
the  irrelevance  of  sound  waves  to  atmospheric  circulation  is  well 
established,  at  least  for  the  troposphere.  The  main  contributor  to 
the  atmospheric  circulation  are  the  inertial  waves,  and  in  fact,  in 
the  early  approximations  of  atmospheric  circulation  only  the  inertial 
waves  were  included.  The  advent  of  greater  computer  capacity  allowed 
inclusion  of  both  gravitational  and  inertial  nodes,  as  in  Ecs.  8.2  - 
8.7. 


It  was  shewn  by  Courant,  Friedrichs  and  Levy  in  1928  that  a 
forward  differencing  scheme,  i.e.,  one  for  which,  in  one  dimension, 
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is  stable  (and  convergent)  as  long  as 


■where  c  is 


the  wave  velocity,  8.14  shows  that  for  a  chosen  oriel  length 
ix ,  stability  requires  that  the  time  step  be  chosen  to  be  less  than 
ix/c,  so  that  the  higher  the  wave  velocities,  the  shorter  is  the  time 
step  that  has  to  be  used.  The  velocities  for  the  three  waves  in 
Eg.  8,12  have  maximum  values  in  the  troposphere  of  about  300  m/sec 
for  the  acoustic  node,  150  m/sec  for  the  gravity  wave  and  10  m/sec 
for  the  inertial  wave.  The  corresponding  maximum,  allowable  time  steps 
given  by  Eq.  8.14  are  15  min,  30  min,  and  8  hr,  respectively.  The 
actual  figures  fer  the  full  linear  fluid  dynamic  equations  are  reduced 
from  these  values  by  two- dimensional  effects  and  by  advection.  The 
latter  results  in  8.14  being  modified  to 

(|v|  +  c)  -g  £  1  (8.15) 

where  |v |  is  the  magnitude  of  the  flow  velocity  at  the  point  in 
question.  The  greatest  reduction  in  bt  that  results  in  using  8.15  in 
place  of  8.14  occurs  for  the  inertial  wave,  since  its  velocity  is 
least. 

In  addition  to  a  stability  criterion  involving  the  wave  nature 
of  the  equations,  a  stability  criterion  can  also  be  formulated  for  the 
diffusive  part,  i.e. ,  for  the  viscous  and  thermal  conductivity  terms. 
The  linear  stability  condition  for  this  case  is 

3  — ~r  £  1  (8.16) 
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in  which  c  is  a  diffusion  coefficient.  The  finite  difference  scheme 
is  again  of  the  forward  type,  as  described  by  Eq.  8.13.  In  applica¬ 
tion  8.16  is  usually  less  severe  than  S.15,  so  that  diffusive  terms 
do  not  control  the  choice  of  grid  size.  An  example  of  the  use  cf  8.16 
is  considered  for  the  complete  atmospheric  model  (discussed  later), 
not  just  for  the  troposphere. 

The  primitive  Eqs.  8.2  -  8.7  exclude  acoustic  -waves  so  .that  the 
controlling  factor  in  determining  the  choice  of  time  step  is  the 
gravity  -wave  speed.  The  time  steps  listed  in  Table  8  for  the  current 


atmospheric  circulation  models  have  been  selected  or.  the  basis  of  a 
gravity  wave  stability  requirement. 

There  are  forms  of  finite  difference  approximations  that  can  re¬ 
duce  or  eliminate  linear  instability.  Backward  and  central  differ¬ 
encing  methods  when  applied  to  wave  and  diffusion  equations  fall  in 
this  category.  However,  such  systems  result  ir.  implicit  rather  than 
explicit  difference  equations,  so  that  a  simple  marching  process  for 
carrying  through  the  finite  difference  solution  is  no  longer  possible, 
and  solution  of  sots  of  simultaneous  equations  is  required.  This  is 
usually  regarded  as  sufficiently  time  consuming  that  explicit  equations 
are  preferred,  even  at  the  cost  of  the  smaller  time  steps  that  are 
needed  for  stability  of  the  latter. 

The  numerical!  simulations  listed  in  Table  3  employ  a  variety  of 
differencing  techniques ,  the  choices  mace  involving  compromises  be¬ 
tween  stability  requirements,  truncation  error  and  computation  time 
for  the  problem  of  interest.  All  these  approaches  car.  be  analyzed  in 
terms  of  linear  approximati ur.s  for  stability,  and  result  in  criteria 
exemplified  here  by  3.15  and  3.16.  Nonlinear  instabilities  can  also 
influence  the  type  c-f  finite  difference  equations  used;  the  factors 
involved  here  are  strongly  linked  to  the  sub-grid  scale  effects  cf 
turbulence  discussed  earlier,  and  the  next  paragraph  'will  discuss 
these  jointly. 

Nonlinear  Effects  and  Stability 

The  first  attempts  to  solve  the  primitive  equations  showed  up 
instabilities  of  an  unanticipated  variety,  net  predicted  by  the 
linearized  equations.  A  possible  cause  for  their  appearance  was 
suggested  by  Phillips  (1959),  and  explained  by  him  in  terms  of  a 
simple  nonlinear  model  equation,  it  being  assumed  that  a  similar 
result  would  apply  also  to  the  complete  equations.  The  essence  of 
Phillips  argument  is  that  nonlinear  mixing  (i.e.  ,  the  product)  cf 
two  different  frequency  -wave  components  of  the  solution  leads  to 
beat  frequencies;  if  the  wavelength  of  one  of  these  should  happen 
to  be  less  than  the  grid  size  (sc  that  it  is  r.ot  resolved  by  the 
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this  v.'ill  be  misinterpreted  by  the  finite 
s  wave  cf  length  greater  than  the  grid  interval 
elution,  The  effect  cf  finite  differencing  on 
non  is  the:  tc  introduce  specious  -wave  codes, 
is  implies  the  raise  introduction  of  energy 
old  in  fact  be  unresolved  within  the  sub-grid 
was  further  able  to  shew  that  the  effect  could 
lead  tc  instabilities  that  could  be  slowed 
by  reducing  the  grid  sise. 


A  method  employed  by  Phillips,  and  also  by  several  of  the  numerical 
schemes  summarised  in  Table  6,  for  controlling  this  form  cf  computa¬ 
tional  instability  is  to  add  smoothing  to  the  system  by  the  introduction 
of  friction  terns  and  dissipation.  The  latter  are  represented  in  the 
nonlinear  fluid  dynamic  Eqs.  3.2  -  3.5  by  the  terms  F  and  ;.  As  dis¬ 
cussed  earlier,  the  magnitudes  assigned  tc  these  quantities  have,  in 
many  instances  of  atmospheric  modeling,  been  chosen,  not  by  physical 
arguments  concerning  the  nature  of  the  sub-grid  scale  stress  effects, 
but  empirically  by  requirements  that  nonlinear  computational  instability 
be  eliminated.  Unfortunately,  in  many  cases  it  has  been  difficult  tc 
decide  which  cf  the  f 3-uctuacicr.s  observed  in  the  numerical  solutions 
were  a  result  of  numerical  instabilities,  and  which  corresponded  to 
natural  instabilities  cf  the  atmosphere.  Since  the  correct  physical 
representation  cf  the  sub-grid  scale  effects  is  not  known  with  cer¬ 
tainty,  but  is  of  the  same  form  as  that  required  for  numerical  sta¬ 
bility,  the  situation  can  become  obscured  very  rapidly.  In  particular, 
the  longer  the  desired  prediction  time  of  the  simulation  the  more  im¬ 
portant  are  the  sub-grid  scale  effects  physically,  and  the  greater  is 
the  need  for  mathematical  damping  of  possible  numerical  instabilities. 

The  resolution  of  this  difficult  problem  requires  not  just  a  better 
method  cf  dealing  with  numerical  instability,  but  a  more  adequate  de¬ 
scription  of  the  sub-grid  scale  dynamics.  The  importance  of  the  latter 
for  long-term  'weather  prediction  is  currently  a  matter  of  highest  in¬ 
terest  and  priority  since  it  involves  the  question  of  the  ultimate 
forecasting  capability  of  the  presently  envisaged  numerical  forecasting 


scherr.es.  The  central  question  here  is  the  rate  at  which  energy  can 
be  transferred  frcrr.  the  small  (sub-grid)  to  the  large  (synoptic) 
scale  notions.  If  this  figure  is  as  large  as  some  estimates  now 
place  it,  the  current  forecasting  grid  mesh  gives  a  time  limit  of 
about  three  days  for  weather  forecasting.  A  resolution  of  this 
particular^  question  is  hoped  for  in  the  near  future  from:  measure¬ 
ments  of  the  energy  spectrum,  of  the  atmospheric  motions.  If  long- 
range  forecasting  is  to  become  a  reality  confirmatory  evidence  for 
a  "gap"  in  the  energy  spectrum:  between  the  m.eso  and  synoptic  scales 
that  is  believed  to  exist  is  needed,  net  just  within  the  atmospheric 
boundary  layer,  where  most  observations  are  made,  but  throughout 
the  troposphere. 

This  characteristic  of  the  energy  spectrum,  is  also  important 
'when  considering  requirements  fer  initial  data  inputs  to  the  numer¬ 
ical  model.  Thus,  'with  conventional  meteorological  sensing  instru¬ 
ments,  measurements  are  made  at  fairly  regular  time  intervals,  of  the 
order  of  an  hour,  say.  but  with  an  output  that  corresponds  to  ar. 
average  value  of  the  meteorological  quantity  oaken  over  both  a  small 
physical  volume  and  a  short-time  interval  (seconds  or  minutes).  When 
such  data  are  interpreted  in  terms  of  much  longer  dimensions  and 
longer  times,  corresponding  to  the  space-time  grid  net  for  the  com¬ 
putational  simulation,  aliasing  occurs.  Thus  the  mathematical  model 
precludes  the  description  of  sub-grid  ...ale  wavelength  ana  frequencies, 
tc  v/h i ch  trh.s  rrscisu^'d-nc  (i, n s t gu p c n tr s  dr*c  clsdi?ly  rs  s poncing .  This 
specific  effect  of  the  use  cf  initial  data  of  this  variety  cn  the 
computer  simulation  depends  on  the  energy  spectrum,  cf  the  motions 
between  the  hourly  and  minute  scales.  If  the  high  frequency  content 
is  low,  aliasing  is  also  lew,  but  in  fact,  as  shown  in  Fig.  37  there 
appears  to  be  a  peak  in  the  energy  spectrum  around  1  minute.  In  an 
analysis  by  Cort  and  Taylor  (1969),  it  is  demonstrated  that  with  a 
lower  period  cut-off  cf  2  hours  in  the  frequency  response  of  the  data 
input  to  the  computer,  an  error  in  the  2-hour  period  energy  spectrum 
content  of  the  order  of  a  factor  cf  10  is  possible.  In  terms  of  the 
computer  simulation  this  type  of  result  implies  very  large  errors  in 
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the  highest  frequency  (grid  scale)  description  of  the  simulation,  and 
a  subsequent  degradation  of  the  prediction  accuracy  as  the  prediction 
time  increases.  If  the  atmospheric  energy  spectra;:,  is  indeed  as  de¬ 
scribed  Vy  Van  dor  Ho  von  (I'J'i?)  and  Oort  and  Taylor  (1969),  the  only 
method  of  dealing  with  this  difficulty  (apart  fro;:  reducing  the  time 
stop  in  the  simulation)  is  to  change  the  T.othoo  cf  taking  meteorolog¬ 
ical  data,  so  that  a  measure  sensitive  to  the  highest  frequencies  but 
averaged  over  space  and  time  steps  corresponding  approximately  to  the 
grid  mesh  sine  of  computation  is  obtained. 


FIGURE  37,  Results  of  an  Attempt  to  Correct  the  Spectrum  at  Caribou,  Maine,  for  the 
Effects  of  Aliasing  from  Periods  Between  1  Min  (the  Basic  Averaging 
Period  of  the  Wind  Reports)  ond  2  Hr  (the  Nyquist  Frequency).  Frequency 
F  in  Cycles/4096  Days  (After  Oort  ond  Taylor,  1969) 


An  alternative  method  aimed  at  eliminating  nonlinear  ccmpututi or 
instability  Is  ti.o  use  of  energy  conserving  finite  difference  schemes 
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(see  Table  8)  was  designed  specifically  to  accomplish  this.  By  the 
placing  of  certain  integral  restraints,  not  just  total  energy,  but 
vorticity  and  zonal  momentum,  are  conserved  exactly  in  this  model  for 
a  frictionless  nen-diss ipative  atmosphere.  The  results  of  numerical 
computations  with  this  system  shew  no  evidence  of  instabilities,  even 
in  the  absence  of  dissipative  terms. 

An  analogous  approach  is  taker,  in  the  Kasahara-Washingtor.  model 
(see  Table  8)  where  a  Lax-Wendroff  energy- conserving  finite  difference 
scheme  is  used.  The  nonlinear  stability  of  this  method,  however,  has 
not  teen  demonstrated. 

Whichever  of  the  finite  difference  schemes  is  used,  there  remains 
a  doubt  as  to  its  mathematical  convergence  and  stability  properties, 
since  they  are  all  too  complex  to  allow  rigorous  mathematical  proof, 
and  basically  the  test  of  any  one  cf  them,  has  been  made. by  a  trial  run 
and  comparison  cf  the  results  with  the  actual  atmosphere.  Unfortunately, 
natural  irregularities  of  the  atmosphere  are  so  prevalent  that  such  a 
comparison  is  never  perfect  and  there  is  always  the  possibility  that 
some  of  the  variations  predicted  by  the  model  are  of  a  mathematical 
and  not  a  physical  variety.  In  the  final  analysis,  and  in  the  absence 
of  mathematical  convergence  proof,  the  only  test  of  the  numerical  solu¬ 
tion  that  approaches  conclusiveness  is  a  systematic  series  of  trail 
calculations  in  which  all  initial  and  boundary  conditions  are  kept  the 
same  while  the  grid  mesh  dimensions  are  changed.  The-  large  running 
times  of  the  presently  formulated  atmospheric  circulation  models  have 
prevented  a  test  of  this  type  fi  v>eing  made  for  any  of  them.. 

Boundary  Conditions 

The  system  of  Eqs.  C.2  -  8.6  or  its  finite  difference  represen¬ 
tation  has  to  be  supplied  with  suitable  boundary  conditions.  The 
simulations  that  are  of  concern  here  are  global,  sc  that  no  lateral 
boundary  conditions  need  be  specified.  The  upper  altitude  limit  for 
tropospheric  models  is  at  the  tropopauso ;  it  is  assumed  that  there  is 
no  influence  of  the  upper  atmosphere  acting  at  the  boundary  on  the 
troposphere.  The  .  ewer  boundary  condition  is  normally  applied  above 

1  '..-7 


the  Ekman  layer,  rather  than  at  the  earth’s  surface.  Hence  the 
boundary  conditions  have  to  allow  for  the  turbulent  transports  of 
heat  and  momentum  that  occur  in  the  EVr-an  and  Prandtl  layers  between 
the  earth  and  the  leanest  horizontal  finite  difference  mesh  surface. 
T'no  semi-empirical  expressions  that  cecide  these  effects  are  con¬ 
tained  in  the  following  lower  boundary  conditions: 


So 

Bo 

5v 

hz 


CD 


"  CD 
=  C~ 


us  +  vs 


Vs  vu|  +  v| 


(Ts  -  Tz)  Vu|  + 


(3.1? 


Here  F^,  is  the  vertical  heat  flux,  Tc  the  ground  temperature,  Cp  the 
surface  friction  coefficient,  ana  the  subscript  s  refers  to  values  at 
the  "surface'1,  i.e,,  at  the  outer  edge  of  the  lower,  Prandtl  layer. 

The  boundary  layer  specification  is  completed  by  adding  the  condition 
w  =  0  to  3.1?. 

The  surface  velocities  u,,  and  vs  are  usually  obtained  by  linear 
extrapolation  from  the  solution  at  higher  levels,  though  Smagorinsky 
uses  a  more  sophisticated  extrapolation  method  that  allows  for  rotation 
of  the  flow  through  the  Ekman  layer  (cf.  Table  S>.  Additional  terms 
for  the  vertical  heat  flux  can  arise  if  humidity  and  water  vaporization 
effects  are  allowed  for.  Indeed,  if  they  are,  their  contributions  also 
must  be  included  in  the  heat  source  term  q  in  Eq.  8.5.  Volume  IV  dis¬ 
cusses  this  briefly,  and  an  indication  of  how  the  humidity  effects  are 
allowed  for  ir.  the  various  models  is  given  by  Table  8. 

Characteristics  of  Current  Tropospheric  Numerical  Simulations 

Table  8  is  a  listing  of  some  of  the  characteristics  of  the  four 
primitive  global  numerical  circulation  simulations  for  the  troposphere 
that  have  beer:  constructed  to  date,  These  models  are  used,  not  for 
operational  forecasting,  but  for  development  of  forecasting  techniques 
for  the  study  of  climatic  problems .  Other  less  comprehensive  models 
are  used  by  the  Weather  Bureau  for  daily  forecasting  in  the  United 
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States;  the  Air  Force  and  the  Navy  have  comparable  models  for  Service 
use.  These  operational  weather  predict! cr.  models  are  not  global  in 
coverage,  and  hence  of  less  interest  for  this  report  t'nsr.  the  larger 
global  simulations,  since  it  is  probable  that  upper  atmospheric  simu¬ 
lations  require  very  large  if  not  global  coverage. 

All  the  Simulations  use  the  primitive  equations,  6.2  -  8,7,  with 
empirical  or  serrd-eo.pirical  values  for  the  stress  terms  Fx  and  Fy ,  and 
i,  and  heat  source  inputs  to  q  from,  radiative  heating  and  from  the 
moisture  field.  The  latter  contributions  can  be  treated  in  great  de¬ 
tail  and  complexity;  indeed  there  have  been  several  large-scale  numer¬ 
ical  simulations  of  radiative  transfer  ir.  the  atmosphere ,  in  which  the 
fluid  motions  are  ignored.  A  brief  description  of  the  type  of  radiative 
sources  and  sinks  taken  for  the  tropospheric  circulation  models  is  in¬ 
cluded  in  Table  8;  a  further  general  discussion  of  this  topic  is  con¬ 
tained  in  Volume  IV. 

Also  shewn  in  the  table  are  the  principal  characteristics  of  the 
finite  difference  schemes  used,  their  computational  properties  and 
machine  running  times.  A  sample  of  the  sequence  cf  operations  run 
through  for  solving  the  finite  difference  equations  is  illustrated 
with  the  Leith  model. 

Extension  to  Upper  Atmospheric  Circulation 

The  survey  cf  the  tropospheric  circulation  simulation  problem, 
given  above  illustrates  the  considerable  magnitude  of  that  task. 

From  the  figures  quoted  in  Table  S  for  machine  requirements  one  might 
guess  that  an  extension  cf  the  tr  ' aspheric  model  to  include  strato¬ 
spheric  circulation  would  place  the  computational  load  outside  the 
capabilities  of  present  day  computers.  Besides  the  increase  in  volume 
coverage  implied  by  increasing  the  altitude  span  of  the  simulation , 
the  greater  ruin, her  of  molecular  constituents  above  the  tropopause  ana 
the  need  for  inclusion  of  their  variation,  through  chemical  reactions 
as  well  as  their  radiative  properties,  lead  to  an  enormous  increase 
in  computer  requirements  for  simulation.  Later  in  this  section  esti¬ 
mates  will  be  made  to  indicate  the  scope  cf  the  problem;  these  indeed 
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not  just  the  effects  of  chemical  changes ,  but  the  electrcdynamie  body 
forces  that  are  needed  re  account  correctly  for  the  dynamics  of  the 
ionised  constituents  of  the  upper'  atmosphere .  With  a  multicomponent 
fluid,  momentum  and  energy  equations  are  required  for  each  species, 
and  there  are  in  addition  'rate  equations  linking  them  and  ohmic  less 
equations  for  the  ionized  components .  It  should  particularly  be  noted 
that  the  hydrostatic  assumption  is  no  longer  applicable  for  the  upper 
atmosphere,  due  to  the  presence  of  large  vertical  winds,  sc  that  in 
place  of  the  Richardson  Eq.  8.8,  there  now  appears  a  fourth  prognostic- 
equation  for  the  vertical  acceleration.  In  general,  the  high  wine 
velocities  are  the  source  of  gravity  wave  motions  that  constitute  an 
important  element  in  upper  atmosphere  dynamics. 
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In  addition  to  those  additional  features  of  the  upper  atmosphere 
which  obviously  can  be  a  source  of  difficulty  in  the  r.unorical  mathe¬ 
matics  of  the  simulation,  the  state  of  knowledge  of  the  physics  of 
the  upper  atmosphere  has  t:  he  considered  in  assessing  the  feasibility 
of  a  model.  Up  to  about  100  km  it  appears  that  this  is  sufficient, 


On  the  positive  side,  however,  competence  to  address  the  problem 
is  improving  steadily  as  evidenced  by  progress  of  the  next  Generation 
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for  a  simulation  of  the  upper  atmosphere  to  an  altitude  of,  say, 

2'jO  km  is  considered.  It  is  first  necessary  specify  in  time  and  ir. 
distance  the  computational  grid  spacing  for  which  sere  criteria  are 
given  in  Table  9,  These  criteria  lead  to  the  conclusion  that  the 
computation  should  be  carried  out  with  a  spacing  which  is  vertically 
less  than  about  1 j  km,  horizontally  less  than  1C 00  km,  and  at  time 
steps  which  are  less  than  IOC  sec.  A  better  simulation  than  the 
minimum  described  above  would  employ  a  grid  specification  of  ICOu  km , 

2  vertical  spacing  of  7  km,  and  a  time  step  of  10  sec. 

Shown  m  Table  1  is  an  estimate  of  the  number  of  variable  points 
required  for  upper  atmosphere  simulation  in  the  altitude  ranges  up  tc 
200  km.  The  number  of  points  which  would  be  required  in  the  vertical 
grid  varies  from  one  altitude  region  to  another  according  to  Nyquist 
criterion  applied  tc  the  temperature  stratification.  By  the  Nyquist 
criterion,  at  least  two  data  points  are  required  rcr  each  cycle  cf  the 
variation  encountered.  In  each  case,  about  SOu  points  are  required  in 
the  horizontal  grid  for  a  100  km  or  9-dec;  spacing. 

The  number  of  dependent  variables  per 
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atmospheric  regions  is  of  course  different 


cependinc  up: 
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leant  variables,  i.e.  ,  components  of  velocity  vectors,  pressure,  density 
'~1  principal  contituents ,  and  heat  absorption  by  principal  atmospheric 
a.iscrbers .  The  total  number  of  dependent  variables  per  horizontal  grid 


TABLE  9.  CRITERIA  FOR  GRID  SPECIFICATION 


TABLE  10.  VARIABLE  POINTS  FOR  UPPER  ATMOSPHERE  SIMULATION 


ALTITUDE  RANGE  (KM)  | 

0-20 

20-50 

50-90 

90-130 

150-200 

INDEPENDENT  VARIABLES 

4 

4 

4 

4 

4 

S',8,X,tl 

Sr,9,X  ,1  ) 

(r,8,X  ,n 

<r,e,X,t) 

(r,6,X,f  1 

NUMBER,  points  in  noriz.  grid 

800 

800 

800 

800 

800 

3 

(  10  km  or  9  6tg  ipocirvg  ) 

NUMBER,  POINTS  IN  VERT.  GRID 
{  7  km  ipocing  ) 

3 

4 

6 

9 

7 

DEPENDENT  VARIABLES  PER  POINT 

’0 

40 

40 

20 

20 

(U,V,W;N,p 

(U.V.W^p 

(  U ,  V ,  W  ,  N ,  p 

( U,V,W,N.p 

(U,V,W,N,p 

nco2'  nH2o- 

nN'nN2' 

nN, '  "Nj' 

nN2‘  nN' 

nN2‘  nN' 

<,h20-  qC02' 

n03P'nOID' 

n03P'nOID' 

n03P'  "OID' 

n03P'NOID' 

qSURF) 

n°2,n°3' 

nO.'nO  ' 
i  3 

"°2 

H  _ 

°2 

X 

c 

X 

O 

c 

nOH'nH' 

■*  3  q'l 

♦  3  q'* 

nH20'°  H.02' 

nH20'nH.02’ 

+  5  K*« ) 

♦  5  K'i  ) 

o 

z 

c 

"  <N 

o 

X 

c 

nH02'nN0' 

nN02 

nN02 

+  6  q'l 

4  6  q*» 

+  13  K‘i) 

*  13  K’t) 

DEPENDENT  VARIABLES  6{10)5 

2A(  10  )3 

128  <  10  )3 

192(IO)3 

1A4  ( 10  )3 

1 12  (10)3 

1L' 


=3 


point  in  the  altitude  range  from  0-200  km  is  approximately  6x1 .  Or. 
the  basis  of  the  data  given  in  Table  10,  there  are  shewn  in  Table  11 
estimates  cf  computer  requirements  for  both  the  complete  simulation 
anu  for  a  minimum  simulation  of  a  more  restricted  nature.  Ey  com¬ 
parison  are  analogous  estimates  for  current  simulations  (NAS -HRS 
Fufclic.au ior.  1350,  1966)  cf  turbulence,  convection,  tropospheric 
circulation,  and  ocean  circulation.  On  the  assumption  that  no  one 
can  justify  830  days  of  computer  time  to  simulate  ICO  cays  of  upper 
atmosphere  circulation,  it  is  plain  from  Table  11  that  current  com¬ 
puters,  represented  by  the  CDC-6600  which  has  an  operation  time  of 
600  nsec  per  operation,  are  clearly  inadequate  for  the  purposes  of 
simulating  the  global  upper  atmosphere.  Such  computers  as  indicated 
by  the  circle  require  approximately  390  hr  of  computer  simulation 
time  for  each  100  days  of  tropospheric  weather  simulated. 

TABLE  11.  ESTIMATES  OF  COMPUTER  REQUIREMENTS 


UPPER  ATMOS.  SIMUL. 
COMPLETE  |  MINIMUM 
SIMULATION  SIMULATION 

TURBULENCE 

CURRENT  S 
(NAS-NRC  p 
CONVECTION 

mulations 
JBL.  1350  1966) 
TROPO.  CIRC. 

OCEAN  CIRC. 

FILTERING  APPROX 

FOURIER 

FOURIER 

HYDROSTATIC 

GEOSTROPHIC 

DIMENSIONALITY 

4 

4 

4 

4 

4 

4 

DEP.  VARIABLES 

10  -  40 

15 

4^  jt  v,  w,  r ) 

5(u,v,w,T,r> 

4(  U/  v,  T,  r ) 

4(u, v/ Tf » ) 

NUMBER,  VERllCAL  LEVELS 

30  (7  KM) 

20  (10  KM' 

30 

100 

10 

10 

NUMBER,  HORIZ.  POINTS 

600(I03KM) 

800  (103KM> 

104 

io4 

io4 

io4 

TIME  STEP  (wcl 

10 

io2 

0.001 

5 

300  (  5  MIN  ) 

1 .  B(  10)4(  5  HO 

SIMULATED  EXP.  TIME  (»c] 

I0?(  lOOd) 

107(100d) 

10 

3,6(10  )4 

107(  lOOd  } 

10%  30  YO 

DEP.  VER.  PER  TIME  STEP 

6(10  )5 

2,4(10  )3 

!04 

(5  H.) 

500  )6 

4(10  )5 

4(10  )5 

TIME  STEP  PER  EXP. 

io4 

!05 

io4 

200  )3 

3(10)4 

6(  10 )“ 

DEP.  VAR  PER  EXP. 

61  10)" 

2.4(10)’° 

io10 

io'° 

1.200)'° 

2.4(10)’° 

COMP.  OPER.  PER  EXP. 

1.2(  IO)'4 

4. BOO)'2 

2(10)'2 

2(  10 )' 2 

2.4(10)12 

4.8O0)'2 

(200  CPER/DEP.  VAR) 

COMP.  TIME- CDC6600 

7. 2(10  )7 

2.9O0)6 

1 . 2(  1 0  )4 

1.2(10)4 

i.4(10r\ 

2.9( 10 16 

(600  mec  PER  OPER) 

(2(  IO)4  H) 

( 8 10H ) 

( 330H ) 

t  330m  ) 

U  390H)  J 

(810M) 

COMP.  TIME-ILIAC  (Kt) 

( 830d 1 

f 1.4  (10)5'b 

5.800  >3 

2.4( 10  )3 

2.4(10  }3 

2,9( 10  )3 

5.8i  10)3 

(1.2  ntftc  PER  OPER ) 

QQ39  H ) 

0.6  H) 

(  40  Mi  N  ) 

(40  MIN) 

(48  MIN  J 

( 1.6  Hr) 

]  J  3 


Tiio  i  ••  n-  i  i  >i.  cf  computers 
expected  to  he  opera ti  oral  in  the  ear 
may  require  only  0?.  m.in  for  ••  simulat 
circulation  of  103  days  duration.  It 


however,  such  as  the  ILIAC  IV, 
y  l'j/Q’s  and  500  times  as  fast, 
or.  of  global  tropospheric 
would  alsc  be  capable  cf  simu¬ 


lating  100  cays  of  upper  atmosphere  circulation  in  33  hr  cf  computer 


-day  allocations  to  ether  problems. 


Smaller  scale  computations  than  considered  in  this  estimate  are 
also  possible.  For  example,  calculations  over  a  restricted  volume, 
either  vertically  or  hori mentally ,  may  be  of  interest  in  some  appli¬ 
cations.  Ir.  such  a  case  the  scope  of  the  computation  is  reduced, 
though  at  the  moment,  knowledge  cf  the  energy  transfer  rates  is  in¬ 
sufficient  tc  enable  prediction  time  to  be  related  to  the  volume  of 
the  atmosphere  that  is  simulated. 


Alternative  numerical  methods  cf  solving  the  equations  perhaps 
should  also  be  considered  and  a  particularly  attractive  approach 
from  the  point  of  view  of  the  aliasing  and  sub-grid  scale  energy 
spectrum  effects  is  the  use  cf  wave  number  rather  than  physical 
space.  Iii  such  a  method  a  large  wave  number  cut-off,  corresponding 
to  the  maximum  resolvability  of  the  rinite  difference  grid  mesh,  can 
be  easily  applied,  thereby  avoiding  some  of  the  confusion  that  results 
with  the  methods  used  tc  date.  Prior  to  the  advent  of  the  fast 
Fourier  transform  technique  such  a  scheme  would  have  implied  excessive 
computational  time ,  but  this  now  appears  no  longer  to  be  the  case. 


mention  should  also  be  made  cf  a  different  variety  of  computer 
simulation  that  could  be  considered.  This  is  the  problem  of  analysis , 
i.e.,  the  estimation  of  current  atmospheric  conditions  from  a  limited 
set  cf  observations,  in  contrast  to  prediction ,  with  which  the  prior 
discussions  in  this  section  have  been  concerned.  Although  common  in 
daily  meteorological  work,  it  appears  that  little  thought  has  been 
given  to  dealing  with  this  question  using  numerical  methods  and  the 
fluid  dynamic  equations.  This  question  appears  to  be  cf  potential  im¬ 
portance  for  'the  upper  atmosphere ,  in  which  for  some  applications  we 
may  think  of  current  satellite  information  on  atmospheric  data  being 
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analyzed  by  computer  to  give  an  instantaneous  picture  of  the  cot  pie  te 
state  of  the  atmosphere .  Limitation  to  analysis,  rather  than  pre¬ 
diction,  eases  the  problem,  frerr.  the  point  of  view  of  computer  capacity 
requirements  and  eliminates  sene  of  the  difficult  questions  of  the 
propagation  of  sub-grid  scale  effects  into  the  synoptic  scales. 


As  was  just  stated,  little  attention  has  beer,  given  to  this  ap¬ 
proach,  and  hero  we  shall  Just  point  cut  a  cross  difference  in  the 
numbericai  calculation  task  from,  that  for  the  prognostic  problem. 

Thus  in  "Typical  Wave  Motions"  of  Volume  IV,  it  -was  shown  that  the 
equations  governing  the  spatial  domain  could  be  of  a  hyperbolic, 
parabolic  or  elliptic  variety,  depending  on  the  region  of  wave  motion 
involved  (inertial,  gravitational  and  acoustic).  Thus,  assuming  these 
considerations  are  transferable  to  the  non-linear  equations,  the  solu¬ 
tion  techniques  have  to  be  adjusted  according  to  the  domain  cf  solu¬ 
tion,  and  this  is  quite  different  from  the  forecasting  problem,  in 
which  a  straightforward  marching  technique  is  always  used.  Further 
elaboration  and  study  cf  the  use  of  this  type  of  analysis  need  further 
attention,  since  it  appears  to  be  cf  potential  significance. 


To  Summarize; 


The  implicit  determination  of  winds  by  computational  simulation 
appears  to  have  potential.  The  justifiable  degree  of  complication 
in  such  simulation  is  questionable  requiring  further  investigation. 
Certainly  th ;  extensive  use  of  computers  for  the  analysis  of  weather 
data  is  amply  .notified.  The  analysis  cf  the  entire  global  atmosphere 
up  tc  an  altitude  of  200  km  seems  feasible  from  the  standpoint  of  the 
capability  of  next  generation  computers  which,  would  require  nc  more 
computer  time  than  the  exercise  of  current  problems.  Accomplishment 
may  be  slowed  by  the  fact  that  present  knowledge  of  the  physics  cf 
the  upper  atmosphere  is  incomplete . 
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1,  Introduction 

The  capability  of  achieving  a  refined  knowledge  of  the  variations 
of  atmospheric  conditions  on  a  global  scale  and  near  real-time  is 
within  the  grasp  of  modern  civilization  through  the  use  of  meteoro¬ 
logical  satellites  and  high  speed  computers.  Such  refined  knowledge 
in  the  above  space  and  time  scales  would  be  derived  from  a  process  that 
involves  two  phases:  (1)  intermittent  global  measurements  in  near 
real-time  cf  atmospheric  conditions  from  mainly  tt  serological  satel¬ 
lites,"'  and  (2)  subsequent  predictions  of  the  variations  cf  the  atmos¬ 
pheric  parameters  in  the  lower  atmosphere  (i.e,,  the  troposphere  or 
altitudes  up  to  about  8  nmi)  during  a  time  that  appears  to  be  limited 
to  about  three  weeks  (Refs .  1,2'?.  The  overall  economical  reward  from  a 
refined  knowledge  of  atmospheric  conditions  that  would  allow  even  a  one 
week  weather  prediction  has  been  estimated  in  billions  of  dollars  as  a 
consequence  of  improved  management  of  water  resources ,  agriculture  , 
surface  transportation,  retail  marketing,  the  lumber  industry,  etc. 

(Ref.  i)  . 

The  need  for  global  measurements  becomes  clear  from  the  rather 
straightforward  formulation  of  the  air  motion  within  the  troposphere, 
i.e.,  the  motion  of  air  cn  a  rotating  earth  receiving  shortwave  radia- 
tion  from  the  sun  and  returning  longwave  radiation  tc  space.  The 
dependent  variables  cf  the  problem  are  the  three  components  cf  the  air 
velocity  (i.e.,  wind),  the  pressure,  temperature ,  air  density,  and 
water  vapor.  The  independent  variables  are  given  by  the  throe  space 

"Meteorological  satellites  can  reach  otherwise  inaccessible  areas  cf 
the  oceans. 


coordinates:  longitude,  latitude,  and  altitude  and  the  time.*  The 
seven  dependent  variables  are  determinable  as  a  function  of  die  four 
independent  variables  from  the  use  of  seven  equations  obtained  as 
follows :  three  momentum  equations ,  the  conservation  of  energy  within 
a  nonadiabatic  troposphere,  the  conservation  of  mass,  the  equation  of 
state  for  a  real  gas  with  a  molecular  weight  of  air,  and  the  diffusion 
equation  for  the  concentration  of  water  vapor  within  the  troposphere 
coupled  with  the  partial  pressure  of  water  vapor  as  a  function  of  pres¬ 
sure  and  temperature.  This  set  of  partial  differential  equations  re¬ 
quires  initial  conditions**  for  each  of  the  seven  variables  and 
boundary  conditions**  that  take  into  account  the  energy  exchange  at  the 
surface  of  the  earth  (between  the  atmosphere  and  botn  the  ground  and 
the  oceans),  the  orographic  effects  on  the  air  velocity  vector,  and  the 
energy  exchange  at  the  top  of  the  troposphere.  While  this  set  of 
partial  differential  equations  is  difficult  to  solve  even  by  numerical 
methods,  the  feasibility  of  obtaining  numerical  solutions  from  them  has 
already  been  demonstrated  under  somewhat  restricted  conditions  (Refs.  2,4 
to  7).  This  feasibility  of  solution  is  further  demonstrated  from  the 
proven  success  in  a  somewhat  similar  numerical  situation,  i.e.,  the 
numerical  solution  of  a  quasi-steady  flow  within  laminar  hypersonic 
boundary  layers  with  mass  injection  from  the  receding  surface  of  an 
ablating  intercontinental  ballistic  missile  (Refs.  8,  9). 

The  need  for  near  real-time  processing  of  the  global  measurements 
stems  from  the  relative  magnitude  of  the  force  terms  in  the  momentum 
equations  for  the  two  horizontal  components  of  the  air  velocity  vector, 
i.e.,  the  time  dependent  term  (which  is  essential  for  time  predictions) 
or  air  acceleration  is  given  by  the  difference  of  two  nearly  equal 
forces,  namely  the  pressure  gradient  and  Coriolis  forces.  This  fact, 
together  with  inherent  errors  of  both  experimental  and  analytical 

*The  pressure  and  altitude  are  usually  interchanged  as  the  independent 
and  dependent  variable,  respectively.  This  result  stems  from  the 
direct  coupling  of  these  two  variables  through  the  hydrostatic  equa¬ 
tion. 

**As  deduced  from  satellite  and  other  available  data. 

124 


nature,*  puts  a  limit  to  the  length  of  time  for  accurate  predictions, 
i.e.,  the  predicted  behavior  of  atmospheric  parameters  as  derived  from 
initial  global  measurements.  The  subsequent  use  of  the  equations  of 
motion  is  estimated  to  be  accurate  during  a  time  that  is  no  longer  than 
about  three  weeks  after  the  measurements  (Refs.  1,  2).  Hence  ,  the 
perishable  data  for  the  initial  and  boundary  conditions  must  be 
utilized  in  the  numerical  integrations  of  the  equations  of  motion  in 
near  real-time  in  order  for  its  utility  to  be  maximized  for  the  pre¬ 
diction  of  the  behavior  of  the  atmospheric  parameters  during  a  rather 
short  subsequent  time. 

While  the  development  of  numerical  techniques  to  describe  the  air 
motion  in  different  scales  within  the  troposphere  is  a  subject  of  great 
interest,  the  subsequent  treatment  is  limited  to  the  development  of 
meteorological  satellites  and  sensor  subsystems  for  the  remote  measure¬ 
ment  of  the  initial  and  boundary  conditions  in  the  global  and  near 
real-time  scales;  that  is  to  say,  meteorological  satellites  operating 
in  medium  altitude  orbits  (e.g.,  in  the  range  from  400  to  800  ami)  and 
in  geosynchronous  orbits  at  19,300  nmi. 

The  subjects  of  interest  here  include  the  following:  (a)  feasi¬ 
bility  of  remote  measurements  of  the  initial  and  boundary  conditions 
within  the  required  degree  of  accuracy  in  a  global  and  real-time  scale; 
(b)  chronological  developments  of  U.S.  meteorological  satellites  and 
sensor  subsystems;  (c)  current  projections  for  the  development  of 
sensor  subsystems  for  U.S.  meteorological  satellites  during  the  forth¬ 
coming  decade;  (d)  brief  description  of  the  international  GARP  (Global 
Atmospheric  Research  Program)  and  WWW  (World  Weather  Watch)  programs, 
and  (e)  conceptual  techniques  for  potential  measurements  from  satel¬ 
lites  . 


:;Errors  in  the  measurement  of  initial  conditions  as  well  as  truncation 
of  the  difference  equations  and  rounding  off  during  the  numerical 
integration . 


2  .  Feasibility  of  Remote  Measurements  With  Appropriate  Accuracies 

The  remote  measurement  cf  the  initial  conditions  in  a  global  scale 
and  real-time  involve  the  following  variables:  geopotential  altitude 
(h)  of  a  given  pressure  level  (p,  ir  millibars),*  temperature  (T), 
density  (p ) ,  water  vapor  Ox),  and  the  three  components  of  the  velocity 
vector  (u ,  v ,  w  in  the  longitude ,  latitude ,  and  vertical  directions , 
respectively).  Since  the  density  is  related  to  the  pressure  and 
temperature  through  the  equation  of  state  for  a  real  gas,**  the  density 
becomes  determined  from  the  pressure  and  temperature  distributions. 
Likewise,  the  determination  of  the  three  velocity  components  requires 
determination  of  the  magnitude  of  the  velocity  vector  (V),  i.e.,  as 
fixed  by  both  the  magnitude  and  the  direction  of  the  wind  vector.  The 
accuracy  for  the  measurement  of  the  initial  conditions  has  been  speci¬ 
fied  by  the  1967  GARP  study  conference  as  follows  (Ref.  11) :+ 

a.  h(x,  y,  t),  as  fixed  by  a  reference  pressure  with  an  r.m.s. 
error  of  2  percent. 

b.  T(x,  y,  p,  t),  with  an  r.m.s.  error  of  ±  1°C. 

c.  x(x,  y>  P,  t),  with  an  r.m.s,  error  of  10  percent. 

d.  V(x,  y,  d ,  t),  with  an  r.m.s.  error  of  ±  3  m/sec. 

where  x  denotes  the  latitude,  y  the  longitude,  and  t  the  time.  The 
horizontal  resolution  of  these  global  measurements  was  specified  at  a 
typical  value  of  about  400  km  (~200  nmi),  while  the  vertical  resolution 
was  set  at  200  mb  in  an  altitude  range  corresponding  up  to  10  mb."r'r 

The  remote  measurements  of  the  boundary  conditions  involve 
measurements  of  the  energy  exchange  at  the  interface  between  the 


*As  consequence  of  the  interchange  of  pressure  and  altitude  as  the 
independent  and  dependent  variables,  respectively. 

**The  gas  constant  can  be  made  to  include  the  effect  of  water  vapor  in 
the  troposphere,  which  is  at  most  of  the  order  of  4  percent  (Ref.  10) 

+See  Ref.  11  to  identify  the  outstanding  membership  of  this  study. 

^+Since  one  bar  or  1000  mb  is  approximately  the  standard  atmospheric 
pressure  at  the  surface  of  the  earth,  the  altitude  range  extends  from 
the  surface  of  the  earth  to  about  100,000  ft  or  16.5  nmi,  approxi¬ 
mately. 
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atmosphere  and  oceans  as  well  as  the  global  radiation  balance  ( i.e.  , 
including  cloud  cove). 

It  is  shown  be  lev.  that  the  geopotential  altitu.k  (h)  and  the  dis¬ 
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titude  orbits,  and  (2;  on  uprated  statistical  .ample  of  data  for  the 
geopocential  altitude  arc  temperature  applicable  to  the  local  region 
of  interest.  Similarly,  the  oistribution  of  ■  .  .v  vapor  ( uu)  through 

the  lower  troposphere  is  also  determinable  fro.  ..emote  radiance  meas¬ 
urements  in  the  infrared.  Potential  techniques  for  the  determination 
of  the  wind  vector  are  as  follows:  (a)  d-. testing  the  motion  of  cloud 
elements  that  move  with  the  wind  in  the  lower  troposphere ,  as  observed 
from  geosynchronous  meteorological  satellites;  and  (b)  using  a  world¬ 
wide  network  of  constant -density  balloons  distributed  at  different 
pressure  (i.e.,  altitude)  levels  in  the  upper  cloudless  troposphere, 
and  tracking  the  balloons  from  meteorological  satellites." 

3.  R.K.S.  Errors  ir.  the  Geopotential  Height 

It  has  beer,  shown  that  the  geo  potential  height  or  altitude  profile 
is  uniquely  related  to  the  temperature  distribution  in  the  atmosphere 
by  the  following  mathematical  relationships  (Ref.  12). 

a.  The  hydrostatic  equation,  which  relates  the  local  drop  ir. 
pressure  with  increasing  altitude  to  the  local  density; 

b.  The  equation  of  state  for  a  real  gas,  which  relates  the  local 
pressure  to  the  local  density  and  temperature; 

c.  Hie  empirical  orthogonal  functions  (Refs.  13,  14),  which 
indicate  that  they  are  statistically  stable  in  that  they  are 
nearly  invariant  with  respect  to  the  latitude,  longitude, 
and  time. 


"Potential  techniques  to  measure  the  wind  vector  ever  the  oceans  in¬ 
volve  determination  of  the  sea  surface  conditions  from  either  micro- 
wave  measurements  or  the  solar  alitter. 
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Thus ,  subject  to  the  validity  of  the  empirical  orthogonal  func¬ 
tions  as  derived  from  the  statistical  behavior  of  the  atmosphere,  the 
geopotential  height  can  be  specified  solely  from  the  temperature  pro¬ 
file;  which  as  treated  subsequently  is  determined  from  radiance  measure¬ 
ments  from  medium  altitude  orbits . 

The  above  method  has  been  tested  by  using  extensive  radiosonde 
data  for  four  consecutive  days.  The  first  three  days  provided  299  sets 
of  observations ,  which  were  used  as  the  required  statistical  or 
dependent  sample  of  observations.  The  fourth  day  provided  data  in  the 
form  of  actual  temperature  profiles  and  calculated  radiance  to  be  used 
as  postulated  measurements ;  the  latter  ones  were  then  used  in  the  pro¬ 
cedure  outlined  above  to  determine  calculated  temperature  profiles. 
Figure  A-l  compares  the  r.m.s.  differences  between  the  radiosonde  and 
calculated  values  for  the  geopotential  height  under  clear  and  cloudy 
conditions  for  two  assumed  values  of  the  random  errors  in  the  measure- 
ments,  i.e.  ,  ce  =  10  and  cg  =  1/4  ergs/cm  -sec-strdn-cm  . "  The 
standard  deviation  of  the  errors  of  the  earth's  surface  radiance 
observations  under  cloudy  conditions  "’S  assumed  to  be  2  ergs/cm  -sec- 
strdn-cm  (Ref.  12).  Figure  A-l  shows  enat  for  the  most  ideal  conditions 
(i.e.,  clear  conditions,  availability  of  ground  pressure  data,  and  the 
smaller  ce  value)  the  r.m.s.  height  differences  are  no  more  than  20  ft 
even  at  the  10  mb  pressure;  while  for  the  worst  conditions  (i.e., 
cloudy  and  unconstrained  or  using  only  statistical  considerations)  the 
r.m.s.  height  differences  are  no  more  than  60  ft  even  at  10  mb.  It  is 
of  interest  to  note  that  for  clear  conditions,  the  geopotential  height 
estimates  show  a  significant  improvement  even  for  the  larger  og  value 
when  using  surface  measurements,  i.e.,  the  r.m.s.  error  is  reduced  to 
no  more  than  20  ft  at  the  200  mb  pressure  level. 


*The  smaller  value  represents  an  unrealistically  small  value,  and  is 
used  to  illustrate  a  limit  for  the  degree  of  accuracy  that  can  be 
achieved  from  satellite  radiance  observations.  The  high  value  is 
only  a  factor  of  two  smaller  from  the  value  that  leads  to  results 
which  do  not  meet  the  requirements  of  meteorology  (Ref.  15). 
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FIGURE  A-1  .  Altitude  Deviation  of  Atmospheric  Pressure  Levels  (in  milibars)  Derived 
from  Statistical  Properties  of  the  Geopotential  Height  and  Satellite¬ 
like  Radiance  Measurements.  The  altitude  deviation  is  given  for  clear 
and  cloudy  conditions,  and  the  RMS  difference  between  radiosonde  ob¬ 
served  and  radiance  calculated  values  (i.e.,  ce/  ergs/cm^-sec-stdn- 
cm-l).  The  smaller  value  of  10”^  for  clear  conditions  defines  a  limit 
for  the  accuracy  obtainable  from  satellite  radiance  observations,  while 
the  larger  value  (ae  =  1/4)  is  the  expected  typical  value  for  the  SIRS 
observations  (Ref.  12). 

4.  R.M.S.  Errors  in  the  Temperature  Distribution 


The  accurate  determination  on  a  global  and  real-time  scales  of 
the  temperature  distribution  in  the  troposphere  from  remote  radiance 
measurements  in  the  infrared  region  of  the  electromagnetic  s^ctrum 
has  already  been  shown  to  be  feasible  (Refs.  IS  to  27)  because  of  the 
following  five  factors : 

a.  The  existence  of  atmospheric  "windows'’  where  gaseous  absorp¬ 
tion  is  minimal.  These  windows  become  useful  for  the  deter¬ 
mination  of  surface  temperatures  and  cloud  conditions  over 
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the  local  regions  of  interest.  Two  windows  of  consequence  to 
the  subsequent  description  of  radiometers  are  one  in  the 
interval  from  3.5  to  4.1a,  and  another  one  in  the  8.0  to  12.5a 
range . 

b.  The  presence  in  the  atmosphere  of  gaseous  tracers  (i.e. ,  water 
vapor  and  carbon  dioxide)  that  are  well  mixed  with  the  air 
(with  a  rather  constant  mixing  ratio)  in  the  altitude  range 
from  the  earth's  surface  to  the  mesopause. 

c.  The  molecular  absorption  bands  (i.e.,  vibration-rotation  and 
rotational)  in  the  infrared  spectrum  of  the  above  gaseous 
tracers.  The  principal  infrared  absorber  in  the  atmosphere  is 
water  vapor  with  strong  vibration- rotation  bands  centered  near 
1.1a,  l-38y,,  1.87a,  2.7u,  and  6.3a  as  well  as  rotational  bands 
from  12  to  about  65a*  A  second  gaseous  absorber  of  importance 
is  carbon  dioxide  with  strong  vibration- rotation  bands 
centered  near  2.7a,  4.3a,  and  15a* 

d.  The  sampling  of  different  depths  in  the  atmosphere  by  using 
different  absorption  bands  and  different  zenith  angles,  i.e., 
as  illustrated  in  Fig.  A~2  for  the  15a  carbon  dioxide  interval 
for  zero  and  60  deg  zenith  angles.  This  figure  shows  that  as 
the  frequency  increases  from  695  cm-1  (i.e.,  14.4a)  to  747  cm"' 
(13.4a)  the  maximum  relative  atmospheric  contribution  to 
radiance  measurements  take  place  at  approximately  the  100  mb, 
400  mb,  700  mb,  and  the  earth's  surface.  A  measurement  made 
at  a  frequency  for  which  atmospheric  absorption  is  intense 
(e.g,,  695  cm-1)  records  radiation  only  from  the  upper  part  of 
the  atmosphere;  the  radiation  coming  from  lower  down  is 
strongly  attenuated  and  does  not  reach  satellite  altitudes. 
Thus  the  radiance  received  at  this  frequency  contains  informa¬ 
tion  about  the  upper  regions  of  the  atmosphere.  Conversely, 
measurements  made  at  a  frequency  of  small  absorption  (e.g., 

747  cm  ^)  detects  radiation  coming  near  or  at  the  surface  of 
the  earth.  The  data  in  Fig.  A-2  also  indicate  the  following: 
(1)  a  given  altitude  usually  contributes  radiation  to  more 
than  one  frequency;  because  of  this  inherent  redundancy,  it  is 
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reasonable  to  expect  that  there  is  a  limit  on  the  amount  of 
information  that  can  be  obtained  from  increases  in  tne  spectral 
resolution  (i.e.,  the  number  of  frequencies);  (2)  no  informa¬ 
tion  can  be  obtained  for  altitudes  that  do  not  make  a  reason¬ 
able  contribution  to  the  radiance  at  any  frequency;  and  (3)  the 
outgoing  radiance  from  a  particular  frequency  is  predominantly 
from  a  layer  approximately  10  km  (~5  nmi)  thick,  which  fixes 
the  vertical  resolution  of  the  deduced  temperature  distribution 
to  about  this  value. 
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RELATIVE  ATMOSPHERIC  CONTRIBUTION 


FIGURE  A-2.  Weighting  Functions  for  the  Chosen  15fi  CC>2  Spectral  Intervals 
for  Two  Different  Zenith  Angles  (Ref.  17) 


e.  The  scanning  of  radiometers  with  small  fields  of  view  (~10  m 
rad);  which  allows  an  increase  in  the  spatial  resolution  of 
the  measurements,  and  therefore  increases  the  probability  for 
radiance  measurements  through  broken  clouds.  The  unique 
advantages  of  deriving  atmospheric  temperature  profiles 
utilizing  data  from  a  High  Resolution  Infrared  R-  ation 
Sounder  (HIRS)  are  as  follows:  (1)  clouds  will  generally 
prohibit  accurate  profiles  to  be  obtained  for  the  troposphere 
except  in  the  central  region  of  storm  systems,  and  (2)  the 
greater  number  of  accurate  independent  measurements  over  a 
given  area  will  increase  the  accuracy  of  the  temperature  pro¬ 
file  for  that  area.  A  proposed  HIRS  subsystem  for  the  Nimbus 
F  satellite  would  utilize  six  different  spectral  intervals : 
four  in  the  15p  carbon  dioxide  band,  one  in  the  6.3p  water 
vapor  band,  and  one  in  the  lip,  window  region.  The  angular 
resolution  of  the  instrument  would  be  8  m  rad,  which  provides 
a  linear  resolution  on  the  earth  of  about  5  nmi  at  the  nadir 
from  a  spacecraft  altitude  of  600  nmi.  This  field  of  view  is 
sufficient  for  resolving  clear  cloud  interstices  as  needed  for 
deriving  accurate  tropospheric  temperature  profiles.  The 
radiometer  would  scan  ±  35  deg  of  nadir  angle,  to  each  side  of 
the  orbital  track  (Ref.  17).  The  areas  completely  free  of 
clouds  would  be  identified  from  a  statistical  frequency 

analysis  of  the  "window"  observations  obtained  over  a  rela- 
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tively  large  area  (i.e.,  150  nmi  ),  a  method  that  is  based  on 
the  following  factors:  (1)  the  major  variations  of  the  window 
radiance  measured  within  such  an  area  would  be  the  result  of 
variations  of  the  properties  of  any  clouds  within  the  field  of 
view  (i.e.,  amounts,  altitude,  opacities),  and  (2)  the  magni¬ 
tude  of  the  radiance  observations  through  the  clear  holes  in 
the  clouds  will  be  in  general  higher  than  those  through  the 
clouds;  as  illustrated  in  Fig.  A-3.  The  figure  illustrates  a 
frequency  histogram  of  the  Nimbus  II  HRIR  (High  Resolution 
Infrared)  radiance  observations  (in  units  of  brightness 


temperature'")  obtained  over  an  area  of  15C  rani 11  and  in  the 
immediate  vicinity  of  a  tropical  storm  (Ref.  17).  The  range 
in  the  brightness  temperature  includes  the  extremes  of  high 
density  cirrus  (210UK)  and  clear  (300°K)  conditions.  The 
observational  frequency  peak  at  296°K  is  associated  with  clear 
conditions,  which  is  recognizable  by  the  following  factors'. 

(i)  the  magnitude  and  repeatability  of  this  radiance  observa¬ 
tion,  and  (ii)  the  standard  deviation  of  the  distribution  on 
the  high  radiance  side  of  this  distinct  mode  is  close  to  the 
r.m.s.  error  of  measurement  (about  1.5  K  for  the  HRIR  obser¬ 
vations).  The  amount  of  dispersion  on  the  high  radiance  side 
of  the  frequency  mode  associated  with  the  highest  radiance 
values  can  be  used  to  discriminate  between  a  mode  associated 
with  totally  clear  observations  and  one  associated  with  either 
thin  overcast  or  low  dense  overcast  cloud  conditions.  This  is 
due  to  the  fact  chat  even  for  these  cloud  conditions,  the  dis¬ 
persion  will  tend  to  be  significantly  greater  than  that  ex¬ 
pected  from  surface  temperature  variations  (at  least  over 
water)  and  random  errors  of  measurements  (Ref.  17). 

The  reduction  of  radiance  data  from  a  given  column  (or  volume)  of 
the  atmosphere  to  a  temperature  profile  within  that  column  of  air  in¬ 
volves  the  solution  of  a  Fredholm  integral  equation  of  the  first  kind 
which  may  be  written  as  (Ref.  15) 

CO 

R(m)  K  (m,Z)  f  (2)  dZ 
0 

where  R  is  the  measured  radiance  ,  m  the  parameter  of  measurement 
(either  frequency  or  nadir  angle,  Fig.  A-2),  f  (Z)  the  atmospheric 
parameter  to  be  deduced  from  the  R  measurement  (e.g.  ,  temperature  or 
water  vapor  concentration),  S  the  altitude,  and  K  a  known  kernel  or 
weighing  function  (e.g.,  like  Tig.  A-2).  Since  the  inversion  procedure 

'•''The  brightness  temperature  is  defined  by  the  product  of  the  emissivity 
and  the  absolute  temperature. 


involved  requires  a  model  of  the  temperature  profile,  available  results 
indicate  the  following:  (1)  methods  that  utilize  statistical  functions 
to  model  temperature  profiles  in  terms  of  infrared  radiances  from  the 
earth's  atmosphere  will  yield  the  most  accurate  results  when  using  rela¬ 
tively  few  spectrally  independent  radiance  measurements  (e.g.,  as  those 
proposed  for  HIRS);  (2)  an  inversion  method  recently  developed  by 
Chahine  (Ref.  18),  enables  the  temperature  profile  to  be  derived  with 
comparable  accuracy  from  a  relatively  large  number  of  spectrally  inde¬ 
pendent  radiances  without  the  use  of  atmospheric  statistical  functions; 
(3)  the  r.m.s.  errors  in  the  temperature  profile  for  clear  and  cloudy 
conditions  are  as  indicated  in  Fig.  A-4  'Ref.  12).  These  results  are 
obtained  from  the  conditions  and  statistical  procedure  used  fc  ■  Fig. 

A-l.  This  figure  indicates  that  an  accuracy  of  ±  3°C  is  achievable  even 
at  the  high  altitudes;  (4)  typical  results  of  the  inversion  method  are 
shown  by  the  solid  lines  in  Figs.  A-S  and  A-5a  (Ref.  15),  which  are  also 
compared  with  the  dashed  liner  obtained  from  radiosonde  soundings.  The 
infrared  temperature  profiles  (solid  lines)  are  based  on  data  obtained 
from  a  high  altitude  balloon,  and  measurements  made  in  seven  narrow 
spectral  intervals  (six  in  the  15u  band  and  one  in  10u  window).  Figure 
A- 5  indicates  mean  deviations  of  2~K  from  the  radiosonde  soundings;  but 
the  errors  are  mainly  near  the  trcpopause,  which  is  systematically 
placed  too  low  because  of  instrument  problems.  Thus,  an  important  re¬ 
sole  which  holds  for  all  inversion  schemes  is  that  the  accuracy  of  the 
inverted  temperature  profile  depends  on  the  accuracy  of  the  measurements 
rather  than  in  the  number  of  frequencies  used  in  the  measurements. 

Figure  A-5a  shows  that  infrared  inversion  provides  no  information  below 
the  top  of  a  cloud  deck,  as  indicated  between  the  300  and  200  mb  pres¬ 
sure  levels;  and  (5)  the  vertical  resolution  and  accuracy  of  the 
temperature  profiles  cannot  be  improved  significantly  by  increasing  the 
spectral  resolution  (bandwidth),  i.e.,  Fig.  A-6  shows  that  the  r.m.s. 
errors  in  the  temperature  profile  deduced  from  radiance  measurements  in 
the  15p,  CO2  spectral  region  are  somewhat  insensitive-  for  half-bandwidths 
less  than  20  cm  1.  This  result  is  apparently  due  to  the  high  interlevel 
autocorrelation  of  temperature  that  exists  for  the  earth's  atmosphere; 
even  though  the  radiances  sensed  in  the  narrower  spectral  regions 


originate  from  slightly  smaller  vertical  layers  of  the  atmosphere,  the 
high  correlation  between  the  temperatures  within  adjacent  layers  pre¬ 
vent  improvements  in  the  vertical  resolution  of  the  temperature  pro¬ 
files  from  increases  in  the  spectral  resolution  (Ref.  17). 

The  above  considerations  have  emphasized  the  infrared  region  be¬ 
cause  such  has  been  the  case  during  the  first  decade  in  the  developments 
of  meteorological  satellites.  Although  the  use  of  microwave  techniques 
is  also  promis'ng  for  the  minimization  of  the  effect  of  clouds  on  the 
inferred  temperature  profiles,  such  techniques  become  of  increasing 
interest  for  the  potential  inference  of  the  wind  vector  over  the  oceans 
from  measurements  of  the  sea  surface  conditions. 


FIGURE  A-4.  Temperature  Deviation  of  Atmospheric  Pressure  Levels  (in  milibars)  Derived 
from  Statistical  Properties  of  the  Geopotential  Height  and  Satellite-like 
Radiance  Measurements.  The  temperature  deviation  is  given  for  clear  and 
cloudy  conditions,  and  the  RMS  difference  between  radiosonde  observed 
and  radiance  calculated  values  (See  Fig.  A-l).  The  clear  conditions  cor¬ 
respond  to  those  of  the  proposed  HIRS,  i.e.,  with  a  field  of  view  small 
enough  to  make  observations  through  broken  clouds  (Ref.  12). 
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FIGURE  A-5  ■  Compari.on  of  fhe  Temperature  Profile  Deduced  from  Remote 
Radiance  Measurements  (in  six  spectral  intervals  in  the  15p 
band)  under  clear  sky  conditions  from  a  high  altitude  balloon 
over  Texas  at  7:57  a.m.  (solid  lines)  and  the  temperature 
profile  corresponding  to  the  noon  Fort  Worth  radiosonde  values 
(dashed  lines)  (Ref.  15) 
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FIGURE  A-5 a  .  Comparison  of  the  Temperature  Profile  Deduced  from  Remote 
Radiance  Measurements  (in  six  spectral  intervals  in  the  I5ji 
band)  under  cloudy  conditions  from  a  high  altitude  balloon  over 
Texas  at  12:17  and  the  temperature  profile  corresponding  to  the 
noon  Forth  Worth  radiosonde  values  (Joshed  lines)  (Ref.  15) 
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FIGURE  A-6  .  RMS  Temperature  Error  Calculations  for  Different  Bandwidths 
and  Standard  Deviation  of  Temperature  of  Statistical  Sample. 
These  Results  are  Derived  from  the  Procedure  U^ed  in  Figs.  1 
and  4  (Ref.  17) 


5 .  R.H.S.  Errors  in  the  Water  Vapor  Distribution 

The  meteorological  parameters  that  can  be  inferred  iron  measure¬ 
ments  of  the  radiation  reflected  or  emitted  by  the  earth's  atmosphere 


-L.  ^  -> 


(e.g.,  E,  in  the  previous  equation)  are  as  follows:  temperature  pro¬ 
files  in  the  atmosphere  including  surface  temperatures  ,‘v  water  vapor 
distributions  including  surface  values,"  ozone  distributions,  and 
aerosols  as  well  as  cloud  characteristics.  Although  the  required 
transmission  or  kernel  functions  (K's)  for  carbon  dioxide  are  known 
with  barely  adequate  precision,  those  for  water  vapcr  need  further 
attention,  while  the  lack  of  adequate  data  for  ozone  is  a  major  ob¬ 
stacle  at  present  (Ref.  15). 


WATER  VAPORS,  w  (g/kg) 

I _ I _ 1 _ 1  1 _ I _ 1 _ I _ 1  1  I _ i  1 

200  220  240  260  280  300  320 

TEMPERATURE  (°K) 

FIGURE  A-/.  Comparison  of  Water  Vapor  (mixing  ratio)  and  Temperature  Profiles 
Deduced  from  Remote  Radiance  Measurements  (using  an  inter¬ 
ferometer  spectrometer)  under  Clear  Sky  Conditions  from  a  High- 
Altitude  Balloon  over  Texas  at  9:06  a.m.  (solid  lines)  and  Corre¬ 
sponding  Profiles  Derived  from  Nearby  Radiosonde  Ascents 
(circles  and  triangles)  (Ref.  15) 

For  clear  or  broken  cloud  conditions. 
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Figure  A- 7  illustrates  the  water  mixing  ratio  (in  g/kg)  and  the 
temperature  profiles  inferred  from  the  interferometer  spectrometer  data 
of  Chaney  (Ref.  28)  taken  from  a  high  altitude  balloon  when  no  clouds 
were  present.  The  circles  and  triangles  are  results  from  nearby  radio¬ 
sonde  ascents  while  the  solid  lines  give  the  inferred  values.  This 
figure  shows  good  agreement  of  the  inferred  water  vapor  profiles  with 
the  radiosonde  soundings.  Similar  results  have  been  found  for  condi¬ 
tions  when  there  were  broken  clouds  below  the  balloon  (Ref.  29).  Thus, 
it  appears  that  the  current  state-of-the-art  allows  the  use  of  infrared 
data  to  infer  the  total  amount  and  at  least  one  or  two  lapse  rates  of 
the  water  vapor  profile  in  the  atmosphere  (Ref.  15). 

6 .  R.M.S.  Errors  in  the  Wind  Field 

Since  the  influence  of  the  wi*  '  and  pressure  on  the  radiation 
field  is  very  small,  the  velocity  vector  of  the  air  motion  (i.e., 
magnitude  and  direction)  is  not  inferable  from  radiation  measurements 
in  the  infrared.  The  development  of  techniques  to  measure  the  wind 
field  from  satellites  ±s  still  in  a  very  early  stage  of  development. 

The  early  C-ARP  program  for  the  1970s  (circa  1973)  plans  to  measure  the 
wind  field  by  using  the  following  two  independent  but  complementary 
techniques  (Ref.  11): 

a.  About  900  constant-density  level  balloons,  used  as  tracers 
for  the  horizontal  air  motion  in  a  global  scale,  are  tracked 
by  radio  from  meteorological  satellites  operating  at  medium 
altitude  orbits.*  The  proposed  system  would  distribute  600 
balloons  on  a  nearly  world-wide  basis  at  the  150  to  200  mb 
pressure  level,  where  the  average  speed  is  high  (about  30  m/sec). 
This  system  would  also  make  use  of  300  balloons  at  the  850  mb 
level  in  the  tropics ,  where  the  average  wind  speed  in  the  lower 
troposphere  is  about  10  m/sec.  The  deployment  of  balloons  at 
additional  levels  in  the  mid-troposphere  is  constrained  by 

"The  balloons  may  also  contain  temperature  and  pressure  sensors  as 
well  as  radio  altimeters  to  measure  the  altitude.  Additional  balloons 
may  be  required  to  correct  for  a  nonuniform  distribution  of  the  bal¬ 
loons  over  the  globe  during  the  estimated  3C  to  40  days  lifetime  of 
each  balloon. 


present  balloon  limitations,  i.e.,  icing  at  night  or  within 
clouds  and  orographic  obstacles.  Thus,  this  constant  density 
level  balloon  system  will  yield  only  1  or  2  deg  of  freedom  of 
the  wind  in  the  vertical,  instead  of  the  required  four  or 
five  levels.  The  accuracy  of  the  wind  measurements  will  be 
well  within  the  desired  limit  of  ±  3  m/sec,  since  the  balloon 
location  can  be  determined  accurately  by  using  on  the  balloons 
either  two-way  transponders  or  stable  transmitters.  The  two- 
way  transponders  (EOLE  type)  allows  a  very  accurate  (1  to  2  km) 
location  of  the  balloon  as  well  as  individual  interrogation 
and  command  capability,  while  the  stable  transmitter  allows 
reasonable  location  accuracies  of  about  10  to  20  km.  Thus, 
observations  of  the  balloon  displacement  in  one  orbital  period 
(i.e.,  about  90  min)  of  the  meteorological  satellite  at  medium 
altitude  orbit  would  yield  average  winds  with  an  error  as  low 
as  0.15  to  0.30  m/sec. 

Using  small  scale  cloud  elements  as  tracers  for  the  horizontal 
wind  field,  and  tracking  such  cloud  elements  from  meteorological 
satellites  operating  at  geosynchronous  altitudes  (Ref.  30).  An 
advantage  of  geosynchronous  meteorological  satellites  is  that 
they  remain  fixed  relative  to  the  motion  of  the  weather  in  a 
global  scale.  Proposals  for  the  early  GARP  program  would 
utilize  two  kinds  of  high  resolution  instruments  for  the  obser¬ 
vation  of  the  cloud  elements:  (1)  an  optical  scanner  similar 
to  the  spin-scan  cameras  used  in  the  NASA  ATS-1  and  ATS-3 
satellites;  it  would  provide  daytime  cloud  tracking  with  a 
horizontal  resolution  of  2  to  3  km  with  an  accuracy  for  hori¬ 
zontal  displacement  of  about  1  km;  (2)  a  scanning  infrared 
radiometer  operating  in  the  8  to  12(j.  window  region;  it  would 
provide  the  altitude  or  cloud  tops  with  an  accuracy  of  ±  0.5  km 
(assuming  a  temperature  profile  and  uniform  and  opaque  cloud), 
a  horizontal  resolution  of  8  to  12  km  with  an  accuracy  for 
horizontal  displacement  of  about  5  km.  This  IR  radiometer 
would  provide  pictures  with  a  spatial  resolution  range  (at 
nadir)  of  about  8  x  8  to  16  x  16  km  (depending  on  the  optics) 
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of  the  disc  of  the  earth,  with  a  full  picture  available  every 
30  min.  This  system  would  yield  daytime  wind  data  from  cloud 
elements  that  are  larger  than  8  to  12  km,  but  which  have  a 
detectable  structure  of  smaller  scale.  Since  the  nighttime 
wind  data  will  not  be  as  good,  this  system  is  limited  to  day¬ 
time  and  cloudy  conditions . 

7.  Chronological  Development  of  U.S.  Meteorological  Satellites 

A  perspective  of  the  current  and  future  potential  state-of-the-art 
in  the  remote  measurement  of  atmospheric  parameters  from  satellites  can 
best  be  established  from  a  brief  review  of  the  chronological  develop¬ 
ment  of  U.S.  meteorological  satellites  and  their  sensor  subsystems. 

The  RDT&E  phase*  of  U.S.  meteorological  satellites  has  consisted 
of  three  steps  TIROS,  NIMBUS,  and  ATS.  Ten  TIROS  (Television  Infra¬ 
red  Observational  Satellites)  were  launched  into  medium  altitude  orbits 
between  April  1,  1960,  and  July  2,  1965,  ay  NASA  (Refs.  31,  32,  33). 
This  TIROS  series  was  then  followed  by  the  TIROS  Operational  Satellites 
(TOS),  which  began  to  carry  out  global  cloud  photography  operationally 
and  routinely  for  the  Environmental  Sciences  Services  Administration 
(ESSA),  The  six  NIMBUS  satellites  were  launched  into  medium  altitude 
orbits  between  1964  and  1971,  The  NIMBUS  series  was  set  up  to  develop 
sensors  and  technology  for  a  "second  generation"  operational  system. 

The  Advanced  Technology  Satellites  (ATS)  were  launched  into  geo¬ 
synchronous  altitude  orbit,  A  chronological  description  of  the  main 
TIROS  and  NIMBUS  series  launched  to  date  is  shown  in  Table  1,  while 
Fig.  A-8  includes  also  the  ATS  launchings.  Each  of  these  three  RDT&E 
series  as  well  as  the  ESSA  operational  phase  is  described  below. 


*The  RDT&E  phase  of  a  new  program  involves  the  research,  development, 
test,  and  evaluation  of  a  prototype  system  for  the  subsequent  opera¬ 
tional  phase  of  the  program. 


TABLE 

A-l.  TIROS  AND 

NIMBUS  RDT&E 

METEOROLOGICAL 

SATELLITES 

PRE LAUNCH 

IN  ORBIT  LAUNCH 

ORBIT 

ORBIT 

DESIGNATION 

DESIGNATION  DATE 

INCLINATION  ALTITUDE  SUESYSTE  o 

REMARKS 

(deg) 

(nmi ) 

TIROS  A-l 

TIROS  I 

APR  1 

•60 

48 

400 

2  V’  vid 

( 1W+1N  lens) 

First  meteor¬ 
ological 
satellite 

TIROS  A- 2 

TIROS  II 

NOV  23 
*60 

48 

400 

2  V  vid 

(1W+1N  lens) 

1  MRIR 

1  wide-angle 
cone  rad 

First  IR 

TIROS  A-3 

TIROS  III 

JULY  12 
•61 

48 

400 

2  V'  (2W  lens) 

1  MRIR 

1  wide  angle  cone 
rad 

1  omnidir.  rad 

TIROS  D 

TIROS  IV 

FEB  8 
■62 

49 

400 

2  V  vid 

( 1W+1M  lens) 

1  MRIR 

1  wide  angle 
cone  rad 

1  omnidir.  rad. 

TIROS  E 

TIROS  V 

JUNE  19 

•62 

58 

400 

2  V  vid 

(1W+1M  lens) 

TIROS  F 

TIROS  VI 

SEP  18 
'62 

58 

400 

2  V  vid 

(1W+1M  lens) 

TIROS  G 

TIROS  VII 

JUNE  19 
•63 

58 

400 

2  V  vid 
(2  W  lens) 

1  MRIR 

1  Omnidir.  rad 

TIROS  H 

TIROS  VIII 

DEC  21 
■63 

58 

400 

1  APT 

1  V  vid 
(W  lens) 

First  APT 

NIMBUS  A 

NIMBUS  I 

AUG  28 
■64 

Sun-sync. 

600 

1  Trimetroyen 
AVCS 

1  APT 

1  HRIR  (56LR) 

First  HRIR 

TIROS  I 

TIROS  IX 

JAN  22 
'  6S 

Sun-sync. 

400/ 

1600 

2  V  vid 
( 2  W  lens ) 

1st  worldwide 
coverage  - -Is 
cartwheel  ‘C 

TIROS  (OT-1) 

TIROS  X 

JULY  2 

•6S 

Sun-sync. 

400 

2  V  vid 
(2  W  lens) 

Axial  mode 

NIMBUS  C 

NIMBUS  11 

MAY  IS 
'66 

Sun -sync. 

600 

1  Trimetrogen 

1  MRIR 

1  HRIR  (S6LR) 

AVCS 

*  The  TIROS  satellite  is  an  18-sided  right  cylinder,  22V  high,  42"  in  diameter, 
and  weight  about  300  lbs  (Ref.  34).  The  NIMBUS  satellite  is  described  in  Fig.  9 
and  weighs  about  700  lbs. 
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FIGURE  A-8  .  RDT&E  and  Operational  U.5.  Meteorological  Satellites  Placed 
in  Orbit  during  Nearly  a  Decade.  Note  that  ten  TIROS  (T)  and 
one  NIMBUS  (N)  RDT&E  Satellites  were  utilized  during  the  first 
half  of  the  decade,  while  eight  ESSA  (E)  account  for  most  of 
the  successful  launchings  so  far  in  the  second  half  of  the  decade 

The  TIROS  series  was  devoted  to  the  identification  and  tracking 
on  a  world-wide  and  synoptic  basis  of  weather  phenomena  such  as  cloud 
cover  mapping,  storm  tracking  (e.g.,  frontal  systems,  hurricanes,  etc.), 
radiation  balance,  cloud  top  heights,  surface  temperatures,  tropo¬ 
spheric  water  vapor,  and  stratospheric  temperature  and  circulations. 
Observations  consisted  initially  of  photographing  cloud  patterns  with 
television  cameras  in  daytime,  which  were  supplemented  later  by  infra¬ 
red  radiometers  to  make  cloud  observations  at  night.  Tju  first  satel¬ 
lite  in  this  series  (TIROS  I)  was  placed  in  a  circular  orbit  at  an 
altitude  of  400  nmi  and  inclination  of  48  deg  to  the  equator.  The 
sensor  subsystem  consisted  of  two  television  (1/2  in.  vidicon)  cameras; 
one  with  a  wide  angle  lens  of  104  deg,  and  the  other  with  a  narrow 


angle  lens  of  12.7  deg.  The  second  meteorological  satellite  (TIROS  II) 
launched  in  the  same  year  added  to  the  camera  system  the  first  infrared 
sensors,*  which  consisted  of  a  Medium  Resolution  Infrared  Radiometer 
(MRIR)  and  one  wide-angle  cone  radiometer.  The  MRIR  provided  a  5  deg 
field  of  view,  a  ground  resolution  of  30  nmi  from  400  nmi  orbit  alti¬ 
tude,  and  scanning  channels  in  the  following  spectral  bands:  0.55  to 
0.75y,,  0.2  to  6p,,  6  to  6.5u,  and  8  to  12u,.  The  wide-angle  cone  had  a 
50  deg  field  of  view,  a  ground  resolution  of  350  nmi,  and  was  used  to 
measure  the  equivalent  local  blackbody  temperature  and  albedo  of  the 
earth.  The  third  meteorological  satellite  (TIROS  III)  was  launched  in 
the  subsequent  year  (July  1961) ,  and  its  sensor  subsystems  added  an 
omnidirectional  radiometer  to  measure  long  and  shortwave  radiation 
over  the  entire  earth’s  disc.  It  also  modified  the  camera  system  to 
use  both  cameras  with  wide  angle  lens.  TIROS  IV,  V,  and  VI  were 
launched  during  1962  and  used  the  same  type  of  sensor  subsystems  as 
in  the  previous  satellites,  except  that  the  two-camera  system  was 
modified  to  use  one  wide-angle  lens  (104  deg)  and  one  medium-angle 
lens  (78  deg).  Likewise,  TIROS  VII  and  VIII  were  launched  in  1963 
with  the  same  type  of  sensor  subsystems ,  except  that  the  two-camera 
subsystem  of  TIROS  VIII  utilized  one  camera  with  wide  lens ,  and  one 
Automatic  Picture  Transmission  (APT)  1  in.  vidicon  camera.  The  TIROS 
series  concluded  with  the  launchings  of  TIROS  IX  and  X  during  1965. 

The  NIMBUS  series  of  six  RDT&E  meteorological  satellites  started 
v/ith  the  launching  of  NIMBUS  I  (August  28,  1964)  and  NIMBUS  II  (May  15, 
1966)  into  a  sun- synchronous  orbit  at  an  altitude  of  600  nmi.**  The 
remainder  of  the  series  consists  of  NIMBUS  III  (launched  on  April  14, 


*The  satellite  Explorer  VII  had  carried  earlier  a  rather  simple  array 
of  hemispheric,  omnidirectional  sensors  to  measure  the  radiation 
balance  of  the  earth. 

‘*A  sun  synchronous  orbit  at  this  altitude  allows  at  least  two  obser¬ 
vations  every  24  hr  of  practically  all  points  on  rhe  globe  and  at 
the  same  local  time.  The  orbit  plane  is  inclined  to  the  equator 
98.7  deg  and  precesses  around  the  center  of  the  earth  at  a  rate  that 
is  synchronous  with  the  revolution  of  the  earth  around  the  sun. 


1969  and  described  in  Fig.  9),*  NIMBUS  IV  (to  be  launched  in  1970), 

NIMBUS  V  and  NIMBUS  VI  (both  being  planned  for  1971).  The  basic  NIMBUS 
satellite  is  a  system  defined  by  the  following  characteristics:  (1)  a 
structure  designed  to  accommodate  a  maximum  number  of  experiments,  (2) 
about  200  watts  of  electrical  power  for  experiments,  (3)  pointing  of  the 
instruments  at  the  earth  at  all  times  with  an  accuracy  of  about  1  deg, 

(4)  a  stable  and  moderate  thermal  environment  required  by  many  instruments, 

Q 

and  (5)  a  data  subsystem  that  can  acquire  and  store  in  the  order  of  10  to 
12 

10  bits  per  orbit,  and  transmit  these  data  to  the  ground  while  passing 
over  Alaska.  The  missions  of  NIMBUS  I  and  NIMBUS  II  included  the 
following  objectives:  (a)  demonstrate  improvements  in  making  high 
resolution  daytime  cloud  observations,  (b)  demonstrate  the  feasibility 
of  high  resolution  nighttime  cloud  mapping,  (c)  transmit  both  day  and 
nighttime  cloud  images  to  local  receivers,  and  (c)  map  emitted  telluric 
radiation  as  well  as  reflected  solar  radiation  in  various  spectral 
bands.  The  three  subsystems  of  NIMBUS  I  consisted  of  the  following:  a 
set  of  three  television  cameras  and  a  tape  recorder  called  the  Advanced 
Vidicon  Camera  System  (AVCS)  for  daytime  cloud  photography,  an  addi¬ 
tional  television  camera  to  transmit  daytime  cloud  pictures  via  the  APT 
to  simple  receivers,  and  the  first  scanning  High  Resolution  Infrared 
Radiometer  (KRIR)  with  Stored  and  Local  Readout  (S&LR)  capability.  The 
HRIR  provided  a  0.46  deg  field  of  view,  and  a  ground  resolution  of  5  nmi 
in  the  spectral  band  from  3.5  to  4  •  l-u,  •  The  subsystems  on  NIMBUS  II  con¬ 
sisted  of  the  AVCS,  the  HRIR  (SS-LR),  and  a  Medium  Resolution  Infrared 
Radiometer  (MRIR)  with  fixed  channels  to  map  and  image  thermal  radiation 
emitted  by  (a)  atmospheric  water  vapor  or  cirrus  clouds  (6.4  to  6.9^), 

(b)  carbon  dioxide  (14  to  16^),  (c)  cloud  tops  or  the  earth's  atmosphere 
(5  to  3C\j. ) ,  and  (e)  reflected  solar  radiation  (0.2  to  4.0U).  The  HRIR 
demonstrated  that  cloud  formations  can  be  observed  globally  at  night 
with  a  spatial  resolution  of  about  4  nmi,  and  that  their  heights  are 
inferable  from  the  measured  cloud  top  temperatures.  Thus,  the  HRIR 
produced  at  night  three  dimensional  cloud  maps  in  a  global  scale^,  which 

*The  NIMBUS  III  satellite  was  placed  in  an  orbit  defined  by  a  perigee 
of  582  nmi,  apogee  of  613  nmi,  and  inclination  of  39.9  deg. 
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are  unattainable  with  television  cameras.  The  10.5  to  11. Sy,  channel 
of  the  MRIR  produced  cloud  cover  and  neight  maps  for  day  and  nighttime 
conditions  with  a  considerable  lower  spatial  resolution  (~25  nmi). 

Based  on  these  results,  the  ESSA's  TOS  will  incorporate  the  use  of  a 
high  resolution  llu  scanning  radiometer  for  day  and  night  cloud  mapping. 
Other  significant  results  from  NIMBUS  I  and  NIMBUS  II  are  as  follows : 

(a)  determination  of  temperatures  of  the  earth’s  surface  with  suffi¬ 
cient  accuracy  to  allow  for  the  tracking  of  ocean  currents  (as  mapped 
by  the  HRIR  in  cloud  free  areas);  (b)  improved  insight  of  the  global 
distribution  of  the  net  radiative  energy  flux  through  the  upper 
boundary  of  the  atmosphere  (Ref.  35);  and  (c)  better  understanding  of 
the  morphology  of  the  circulation  in  the  lower  stratosphere  (Ref.  36). 
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FIGURE  A-9  .  The  NIMBUS  III  Spacecraft 


Besides  the  uwo  main  series  of  RDT&E  meteorological  satellites 
(TIROS  and  NIMBUS)  in  medium  altitude  orbit,  a  third  RDT&E  phase  in¬ 
volves  three  NASA  Application  Technology  Satellites  (ATS)  for  operation 
in  geosynchronous  orbit.*  Meteorological  satellites  in  synchronous 
altitude  will  allow  nearly  continuous  observations  of  the  same  areas  of 
the  earth,  since  the  meteorological  application  of  the  ATS  observations 
in  the  form  of  time  lapse  movies  have  already  shown  the  weather  in 
motion  over  a  given  geographical  area.**  While  the  infrared  signals 
become  very  low  at  geosynchronous  altitude ,  the  ATS  system  allows  the 
observation  of  the  development  of  local,  short-time  scale  meteorological 
phenomena  (e.g.,  growth  of  severe  storms,  the  tracking  of  small  cloud 
elements  for  wind  measure  events,  etc.).  The  first  ATS  system  (ATS-I) 
was  launched  into  geosynchronous  orbit  on  December  7,  1966,  the  second 
one  (ATS-III )  on  November  5,  1967,  and  a  third  one  (ATS-E)  is  planned 
for  launch  in  August  1969  (Ref.  31).  The  characteristics  and  perform¬ 
ance  of  the  meteorological  sensor  subsystem  on  ATS-I  and  III  are 
summarized  in  Table  2  (Ref.  37).  The  sensor  subsystem  of  ATS-I  con¬ 
sisted  of  a  spin  scan  cloud  camera  in  brack  and  white  for  the  Spin  Scan 
Cloud  Cover  Experiment  (SSCEE).  This  instrument  is  really  a  scanning 
telescopic  photometer  with  a  field  of  view  of  0.1  milliraaian,  which 
utilizes  the  spin  of  the  satellite  to  scan  the  scene  from  west  to  east, 
and  which  tilts  the  telescope  slowly  through  a  small  angxe  to  scan  the 
scene  from  north  co  south.  The  transmission  of  the  pictures  was  accom¬ 
plished  by  satellite  facsimile  or  Weather  Facsimile  (WEFAX).  The  ATS- 
III  utilized  a  spin  scan  cloud  camera  using  three  in-line  pinholes  to 
generate  three  color  components:  red,  blue,  and  green.  Fiber  optics 
are  used  to  carry  the  color  information  to  three  separate  photo¬ 
multipliers.  The  first  tube  is  equipped  with  a  green  filter,  the 


*The  ATS  systems  are  nonmeteorological  satellites  that  have  allowed 
limited  experiments  to  detect  mainly  cloud  motions. 

**The  time  lapse  technique  speeds  up  the  motion  of  the  weather ,  i.e., 
one  can  see  what  took  place  in  IB  hr  of  a  day  (time  for  the  sun  to  go 
from  horizon  to  horizon)  in  about  3  or  4  sec  (Ref.  30). 

+Two  other  ATS  launch  attempts  took  place  on  April  6,  1967  (ATS-II), 
and  August  1C,  1963  (ATS-IV). 


second  with  a  red  filter,  and  the  third  with  a  blue  filter.  The 
signals  thus  generated  are  transmitted  to  the  ground.  The  information 
is  sent  in  groups  of  four  signals.  The  first  signal  carries  informa¬ 
tion  concerning  the  intensity  of  green,  the  next  signal  represents  red, 
and  the  third  one  blue.  The  last  signal  in  the  sequence  is  a  fixed 
reference  signal  to  help  balance  the  colors  when  all  the  signals  are 
integrated  into  a  complete  picture  at  the  receiving  stations.  The  ATS- 
III  also  included  one  Image  Dissector  Camera  System  (IDCS)  and  used 
WEFAX  and  Omega  Position  Location  Equipment  (OPLE)  for  the  transmission 
and  collection  of  data  from  remote  sensors.  The  use  of  WEFAX  through 
the  ATS-I  and  ATS-III  satellites  has  demonstrated  the  practicality  of 
transmitting  data  collected  and  analyzed  at  a  central  location  to  users 
in  remote  locations . 


TABLE  A-2.  SPIN-SCAN  SYSTEM  CHARACTERISTICS  AND  PERFORMANCE 


Spacecraft 

ATS-I 

ATS-III 

Size,  in. 

7  x  11  x  10 

7  x  11  x  12 

Weight,  lb 

<16 

23.5 

Primary  Power 

-24  vdc,  0.6  A 

-24  vdc,  0.9  A 

Resolution 

0.1  m  rad  (2  mi) 

0.1  m  rad  (2  mi) 

Dynamic  Range 

>400:1 

>400:1 

Scan  Lines 

2000 

2400 

Earth  Coverage 

~50  deg  N  to  50  deg  S 

Full  earth 

Scan  Rate 

100  lines /min 

100  lines /min 

Picture  Time 

20  min 

24  min 

Retrace  Time 

2  min 

2  min 

IFOV3  Dwell  Time 

9.6p  sec 

9.6p  sec 

Video  Bandwidth 

100  KHz 

100  KHz/channel 

Total  Bandwidth 

100  KHz 

4  MHz 

Spectral  Response 

0.48-0. 63p 

Blue  0.38-0.48u, 
Green  0.48-0. 58^ 
Red  0.55-0.63^ 

aInstantaneous  Field  of  View. 


ISO 


The  operation  phase  to  date  of  meteorological  satellites  consists 
of  eight  ESSA  satellites  launched  into  medium  altitude  orbit  between 
February  3,  1960,  and  December  15,  1968.  The  ESSA  satellites  have 
utilized  sensor  subsystems  developed  in  the  earlier  NASA  TIROS  series ; 
except  for  the  use  of  a  Flat  Plate  Radiometer  (FPR),  which  is  a  varia¬ 
tion  of  the  omnidirectional  radiometer  for  measurements  of  longwave 
emission  of  radiation.  The  ESSA  meteorological  satellites  have  been 
launched  at  the  approximate  average  rate  of  three  satellites  each  year 
(Table  A-3). 


The  subsystems  on  the  meteorological  satellites  flown  to  date  were 
designed  to  map  cloud  cover  by  using  visual  and  radiometric  sensors , 
and  to  map  telluric  and  reflected  solar  radiation  in  the  various  spec¬ 
tral  bands.  Most  of  the  sensors  being  developed  for  future  use  repre¬ 
sent  improvements  of  present  sensors  to  provide  better  spatial  and 
spectral  resolution.  Some  sensors  will  provide  observations  in  several 
spectral  bands  so  as  to  allow  inference  of  temperature  distributions 
and  water  vapor  profiles  through  the  troposphere.  A  precise  perspec¬ 
tive  of  the  current  projections  for  the  development  of  meteorological 
satellite  subsystems  is  best  derived  from  a  brief  review  of  the 
experiments  planned  for  the  research  satellites  NIMBUS  III  (or  NIMBUS 
B-2 ,  launched  on  April  14,  1969)  and  NIMBUS  IV  (or  NIMBUS  D,  to  be 
launched  in  March  1970).* 

The  experiments  on  NIMBUS  I  and  II  were  devoted  exclusively  to 
improving  the  synoptic  mapping  of  meteorological  phenomena  or  to  making 
atmospheric  measurements  of  general  scientific  interest  from  medium 
altitude  orbit  (600  iuni).  The  experiments  on  NIMBUS  III  and  IV  put 
emphasis  on  the  development  of  subsystems  for  the  determination  from 
medium  altitude  orbrt  (600  rimi)  of  initial  and  boundary  conditions  for 
the  numerical  integration  of  the  equations  of  motion,  i.e.,  the 

'•'•'The  subsystems  for  the  remaining  two  NIMBUS  satellites  (l.e.,  N.MBUS 
E  and  F)  are  as  yet  unspecified. 
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measurement  of  temperature  and  water  vapor  profiles  through  the  tropo¬ 
sphere,  the  measurement  of  wind  fields  by  tracking  free  floating 
balloons  distributed  globally  at  various  height  levels,  and  the  real¬ 
time  recovery  of  surface  data  from  buoys.* 

The  NIMBUS  III  carried  out  seven  experiments,**  three  of  which 
put  emphasis  on  the  requirements  for  the  numerical  integrations  of  the 
equations  of  motion,  one  is  for  scientific  investigations,  and  three 
continue  the  mapping  observations  of  NIMBUS  I  and  II.  The  first  group 
of  three  experiments  consists  of  two  spectrometer  experiments  to  yield 
data  for  the  determination  of  temperature  and  water  vapor  profiles  in 
the  troposphere  and  one  experiment  capable  of  tracking  balloons  and 
relaying  measurements  from  automatic  stations.  These  experiments  are 
as  follows: 

a.  One  spectrometer  is  a  Micheison  Interferometer  called  Infra- 
Red  Interferometer  Spectrometer  (IRIS),  which  operates  between 
6  and  20p  at  a  relative  spectral  distribution  of  1:200.  In¬ 
cluded  within  this  spectral  interval  are  the  water  vapor 
absorption  band  centered  at  6.3p,  the  9.6p  ozone  band,  and  the 
15p  carbon  dioxide  band. 

b.  The  second  spectrometer  is  a  modified  Fastie-Ebert  grating 
spectrometer  called  Satellite  Infra-Red  Spectrometer  (SIRS). 
Radiant  energy  is  detected  in  seven  spectral  intervals  of  the 
15p  carbon  dioxide  band.  The  spectral  intervals  are  5  cm-1 
wide.  An  eighth  channel  senses  radiation  in  the  atmospheric 
window  centered  at  11. Ip,. 

c.  The  data  relay  and  balloon  tracking  experiment  consists  of  a 
satellite  borne  transmitter,  receiver,  and  computer.  This 
experiment  is  called  Interrogation,  Recording,  and  location 
System  (IKLS).  By  using  communication  with  a  given  automatic 

*An  operational  system  of  meteorological  satellites  operating  at 
medium  altitude  orbit  will  require  several  satellites  to  yield 
global  coverage  in  real-time. 

**These  experiments  were  described  in  August  1968  (Ref.  33)  and,  there¬ 
fore,  exclude  possible  subsequent  changes  prior  to  the  actual  flight. 
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ground  station,  the  IRLS  can  determine  the  location  of  a 
station  within  about  1  runi,  and  can  also  interrogate  a  set  of 
sensors  (such  as  a  thermistor)  contained  in  the  station.  This 
IRLS  experiment  is  an  engineering  test  of  a  system  involving 
observations  of  six  balloons  over  the  United  States  and  utiliz¬ 
ing  20  ground  stations. 

The  scientific  experiment  consists  of  a  Monitor  of  Ultraviolet 
Solar  Energy  (MUSE),  which  will  expand  the  capability  of  NIMBUS  to 
perform  scientific  investigations  relevent  to  meteorology.  MUSE  will 
measure  solar  radiation  in  five  spectral  intervals  (each  100  A  wide) 
ranging  from  1200  A  to  2600  A.  The  three  experiments  that  will  continue 
the  mapping  observations  of  NIMBUS  I  and  II  are  as  follows:  The  HRER 
(S&LR)  with  day  and  night  capability,  the  MRIR,  and  the  Image  Dissector 
Camera  System  (IDCS). 

The  NIMBUS  IV  plans  to  carry  out  nine  experiments,  five  of  which 
emphasize  the  requirements  for  mathematical  prediction  models,  two  will 
make  basic  scientific  investigations,  and  the  remaining  two  are  for  the 
mapping  of  horizontal  fields.  The  first  group  of  five  experiments  con¬ 
sists  of  four  experiments  to  yield  temperature  (all  four  of  them)  and 
water  vapor  (three  of  them)  profiles  in  the  troposphere,  and  one 
experiment  for  the  tracking  of  balloons .  These  experiments  are  as 
follows : 

a.  The  IRIS  experiment  (as  in  NIMBUS  III)  with  its  spectral  range 
extended  to  the  range  between  8  and  40^,. 

b.  The  SIRS  experiment  (as  in  NIMBUS  III)  with  the  addition  of 
six  spectral  intervals,  mostly  in  the  rotational  water  vapor 
bands  between  18  and  30^. 

c.  A  Filter  Wedge  Spectrometer  ( EWS )  operating  between  1.2  and 

6 . 4p, . 

d.  A  Selective  Chopper  Radiometer  (SCR)  operating  in  the  15yj, 
carbon  dioxide  band,  which  is  expected  to  achieve  a  very  high 
spectral  resolution  by  filtering  the  radiation  through  carbon 
dioxide  absorption  cells. 


e.  The  IRLS  experiment  (as  in  NIMBUS  III)  will  be  expanded  to  be 
capable  of  interrogating  hundreds  of  stations  instead  of  20  on 
NIMBUS  III.  Also  a  balloon  experiment  will  be  performed  to 
measure  the  wind  field  over  a  large  area  (such  as  the  tropics) 
within  at  least  one  height  level. 

The  two  scientific  investigations  will  use  the  MUSE  experiment  (as 
in  NIMBUS  III)  and  a  Backscatter  Ultraviolet  (BUV)  spectrometer  to 
measure  the  intensity  of  solar  radiation  reflected  by  the  atmosphere  in 
14  intervals  (each  10  A  wide)  over  the  spectral  range  from  2500  to 
3400  A.  The  BUV  experiment  will  yield  the  global  vertical  distribution 
of  ozone  between  7  and  25  nmi,  approximately.  The  two  experiments  for 
the  mapping  of  the  horizontal  fields  consists  of  the  IDCS  experiment 
(as  in  NIMBUS  III)  to  provide  television  pictures,  and  a  scanning 
Temperature  and  Humidity  Infrared  fcidiometer  (TffIR)  operating  in  the 
6.3p,  water  vapor  and  the  11^  window  bands  with  a  spatial  resolution  of 
the  HRIR.  The  THIR  experiment  will  yield  measurements  of  the  horizontal 
water  vapor  and  cloud  fields. 

Besides  the  main  projected  RDT&E  activities  of  the  NIMBUS  II  and 
NIMBUS  III  satellites,  it  is  also  of  interest  to  review  the  development 
of  three  radiometers  for  the  new  ESSA  Improved  TIROS  Operational  System 
(ITOS)  as  well  as  a  proposed  very  high  resolution  radiometer  for 
synchronous  satellites."'  The  three  radiometers  for  the  ITOS  are  as 
follows : 

a.  ITOS  high  resolution  radiometer  with  two  channels  in  the  10.5 
to  12. 5p.  and  0.52  to  0.73^  intervals.  The  field  of  view  of 
these  channels  are  0.30  deg  and  0.16  deg,  respectively,  and 
yield  corresponding  ground  resolutions  from  the  nominal  orbit 

*The  projected  activities  for  1969  of  the  operational  ESSA  ser’ies  in¬ 
volve  two  additional  launchings  of  the  current  TOS  (TIROS  Operational 
System),  and  a  subsequent  shift  to  the  new  ITOS.  Four  projected  ITOS 
launchirgs  (TIROS  M,  ITOS  A,  B,  and  C)  will  utilize  the  following 
eight  subsystems:  two  AVCS,  two  APT,  two  HRIR  (S&LR),  one  FPR,  and 
one  SPM  (Solar  Proton  Monitor).  'The  SPM  will  measure  the  flux  of 
protons  in  the  10.30  and  60  M-ev  range,  and  electrons  in  the  100  ara 
750  E-ev  range. 


altitude  of  the  ITOS  (775  nmi)  of  3.7  and  2  nmi.  Because  of 
this  high  spatial  resolution,  the  quality  of  the  photo-images 
should  approach  that  of  television  cameras.  Some  other  advan¬ 
tages  of  this  radiomecer  are:  (1)  its  capability  to  measure 
surface  temperatures  ir.  the  20.5  to  12. 5p,  interval  equally 
well  during  day  and  nighttime  conditions,  (b)  the  identifica¬ 
tion  of  cloud  free  areas  in  the  daytime  by  means  of  the  coor¬ 
dinated  visible  Channel  measurements,  and  (c)  its  adaptation 
to  the  APT  system  to  yield  a  Direct  Readout  of  InfraRed  (DRIR) 
system.  The  DRIR  system  will  operate  during  the  entire  orbit, 
reading  out  infrared  imagery  to  APT  stations  all  over  the 
world  in  the  daytime  (when  the  infrared  data  will  be  inter¬ 
leaved  with  the  standard  APT  pictures)  as  well  as  in  the  night¬ 
time  . 

ITOS  Vertical  Temperature  Profile  Radiometer  (VTPR)  with  eight 
channels,  six  of  which  will  be  in  the  lSo  carbon  dioxide  band, 
one  in  the  rotation  water  vapor  band,  and  one  in  the  11^ 
window.  The  spectral  resolution  for  all  channels  will  be 
about  5  cm  x.  The  VTPR  is  intended  solely  for  the  probing  of 
the  temperature  profile  from  operational  meteorological  satel¬ 
lites,  and  it  takes  advantage  of  the  latest  advances  in  tech¬ 
nology,  i.o.,  the  weight  (10  lb)  and  power  (2  watts)  are 
significantly  smaller  than  chose  of  earlier  instruments 
designed  for  vertical  temperature  probing  (e.g.,  SIRS:  92  lb 
and  20  watts;  IRIS:  28  Id  and  12  watts). 

ITOS  Very  High  Resolution  Radiometer  (VHRR)  with  two  channels 
in  the  .10.5  to  12. Sp  and  0.6  to  C.7g  intervals.  The  field  of 
view  c>f  each  channel  is  0.034  deg,  which  yields  a  ground  reso¬ 
lution  of  0.4  nmi  from  the  nominal  altitude  orbit.  Thus,  the 
ground  resolution  of  the  VHRR  represents  an  order  of  magnitude 
improvement  over  that  of  the  NIMBUS  HRIR  and  ITOS  High  Resolu¬ 
tion  Radiometer,  and  nearly  two  orders  of  magnitude  improvement 
over  that  of  the  TIROS  and  NIMBUS  medium  resolution  radio¬ 
meters.  -As  such,  the  resolution  of  the  VHRR  is  comparable  to 
that:  of  television  systems  currently  flown  on  meteorological 
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satellites.  This  level  of  resolution  is  achieved  by  using  a 
mercury-cadmium-telluride  detector  cooled  radiatively  to  80°K. 
This  high  resolution  will  allow  the  study  of  smaller  scale 
phenomena  such  as  the  three-dimensional  structure  of  cloud 
patterns  and  frontal  systems.  Since  the  bit  rate  of  the  in¬ 
strument  exceeds  the  capacity  of  any  available  tape  recorder 
system,  it  may  be  necessary  to  readout  the  data  directly  via 
a  relay  geosynchronous  satellite  (e.g.,  see  Ref.  38). 

The  proposed  Very  High  Resolution  Radiometer  for  Geosynchronous 
Altitude  (not  yet  approved)  with  one  channel  in  the  10.5  to  12. 5p, 
interval  has  a  detector  similar  to  the  VHRR.  The  field  of  view  of  this 
instrument  is  0.023  deg,  which  yields  a  ground  resolution  of  8  nmi  from 
geosynchronous  altitude.  The  advantages  or  a  very  high  resolution 
radiometer  over  a  camera  that  senses  visible  radiation  are  as  follows: 
(1)  the  radiometer  system  will  take  "full  earth"  pictures  continuously 
and  independently  of  the  position  of  the  terminator.  The  familiar 
"full  earth"  pictures  made  with  the  ATS  cameras  can  only  be  taken  during 
a  small  fraction  of  the  day,  since  pictures  taken  during  a  24  hr  period 
show  varying  smaller  portions  of  the  earth  (depending  on  the  location  of 
the  terminator);  and  (2)  unlike  the  camera  system,  the  radiometer  has  a 
potential  for  temperature  measurements;  which  can  in  turn  be  used  to 
infer  cloud  top  heights,  sea  surface  temperature  variations  (in  cloud- 
free  regions),  ana  the  altitude  of  winds  deduced  from  time  lapse  move¬ 
ments  of  cloud  patterns . 

Table  A-4  shows  a  summary  of  the  main  characteristics  of  all  the 
radiometric  instruments  that  have  been  developed  during  the  last  decade 
for  the  remote  probing  of  the  troposphere  and  stratosphere  from  meteoro¬ 
logical  satellites.*  The  first  13  instruments  in  this  table  have  been 
either  flown  or  approved  for  the  RDTS-E  (i.e.,  TIROS  and  NIMBUS  HAS  A 
series)  and  operational  (TOS  and  ITOS ,  ESSA  series)  programs,  while  the 
last  two  have  been  proposed  for  future  RDTSE  meteorological  satellites 

"The  mechanical  details  of  the  first  14  instruments  are  described  in 
Ref.  32  (and  the  references  thereto),  while  the  proposed  HIRE  design, 
is  treated  in  some  detail  in  Ref.  17, 
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operating  either  at  geosynchronous  (ATS)  or  medium  (NIMBUS  F)  altitude 
orbits .  The  characteristics  of  these  IS  instruments  indicate  the 
following  two  trends  during  the  first  decade  of  remote  probing  of  the 
lower  atmosphere  from  meteorological  satellites:  (1)  an  exclusive 
emphasis  on  the  development  of  radiometric  instruments  in  the  infrared, 
and  (2)  an  evolution  process  toward  improved  spectral  and  spatial  reso¬ 
lutions.  The  improved  spectral  resolutions  yield  temperature  and  water 
vapor  vertical  profiles  through  the  troposphere,  including  the  surface 
values  for  cloudless  conditions;  while  the  simultaneous  use  of  high 
spectral  and  spatial  resolutions  (i.e.,  HIRS)  will  allow  the  determina¬ 
tion  of  the  complete  temperature  and  water  vapor  vertical  profiles  over 
a  large  area  under  broken  cloud  conditions . 

9.  The  Global  Atmospheric  Research  and  World  Weather  Watch  Program 

A  brief  outline  of  the  international  Global  Atmospheric  Research 
and  World  Weather  Watch  Programs  (GARP  and  WWW ,  respectively)  is  neces¬ 
sary  in  order  to  gain  a  complete  perspective  of  the  planned  activities 
for  the  next  decade  concerning  efforts  for  the  prediction  of  the  global 
variations  of  atmospheric  parameters  mainly  within  the  troposphere. 

The  GARP  program  is  essentially  the  international  RDTc-E  phase  for 
an  eventual  operational  WWW  in  a  global  scale.  The  GARP  program  was 
proposed  in  1967  jointly  by  relevant  committees  of  the  International 
Council  of  Scientific  Unions  (ICSU)  and  the  World  Meteorological 
Organization  (WMO )  of  the  United  Nations.  'The  main  objective  of  GARP 
is  to  provide  research  and  development  efforts  directed  toward  an 
experiment  that  will  measure  the  large  scale  motions  of  the  entire  lower 
atmosphere  for  a  limited  period  of  time  (Refs.  11,  39,  40).  This  plan 
set  1972  as  the  year  for  intensive  observations  on  a  global  scale,*  and 
outlines  problems  such  as  relevant  scientific  areas  (e.g. ,  energy 
transfer-motions  of  the  tropical  troposphere,  interactions  at  the 
earth’s  surface,  the  dynamics  of  the  boundary  layer,  etc.)  as  well  as 
new  methods  of  observation  and  data  handling.  Since  the  ICSU  has  no 

*Current  estimates  indicate  now  that  1975-1976  is  a  more  realistic  date 
for  such  effort. 


resources  to  fund  such  a  program,  the  tempo  of  activities  in  GARP  will 
be  set  by  voluntary  efforts  within  the  world  community.  Besides  the 
theoretical  efforts  concerning  the  conceptual  design  of  an  RDTS-E  global 
system  (Ref.  11),  some  major  experiments  compatible  'with  GARP  are  as 
follows:  (a)  the  RDT&E  meteorological  satellites  such  as  the  NIMBUS 
and  ATS  series  described  above;  (b)  the  U.S.  GHOST  (Global  Horizontal 
Sounding  Technique)  project,  which  has  demonstrated  the  capability  for 
constant-density  balloons  to  fly  in  the  stratosphere  for  average 
periods  of  six  months;*  and  (c)  the  French  EOLE  balloon  project,  which 
will  utilize  a  satellite  for  the  interrogation  of  balloons. 

The  World  Weather  Watch  program  was  approved  in  1967  by  the  World 
Meteorological  Organisation  (a  specialized  agency  of  the  United  Nations). 
The  main  objective  of  the  WWW  program  is  data  collection  on  a  global 
scale  to  be  processed  in  world  and  regional  weather  centers.  The  data 
collection  would  utilize  conventional  observing  techniques ,  meteoro¬ 
logical  satellites ,  coordinated  meteorological  rocket  soundings ,  tele¬ 
metering  ocean  buoys  and  automatic  land  stations,  telemetering  free- 
flying  constant  density  balloons,  radi"  "ndes  launched  from  merchant 
ships,  and  other  observations  from  commercial  aircraft.  The  data  so 
obtained  would  be  handled  by  an  augmented  meteorological  communications 
network,  including  geosynchronous  satellites  with  sensor  and  relay  sub¬ 
systems.  Again,  since  the  WHO  cannot  fund  such  programs,  the  tempo  of 
activities  in  the  WWW  program  will  be  set  by  the  voluntary  efforts  of 
the  individual  Member  States .  The  significant  activities  within  the 
WWW  program  are:  (a)  a  four  year  suggested  program  for  1968-1371  adopted 
by  the  Fifth  World  Meteorological  Congress  in  1967  for  a  Global  Observing 
System.  This  program  includes  the  use  of  meteorological  satellites,  the 
establishment  of  a  large  number  of  new  upper  air  stations,  expanded 
programs  of  observations  at  existing  stations ,  an  increased  number  of 
selected  ships  for  making  surface  observations  over  the  oceans ,  and  be¬ 
tween  five  and  ten  additional  ocean  weather  stations  (Ref.  41);  (b)  a 
plan  for  the  addition  of  Regional  Meteorological  Centers  to  the  World 


*The  life  of  che  GHOST  balloons  in  the  troposphere  is  coo  short  to  be 
useful. 


Centers  in  Melbourne,  Moscow,  and  Washington;  and  (c)  the  role  of  the 
WWW  through  the  following  four  main  aspects:  (1)  serving  as  a  catalytic 
and  coordinating  body  in  such  matters  as  standardization  of  methods  and 
hours  of  global  observation,  international  coding  procedures,  merchant 
ship  voluntary  observing  system,  telecommunications  systems  for  the 
international  exchange  of  meteorological  data,  etc.;  (2)  providing 
direct  assistance  to  needy  countries  regarding  meteorological  equipment 
and  training;  (3)  providing  the  machinery  for  continuous  planning;  and 
(4)  providing  a  link  at  the  governmental  level  with  the  United  Nations 
on  such  programs. 

10.  Conceptual  Techniques  for  Potential  Measurements  from 

Meteorological  Satellites 

A  potential  technique  for  the  global  measurement  (in  near  real-time) 
of  the  boundary  condition  involving  the  wind  field  over  the  oceans  con¬ 
sists  in  the  measuring  of  the  thermal  radiation  emitted  by  the  oceans  in 
the  microwave  region  of  the  electromagnetic  spectrum.  This  use  of 
passive  microwave  measurements  would  also  yield  the  temperature  and 
water  vapor  vertical  profiles  over  oceans  under  cloudy  conditions,  which 
would  be  important  to  supplement  the  measurements  by  high  resolution 
radiometers  in  the  infrared  such  as  HIRS  (Table  A-4).  For  these  reasons, 
a  brief  review  is  given  below  of  the  following  two  topics:  (1)  a  pre¬ 
liminary  microwave  experiment  involving  remote  measurements  from  an  air¬ 
craft  at  altitudes  up  to  about  6  nmi  of  the  radiation  emitted  by  the 
ocean  at  a  wavelength  of  1.55  cm  (Ref.  42);  and  (2)  the  feasibility  of 
inferring  temperature  and  water  vapor  vertical  profiles  over  the  oceans 
under  cloudy  conditions  by  passive  microwave  measurements  from  meteoro¬ 
logical  satellites  at  medium  altitude  orbits." 


,VA  contrast  must  now  be  made  between  sensor  subsystems  for  the  global- 
passive  (remote)  measurements  from  satellites  and  those  for  local- 
active  (remote)  measurements  from  the  ground,  fc.g.,  (1)  ultrasensitive 
radars  to  measure  winds  (Ref.  43),  (2)  doppler  radar  to  measure  atmo¬ 
spheric  motions  in  virtually  all  scales  (Ref.  43),  (3)  LIDAR  (Light 
Detection  and  Ranging)  to  map  aerosol  distributions  (Refs.  44,45),  etc. 


The  microwave  experiment  from  the  aircraft  is  planned  to  demon¬ 
strate  the  feasibility  of  such  measurement  from  the  NIMBUS  E  meteoro¬ 
logical  RDT&E  satellite  (Refs.  15,  42).  It  utilizes  an  electrically 
scanning  radiometer  that  consists  of  the  following  three  elements. 

A  phased-array  antenna  accepts  radiation  in  a  beam  about  2.8  deg  wide. 
The  beam  electrically  scans  in  a  direction  perpendicular  to  the  air¬ 
craft’s  flight  direction  and  through  an  angle  ±  50  deg  from  nadir.  A 
receiver  measures  the  detected  radiation  with  an  accuracy  correspond¬ 
ing  to  about  ±  1°C  of  the  equivalent  brightness  temperature. 

Associated  calibration  and  recording  equipment  are  used.  Since  the 
spectral  intensity  of  the  radiation  emitted  at  1.55  cm  is  propor¬ 
tional  to  the  brightness  temperature  (i.e.,  the  product  of  the 
emissivity  and  temperature  of  the  water  surface),  the  radiometer  was 
calibrated  with  reference  to  the  radiation  source  with  a  known  bright¬ 
ness  temperature.  In  addition  to  the  microwave  radiometer,  the 
experiment  utilized  an  MRIR  radiometer  in  the  10.5  to  11. 5p,  spectral 
range  to  measure  the  equivalent  blackbody  temperature  of  the  sea  sur¬ 
face.  Wind  speeds  over  the  ocean  surface  were  obtained  from  aircraft- 
drift  measurements  and  estimates  of  the  wave  heights.  Some  results  of 
pertinent  interest  from  these  experiments  are  shown  in  Fig.  A-10,  where 
the  brightness  temperature  is  plotted  as  a  function  of  the  angle  from 
nadir  for  smooth  (i.e.,  negligible  wind  of  less  than  3  m/sec)  and  rough 
(i.e.,  a  wind  speed  of  about  15  m/sec)  sea  surface  conditions.  This 
figure  also  shows  similar  theoretical  results  as  derived  by  Stogryon 
(Ref.  46).  The  experimental  results  show  a  higher  brightness  tempera¬ 
ture  (by  about  20  deg)  for  the  rough  sea  condition  regardless  of  the 
angle  from  nadir  in  the  vertical  plane  of  observation,  a  result  that 
stems  from  the  increase  in  the  emissivity  from  (a)  an  increase  in  the 
roughness  scale  (i.e.,  foam,  spray,  white  caps,  and  heavy  swells),  and 
(b)  a  decrease  in  the  water  temperature  of  the  rough  relative  to  the 
smooth  surface.  The  theoretical  lines  in  Fig.  A-10  indicate  good  agree¬ 
ment  with  the  experimental  data  only  for  the  smooth  sea  condition,  a 
result  that  is  a  consequence  of  the  rather  poor  understanding  of  the 
relationship  between  the  roughness  scale  and  emissivity.  This  effect 
of  sea  surface  roughness  on  the  brightness  temperature  of  microwave 
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emission  was  also  verified  by  the  KRIR  radiometer,  which  measured 
higher  radiances  for  the  rough  sea  surface  relative  to  those  for  the 
smooth  sea  surface.  These  results  then  indicate  that  measurements  of 
microwave  radiation  emitted  at  1.55  cm  can  be  interpreted  in  terms  of 
sea  surface  roughness.  Since  the  roughness  scale  of  the  sea  surface 
is  intimately  related  to  the  prevailing  winds ,  these  results  indicate 
that  a  worldwide  mapping  of  the  sea  state  from  meteorological  satel¬ 
lites  can  yield  the  global  wind  field  over  the  oceans.* 

The  inversion  of  temperature  and  water  vapor  profiles  from 
measurements  of  emitted  radiation  in  the  microwave  region  of  the 
electromagnetic  spectrum  depends  on  two  factors:  (1)  the  brightness 
temperature  of  the  earth’s  surface,  and  (2)  the  transparency  of  the 
atmosphere  (Ref.  15).  Atmospheric  emissions  in  the  microwave  spectrum 
(e.g.,  the  weak  resonance  of  water  vapor  at  1.35  cm  wavelength,  and 
the  strong  complex  of  oxygen  resonances  at  0.5  cm)  can  bo  measured  with 
microwave  radiometers ,  which  yield  the  received  radiance  in  terms  of 
the  brightness  temperature  (e.g.,  Fig.  A-10).  A  microwave  radiometer 
looking  down  at  the  surface  of  the  earth  can,  therefore,  distinguish 
between  atmospheric  and  surface  emissions  only  when  the  brightness 
temperature  of  the  surface  is  low  (i.e.,  over  the  water,  where  the 
emissivity  is  low)  relative  to  the  atmospheric  temperature.  The  trans¬ 
parence  of  the  atmosphere  allows  then  inversion  of  temperature  profiles 
in  the  nearly  opaque  regions  of  the  spectrum  (in  a  manner  exactly 
analogous  to  that  in  the  infrared  spectral  region),  and  of  water  vapor 
profiles  in  the  nearly  transparent  regions  in  some  circumstances.  Since 
microwave  measurements  from  meteorological  satellites  could  be  used  to 
infer  temperature  profiles  below  clouds,  microwave  techniques  offer  the 
potential  to  yield  global  wind  fields  over  the  oceans  and  complete 
temperature  profiles  over  the  major  water  fraction  of  the  earth’s 
surface . 

•''Scanning  microwave  measurements  from  meteorological  satellites,  at 
medium  altitude  orbit  would  be  made  at  wavelengths  longer  than  3  cm  to 
minimize  the  effect  of  cloud  interference,  provided,  of  course,  that 
the  effect  of  roughness  on  the  brightness  temperature  (i.e.,  Fig. 

A-10)  persists  at  these  wavelengths . 
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11.  Concluding  Remarks 


A  review  of  the  RDT&E  activities  over  nearly  the  last  decade  con¬ 
cerning  satellite  measurements  of  atmospheric  parameters  indicates  that, 
in  general,  the  use  of  meteorological  satellites  operating  in  medium 
altitude  (e.g.,  400  to  750  nmi)  and  geosynchronous  (13,300  rani)  orbits 
can  provide  remote  measurements  with  adequate  accuracy  in  a  global  and 
near  real-time  scales  of  the  atmospheric  parameters  that  are  required 
for  the  numerical,  long-range  prediction  (maximum  of  about  3  to  5  weeks) 
of  the  global  atmospheric  circulation  within  mainly  the  troposphere  (up 
to  about  7  to  10  nmi). 

Some  specific  results  concerning  the  satellite  measurement  of 
atmospheric  parameters  are  as  follows : 

1.  Radiance  measurements  in  the  infrared  region  of  the  electro¬ 
magnetic  spectrum  from  meteorological  satellites  operating  in 
medium  altitude  orbit  can  determine  the  geopotential  height 
and  vertical  temperature  profile  to  within  60  meters  and  3'K, 
respectively,  in  an  altitude  range  up  to  about  10  mb,  or 
approximately  16.5  nmi  (Figs.  A-l,  A-4).  'The  advent  of  high- 
spatial-resolution  scanning  radiometers  (e.g.,  HIRS)  would 
provide  even  higher  accuracies  in  the  geopotential  height  and 
vertical  temperature  distribution  (Fig.  A-4)  over  most  of  the 
globe.  These  results  are  based  on  (a)  the  use  of  updated 
statistical  data  for  the  geopotential  height  and  temperature 
in  the  troposphere  and  lower  stratosphere,  and  (b)  the  solu¬ 
tion  of  the  inversion  problem  to  deduce  the  temperature  pro¬ 
file  from  radiance  measurements  in  different  spectral  bands 
(Fig.  A-2)  through  the  use  of  an  iterative  method  that  utilizes 
statistical  atmospheric  temperature  profiles. 

2.  The  main  emphasis  in  the  RDT&E  of  sensor  subsystems  for 
meteorological  satellites  over  nearly  the  last  decade  and  for 
the  near  future  (Tables  A-l,  A-2,  A-4)  has  been  toward  the 
development  of  high-spatial-resolution  scanning  radiometers; 
i.e.,  sensor  subsystems  that  would  yield  accurate  determina¬ 
tion  of  the  geopotential  height  and  vertical  temperature 


distributions  to  sea  level  altitude  over  most  of  the  globe 
(Fig.  A-4). 

3.  Radiance  measurements  in  the  infrared  region  of  the  electro¬ 
magnetic  spectrum  from  meteorological  satellites  operating  in 
medium  altitude  orbit  can  also  yield  the  vertical  temperature 
profile  through  the  troposphere  (Fig.  A-7). 

4.  Potential  techniques  for  measuring  the  velocity  vector  of  the 
air  motion  or  wind  within  troposphere  in  a  global  and  near 
real-time  scales  are  as  follows:  (1)  microwave  measurements 
from  meteorological  satellites  operating  at  medium  altitude 
orbit  of  the  thermal  radiation  emitted  by  the  oceans  to  pro- 

"  vide  data,  of  brightness  temperature  as  a  function  of  the  sea 
surface  condition  (Fig.  A-10),  which  is  controlled  by  the 
magnitude  of  the  wind  over  the  sea  surface;  (2)  the  use  of 
meteorological  satellites  in  geosynchronous  altitude  (e.g., 

ATS- I  and  ATS-III  satellites  provided  with  a  spin-scan  camera 
subsystem  as  indicated  in  Table  A-2)  to  provide  data  on  the 
motion  of  cloud  elements  over  fixed  geographical  locations. 

A  correlation  of  the  motion  of  certain  cloud  elements  with 
wind  measurements  would  yield  the  wind  vector  in  cloudy  areas; 
and  (3)  the  tracking  of  a  worldwide  network  of  constant  density 
balloons  from  meteorological  satellites  would  yield  the  wind 
vector  over  geographical  areas  and  atmospheric  levels  without 
clouds . 
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Techniques  which  have  potential  uae  for  e*pl  c  H term nail; o^t> r£eas«e£e«  The 
by  observations  from  satellites  are  evaluated  <°r  h  =  •  :  “  \  and  it5 

discussion  describes  the  instrumental  technique  f  Md  i  a  tier. ,  de¬ 
constituents  by  atmospheric  million  sateliite-to-sateilite  microwave 

celeration  of  low  altitude  sate-.t  ’  utgil  -.aser  sionals  and  wind  by  the  motion  of 
transR’-i'ii  ions ,  backscatte:  by  du3u  w  P  -  ~~  7  c,At elites  and  the  displacement  of 
clouds  determined  by  time-sequenced  pict a.  -  of  0p?l-aticr.,  uncertainty 

balloons  cracked  in  position  by  satellites^  I* ^development  0f  each  of  the  techniques 
of  measurement,  coverage  in  vo.Ji,  an-  -ta.us  ..  ue/...r 

considered  are  detail*.  J .  .  .  „.,CMrn-«r*-  of  a^Tiosoheric  pararrecers,  and  im- 

in  addition  ro  the  methods  of  *tation  of  che  aqua- 

plicit  method  of  determining  winds  a,  a-  a-titude,  y  .  derived  inputs  of 

tions  of  continuity  of  momentum,  4"*  “  ‘!n  application  of  the 

densities,  temperature  and  heating,  IS  sue t^S‘; -  ^r  than  a  dav  ahead  is 

method  for  purposes  of  numerical  prediction  of  weathe^  .0  tropospheric  simula- 

dus.-ribed.  Special  problems  in  extrapolating  the  exper.enc-  ^  J  data  are 

tier,  tor  the  determination  of  winds  of  the  upper  a  _  p..  •  d  Esri„dtos  are  giver, 

potentially  available  by  measurement  s  non'  ,  .e  -/ ,  I>pQcted  capacity  of  next- 
indicating  that  the  upper  atmosphere  simulation  Is  k-t.ar.  t.ic  -  P-  P 

ge.-.'rat  ion  computers  SU'  ii  as  ILLlnC  b  . 
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